## Abstract

We study the mechanical behavior of two-dimensional cellular tissues by formulating the continuum limit of discrete vertex models based on an energy that penalizes departures from a target area A0 and a target perimeter P0 for the component cells of the tissue. As the dimensionless target shape index s0=(P0/A0) is varied, we find a transition from a soft elastic regime for a compatible target perimeter and area to a stiffer nonlinear elastic regime frustrated by geometric incompatibility. We show that the ground state in the soft regime has a family of degenerate solutions associated with zero modes for the target area and perimeter. The onset of geometric incompatibility at a critical s0c lifts this degeneracy. The resultant energy gap leads to a nonlinear elastic response distinct from that obtained in classical elasticity models. We draw an analogy between cellular tissues and anelastic deformations in solids.

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

Article number | 268105 |

Journal | Physical Review Letters |

Volume | 120 |

Issue number | 26 |

DOIs | |

State | Published - Jun 29 2018 |

## ASJC Scopus subject areas

- Physics and Astronomy(all)

## Fingerprint Dive into the research topics of 'Geometric Frustration and Solid-Solid Transitions in Model 2D Tissue'. Together they form a unique fingerprint.

## Cite this

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

*Physical Review Letters*,

*120*(26), [268105]. https://doi.org/10.1103/PhysRevLett.120.268105

**Geometric Frustration and Solid-Solid Transitions in Model 2D Tissue.** / Moshe, Michael; Bowick, Mark J.; Marchetti, M. Cristina.

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

*Physical Review Letters*, vol. 120, no. 26, 268105. https://doi.org/10.1103/PhysRevLett.120.268105

**Geometric Frustration and Solid-Solid Transitions in Model 2D Tissue**. In: Physical Review Letters. 2018 ; Vol. 120, No. 26.

}

TY - JOUR

T1 - Geometric Frustration and Solid-Solid Transitions in Model 2D Tissue

AU - Moshe, Michael

AU - Bowick, Mark J.

AU - Marchetti, M. Cristina

N1 - Funding Information: Moshe Michael 1,2 ,* Bowick Mark J. 2,3 Marchetti M. Cristina 2 Department of Physics, 1 Harvard University , Cambridge, Massachusetts 02138, USA Department of Physics and Soft and Living Matter Program, 2 Syracuse University , Syracuse, New York 13244, USA Kavli Institute for Theoretical Physics, 3 University of California , Santa Barbara, California 93106, USA * Corresponding author. mmoshe@fas.harvard.edu 29 June 2018 29 June 2018 120 26 268105 5 March 2018 31 August 2017 © 2018 American Physical Society 2018 American Physical Society We study the mechanical behavior of two-dimensional cellular tissues by formulating the continuum limit of discrete vertex models based on an energy that penalizes departures from a target area A 0 and a target perimeter P 0 for the component cells of the tissue. As the dimensionless target shape index s 0 = ( P 0 / A 0 ) is varied, we find a transition from a soft elastic regime for a compatible target perimeter and area to a stiffer nonlinear elastic regime frustrated by geometric incompatibility. We show that the ground state in the soft regime has a family of degenerate solutions associated with zero modes for the target area and perimeter. The onset of geometric incompatibility at a critical s 0 c lifts this degeneracy. The resultant energy gap leads to a nonlinear elastic response distinct from that obtained in classical elasticity models. We draw an analogy between cellular tissues and anelastic deformations in solids. National Science Foundation 10.13039/100000001 DMR-1435794 DMR-1609208 DMR-1435999 Kavli Institute for Theoretical Physics, University of California, Santa Barbara 10.13039/100005956 PHY-1125915 United States - India Educational Foundation 10.13039/100006222 Simons Foundation 10.13039/100000893 Living tissues are far-from-equilibrium materials capable of spontaneously undergoing large-scale remodeling and adapting their mechanical behavior in response to internal and external cues. The experimental observation of glassy dynamics in epithelia [1–3] has motivated interest in quantifying the relation between rheological and structural properties of tissue. The aim is to provide a framework for organizing biological data by describing tissue as a material, with mechanical behavior tuned by effective parameters that provide a coarse-grained description of both intra- and intercellular interactions. Significant progress has been made using a class of models that describe a confluent cell monolayer, in which there are no gaps or overlaps between cells, as a tiling of space. Each cell is a polygon [vertex models (VMs)] [4–8] or a Voronoi area (Voronoi models) [9,10] , with the polygons’ vertices or the Voronoi centers taken as the degrees of freedom. For uniform cell edge tensions, both models are based on a tissue energy that penalizes deviations of the cell area A and perimeter P from prescribed target values A 0 and P 0 . Numerical solutions of vertex and Voronoi models with disordered polygonal configurations reveal rich behavior. Most interesting is the prediction of a rigidity transition tuned by cellular shape, as quantified by the dimensionless target cell shape index s 0 = P 0 / A 0 , which is in turn controlled by cell-cell adhesion and cortex contractility [9–13] . Vertex and Voronoi models have also been used successfully to quantify the relation between cellular shape and rearrangements and tissue mechanics in the developing Drosophila embryo [5,6] . A powerful tool for describing the mechanical properties of matter is continuum elasticity. While continuum models of epithelia have been developed and used to describe biological processes, such as wound healing and morphogenesis [14–16] , the development of a continuum theory that incorporates the rich behavior of the tissue energy used in discrete vertex and Voronoi models remains an open challenge [17] . Here, we tackle this challenge by considering a regular polygonal tiling and develop a geometric formulation of VM elasticity similar to that used to describe disordered solids or materials that exhibit nonuniform differential growth, such as plant leaves [18–22] . Our model is a coarse-grained version of a vertex model, albeit with uniform cell edge tensions, because it allows for shape changes through variations in the position of the vertices of the polygons. At a critical s 0 ( s 0 c ), corresponding to the isoperimetric quotient [23] , we find a transition between a soft and a stiff solid. In the soft solid ( s 0 > s 0 c ), the target area and perimeter of individual cells are simultaneously satisfiable (compatible), whereas in the stiff solid ( s 0 < s 0 c ) they are not (incompatible). This geometric frustration is associated with a sharp rise in the effective stiffness of the tissue and the onset of residual stresses. The critical value of the target shape index depends on the geometry of the unit cell, with s 0 c = 8 3 ≃ 3.722 for hexagons. For s 0 < s 0 c no hexagonal polygon exists, supporting a geometric origin for the transition. Earlier work has shown that the hexagonal ground state of the VMs is linearly unstable for s 0 > s 0 c , where it is replaced by a soft network of irregular polygons [8] . Here, we show that the ground state of this soft solid is a family of degenerate area-preserving and perimeter-preserving states with a band of zero modes. The onset of geometric incompatibility below s 0 c lifts this degeneracy and results in an energy gap, leading to a finite residual stress or prestress, as commonly seen in living tissues [24] . Although our starting point is a continuum elasticity seemingly quadratic in the strain, which would suggest a linear response at all values of imposed deformations, the force extension curves are nonlinear due to the degeneracy of the target configurations. The incompatible tissue shows strain stiffening, which is observed ubiquitously in living cells [25,26] . The solid-solid (SS) stiffening transition obtained here for regular polygonal tilings is distinct from the solid-liquid (SL) rigidity transition predicted earlier for disordered tilings [11,12] . The latter is driven by the vanishing of the energy barriers for T 1 transformations, which are forbidden in our model. Our work suggests that the SL rigidity transition may be facilitated by a SS transition, where the effective Young’s modulus of the soft solid phase becomes very small, thus easing T 1 transformations. This is supported by the recent suggestion that the SL rigidity transition in a disordered Voronoi model of 3D cellular agglomerates may also be associated with underlying geometric incompatibility [27] . The geometric formulation of elasticity used here highlights the underlying geometric nature of SS and SL transitions in vertex and Voronoi models and allows for the analytical calculation of stress-strain curves for regular lattices. The formalism can also be extended to incorporate disordered structures. Finally, our model is appropriate for rigid, solidlike tissue with tight adherent junctions between the cells and does not incorporate the “giant number fluctuations” observed in the fluid regime of epithelial tissue [28–30] or cell division. The latter is likely to fluidify the tissue at long times [13] and could be included in our continuum model through a target area metric that grows in time. The study of these effects is left for future work. In VMs, cells are modeled as polygons that can independently adjust their area A i and perimeter P i according to the energy [4,5,7,8] E T = 1 2 ∑ i [ κ A ( δ A i / A 0 ) 2 A 0 + κ P ( δ P i / P 0 ) 2 P 0 ] , (1) with δ A i = A i - A 0 and δ P i = P i - P 0 . The stiffnesses κ A and κ P have dimensions of energy per unit area and perimeter, respectively. The first term in Eq. (1) arises from bulk elasticity as well as the ability of cells to adjust their area by changing their height. The second term describes the interplay of cortical tension and cell-cell adhesion. It has been shown numerically that the disordered vertex model exhibits a solid-liquid transition at s 0 * ≈ 3.81 [11,12,31] , where the energy barrier for bond-flipping T 1 transitions that remodel the local cell neighborhood vanishes. This prediction for the disordered case has been validated in bronchial cells [32] . A generalization that includes cell motility has yielded a surface of solid-liquid transitions tuned by cell speed and the persistence time of single-cell dynamics [9] . Here, we do not include T 1 transitions, or any other topological excitations such as cell divisions. We show that even this simple limit exhibits unusual elastic behavior. Geometric formulation of tissue energy.— In the geometric approach to linear elasticity, a thin planar sheet is described as a surface equipped with a metric, a 2 × 2 symmetric tensor that locally specifies the distance between points on the surface [33–35] . Simple elastic solids are characterized by a global reference configuration, or target metric g ¯ , that is stress-free in the absence of external constraints or loads [36] . The strain tensor for a deformed state with actual metric g is defined as u = 1 2 ( g - g ¯ ) . The elastic energy of an isotropic Hookean solid spanning a region Ω is then given by E HS = 1 2 ∫ Ω A α β γ δ u α β u γ δ d S g ¯ , (2) where A α β γ δ = λ g ¯ α β g ¯ γ δ + μ ( g ¯ α γ g ¯ β δ + g ¯ α δ g ¯ β γ ) is the elastic tensor, λ and μ are Lamé constants, and d S g ¯ = det g ¯ d 2 x . While the geometric formulation of elasticity may appear to be unnecessarily formal, it is useful when describing materials laden with defects and solids with nonuniform differential growth that do not possess stress-free target configurations [37] . In such cases the material is prestressed , meaning that there is a residual stress even without an external load. As a result there is no global stress-free target or reference configuration, and displacement fields are consequently ill defined. The definition of strain as a deviation of the actual metric from a target one is, however, still valid and reflects the existence of local stress-free configurations [35] . This formulation naturally captures the so-called incompatible elasticity of such materials [19,22,38,39] . Our goal is to obtain a coarse-grained form of the tissue energy of Eq. (1) , and to express it in terms of a local measure of strain. Since area and perimeter can be tuned independently, different target metric tensors g ¯ A and g ¯ P are needed to characterize A 0 and P 0 : this requires two measures of strain. A single metric tensor g , however, characterizes the deformed state. Defining u A , P = 1 2 ( g - g ¯ A , P ) , the simplest energy functional quadratic in strains is E = E P + E A , with E A , P = ∫ Ω 1 2 A A , P α β γ δ u α β A , P u γ δ A , P d S g ¯ A , P , (3) where A x α β γ δ = λ x g ¯ x α β g ¯ x γ δ + μ x ( g ¯ x α γ g ¯ x β δ + g ¯ x α δ g ¯ x β γ ) , (4) with x ∈ { A , P } . A derivation of Eq. (3) from the discrete equation (1) is carried out in the Supplemental Material (SM) [40] following the procedure used in Ref. [41] for flexible membranes. It confirms Eq. (3) with λ P = μ P = 3 κ P P 0 / 8 and λ A = 2 μ A = κ A A 0 / 3 . In the remainder of this Letter, we consider for simplicity a spatially uniform system. An example of a calculation for a tissue with a nonuniform shape parameter is shown in a Mathematica notebook attached as SM [40] . Ground states.— Now comes an important subtlety—a given target area or perimeter can be realized by a family of target metrics, rather than just a single one. For a lattice of quadrilaterals, for instance, two families of metrics G A = { g ¯ A ( ε A ) } and G P = { g ¯ P ( ε P ) } corresponding to a given target area and perimeter are, respectively, g ¯ A = α A [ ε A 0 0 ε A - 1 ] , g ¯ P = [ ( α P + ε P ) 2 0 0 ( α P - ε P ) 2 ] , (5) with P 0 = 2 ( g ¯ 11 P + g ¯ 22 P ) = 4 α P and A 0 = det g ¯ A = α A , both in units of the lattice constant. This is illustrated in Fig. 1 for a tissue of quadrilateral cells with κ A = 0 . Figures 1(a) and 1(b) display two “undeformed” configurations with P = P 0 that can be exchanged, with no work, by applying a uniaxial strain in the x direction. This is commonly called a zero mode. Figures 1(c) and 1(d) show deformed configurations with P ≠ P 0 . It is evident that Fig. 1(c) deviates only slightly from Fig. 1(a) but is highly deformed compared to Fig. 1(b) . Similarly, Fig. 1(d) is close to Fig. 1(b) but highly deformed compared to Fig. 1(a) . In terms of strains, the elastic energy corresponding to small deviations of each configuration from the target one should then be calculated by comparing Fig. 1(c) to Fig. 1(a) , and Fig. 1(d) to Fig. 1(b) . The elastic energy of a configuration of a tissue, characterized by an actual metric g , is then given by E T ( g ) = min g ¯ A ∈ G A min g ¯ P ∈ G P [ E P ( g , g ¯ P ) + E A ( g , g ¯ A ) ] . (6) In other words, we obtain the energy of a given configuration g by finding the target metric that minimizes the shape energy given in Eq. (3) . In practice, we do this by minimizing the explicit expressions given in Eq. (5) with respect to ε A and ε P . This leads to algebraic equations, rather than the Euler-Lagrange differential equations arising from direct minimization with respect to g . It is important to note that since external loads change the preferred target metrics, the Lamé coefficients in Eq. (4) do not alone determine the elastic response of the tissue. 1 10.1103/PhysRevLett.120.268105.f1 FIG. 1. Strain measurement with respect to a family of reference configurations in a tissue that penalizes only perimeter discrepancies. (a) and (b) are two configurations with the same perimeter P 0 ; (c) and (d) are two deformed states. Strain should be measured by comparing each deformed configuration to the closest reference configuration; thus, (c) should be compared to (a), and (d) compared to (b). There are two classes of solutions. The first class corresponds to configurations for which both the area and the perimeter can obtain their target values. In this case the ground state energy vanishes and the tissue behaves like an anomalously soft material. The second class of solutions corresponds to the case where there are no configurations that simultaneously satisfy the target area and perimeter, which are then said to be incompatible. The energy of the ground state in this case is finite. The tissue has a finite prestress and is “stiff” in its response to external loads. The transition between these two classes of solutions corresponds to the stiffening of tissue. It is controlled by a purely geometric incompatibility and is independent of the specific form of the energy functional or the specific measure of strain. To demonstrate this, we now specialize to a lattice of hexagonal cells, as shown in the inset of Fig. 2(a) . The calculation is easily extended to other lattices of regular polygons (see the SM [40] ). To find the ground state, the symmetry of the problem allows us to neglect the off-diagonal elements of g ¯ P , g ¯ A associated with shear zero modes. The parametrization of the target metrics for the area does not depend on polygonal shape: it is given by Eq. (5) . For a hexagon whose base is oriented along the x direction, the metric g ¯ P is (see the SM [40] ) g ¯ P ( ε p ) = α p 2 [ ε p 2 0 0 3 - 2 ε p ] . (7) The unknowns characterizing the ground state are the two components g 11 and g 22 of the actual metric and the minimizers ε A and ε P of the target metrics. For hexagonal cells we find that area and perimeter are compatible for s 0 > s 0 c = 8 3 ≃ 3.722 and that the ground state energy vanishes. For s 0 < s 0 c there are no compatible solutions and the ground state energy is finite, indicating that the lattice is prestressed. Note that s 0 c is simply the isoperimetric value of a regular hexagon. The nature of the solution depends only on the target shape index s 0 [9] , but not on κ A and κ P . The ground state metric, denoted g * , as well as the prestress, depend on the model parameters. The solutions are always compatible if either κ A or κ P vanish. We emphasize that while prestress is commonly associated with incompatibility between adjacent material elements, here it reflects incompatibility at the level of a single tissue element. 2 10.1103/PhysRevLett.120.268105.f2 FIG. 2. (a) Dimensionless energy E ( Δ ) versus strain on a linear scale for s 0 = 3.874 (compatible, green), s 0 = 3.722 (at threshold, red), and s 0 = 3.577 (incompatible, blue). (Inset) Two tissue configurations of equal minimal energy for s 0 = 3.874 . The incompatible tissue has an energy gap at zero strain, corresponding to the residual stresses associated with the onset of strain stiffening. This is shown in (b), which displays the magnitude of the force F = ∂ E / ∂ Δ as a function of strain Δ on a log-log scale for ζ = 1 and s 0 = 3.577 . The dashed line has the slope 1. Solid-solid transition.— We now consider the mechanical response of tissue to an externally applied uniaxial deformation along the x direction, with the y direction left free. The strain is defined as u = 1 2 ( g - g * ) . We let g 11 = g 11 * + Δ , and determine g 22 ( Δ ) by minimizing the energy for a fixed Δ (this includes minimizing with respect to ε A and ε P ). We then evaluate the energy of this configuration to obtain E T ( Δ ) , shown in Fig. 2(a) . Here and below, energies are measured in units of κ P P 0 , and g is rescaled by α A , corresponding to lengths measured in units of A 0 . The model is then fully characterized by two dimensionless parameters: s 0 and the ratio ζ = ( κ A A 0 / κ P P 0 ) (see SM [40] ). The deformation energy of the compatible tissue, displayed in Fig. 2(a) (green curve), is a nonmonotonic function of Δ , with two degenerate minima. This can be understood by noting that, for these parameters, the minimization described in Eq. (6) yields the two degenerate target configurations shown in the inset of Fig. 2(a) (labeled V and H ). These can be transformed into each other via a uniaxial strain. The energy E ( Δ ) shown in Fig. 2(a) is calculated by measuring the deformation of the lattice relative to the V configuration. When Δ = 0 , the system is in the V ground state and has zero energy. As Δ is increased, the lattice deforms relative to the V configuration and the energy increases, eventually reaching a new zero when the deformed configuration becomes identical to the H ground state. The energy of the incompatible tissue [Fig. 2(a) , blue curve] is also quadratic at very small Δ but has an energy gap at Δ = 0 . The red curve corresponds to the critical state at s 0 = s 0 c = 3.722 . In Fig. 2(b) we show that both compatible and incompatible tissues respond linearly at small strain, albeit with different stiffnesses. Tissue stiffness may be quantified by defining an effective Young’s modulus Y that measures the response to stretching by fitting the energy just beyond the minimum to a quadratic form (see the SM [40] ). In compatible tissue, Y is very small for small strain, but it becomes appreciable once the tissue has settled in the minimum at finite Δ . The effective Young’s modulus is then calculated by a quadratic fit in the region beyond this second minimum. The effective Young’s modulus shown in Fig. 3(a) shows the onset of stiffening at s 0 = s 0 c . 3 10.1103/PhysRevLett.120.268105.f3 FIG. 3. (a) Effective Young’s modulus, as defined by the response to stretching, versus target shape index s 0 for ζ = 0.125 , 0.25, 0.5, 1, 2, 4, 8 of a uniform tissue of hexagonal cells. A transition occurs at the critical value s 0 c = 8 3 (vertical dashed line), below which the perimeter and area are incompatible. (b) Energy as a function of strain Δ x x for various values of fixed transverse strain Δ y y = 0.15 , 0.6, 0.75, 0.9, 1.05, 1.35, for an incompatible tissue with s 0 = 3.438 and ζ = 1 . The essential minimization in Eq. (6) renders our model tissue nonlinear at large Δ . This nonlinearity is further highlighted by noting that straining the system along a specific direction affects the mechanical response both along that direction and in the transverse direction. This is shown in Fig. 3(b) , where we plot the elastic energy as a function of a strain Δ x x along the x direction for various fixed strains Δ y y along the y direction. The corresponding effective moduli and their dependence on Δ y y , as well as similar figures for compatible and critically compatible tissue, are shown in the SM [40] . The connection between geometric incompatibility in cellular tissues and the emergence of stiffness is made even clearer by rewriting Eqs. (3) and (6) in terms of a single effective target metric g ¯ 0 . Completing the square gives E eff = min G A , G P 1 2 ∫ Ω A 0 α β γ δ u α β 0 u γ δ 0 | g ¯ 0 | d 2 x + E res . (8) The full expression for g ¯ 0 in terms of g ¯ P and g ¯ A is given in the SM [40] . The residual energy E res ( g ¯ 0 , g ¯ P , g ¯ A ) is independent of the actual configuration g and depends only on the families of target metrics and elastic moduli. In the absence of external loads, the actual metric that minimizes the energy is g * = g ¯ 0 , and the optimal target metrics g ¯ A * , g ¯ P * , are found by minimizing the residual energy. If the area and perimeter are compatible, E res = 0 , and the equilibrium target metrics are degenerate. For an incompatible area and perimeter, deformations from the equilibrium target metrics are no longer zero modes, as shown by expanding E res in the vicinity of the minimizers g ¯ A * and g ¯ P * . In this case variations of the reference metrics within the families G A , G P cost a finite energy, resulting in a gapped ground state. Stiffness emerges purely as a result of this geometric frustration, as suggested by previous numerical analysis of VMs [27] . Discussion.— Using the geometric formulation of elasticity, we have proposed a continuum energy for a two-dimensional tissue that incorporates the physics of well-established cellular tissue vertex models, excluding T 1 transformations, and accounts for zero modes associated with area- and perimeter-preserving deformations. We have shown that this energy yields two classes of ground states tuned by the target cell shape index s 0 . For s 0 > s 0 c the tissue is soft, with zero modes associated with a family of degenerate target metrics. For s 0 < s 0 c one obtains a stiffer nonlinear solid with residual stress at zero external deformation as a result of geometric incompatibility. An onset of stiffness accompanied by the appearance of a residual stress was also recently demonstrated numerically in a disordered Voronoi model in 3D [27] . Our model is purely elastic and considers a regular lattice. The increase in stiffness is distinct from the SL transition previously reported in the literature for disordered tilings and associated with the onset of finite energy barriers for T 1 transitions. The softening of the hexagonal lattice at s 0 = 3.722 can effectively lower the energy barriers for T 1 transformations. If these are allowed, the tissue may then melt. We expect this melting to occur at s 0 = 3.722 for a regular hexagonal lattice. The finite energy cost of local deformations of the target geometry resulting from incompatibility is directly linked to the extensive literature on the plasticity of solids (Refs. [42,43] and references therein). In that setting, changes in the target geometry are interpreted as plastic deformations [44] . Fixed isotropic inclusions in amorphous solids, for example, are known to strengthen the material by increasing the yield strain required for the formation of system-spanning shear bands [22] . The analogy between variations of the reference geometry and anelastic deformations in amorphous solids suggests that introducing isotropic sources of stresses, such as inhomogeneities in the target area, may strengthen cellular tissue. This could be tested numerically. Finally, the formalism presented here is general and can be extended to spatially inhomogeneous target metrics to describe disordered cellular structures, or even time-varying metrics to allow for local growth. An open question is whether the energy gap obtained here for ordered lattices will persist in disordered tilings. While a direct extension of our model to disordered tissue is challenging, exact and approximate analytical solutions for certain realistic problems of spatially varying shape-parameter and nonhomogenous deformations are tractable, and they will be presented in a future publication. We thank Max Bi, Matthias Merkel, and Lisa Manning for valuable discussions, and Yohai Bar-Sinai and Daniel Sussman for a critical reading of the manuscript. We acknowledge support from the National Science Foundation at Syracuse University through Grants No. DMR-1435794 (M. M., M. J. B.) and No. DMR-1609208 (M. C. M.), at Harvard through Grant No. DMR-1435999 (M. M.) and at KITP through Grant No. PHY-1125915 (M. M., M. J. B.). M. M. acknowledges the USIEF Fulbright program. M. M. and M. C. M. acknowledge support from the Simons Foundation Targeted Grant in the Mathematical Modeling of Living Systems 342354. All of the authors thank the Syracuse Soft Matter Program for its support and the KITP for hospitality during the completion of this work. [1] 1 T. E. Angelini , E. Hannezo , X. Trepat , J. J. Fredberg , and D. A. Weitz , Phys. Rev. Lett. 104 , 168104 ( 2010 ). PRLTAO 0031-9007 10.1103/PhysRevLett.104.168104 [2] 2 T. E. Angelini , E. Hannezo , X. Trepat , M. Marquez , J. J. Fredberg , and D. A. Weitz , Proc. Natl. Acad. Sci. U.S.A. 108 , 4714 ( 2011 ). PNASA6 0027-8424 10.1073/pnas.1010059108 [3] 3 K. D. Nnetu , M. Knorr , J. Käs , and M. Zink , New J. Phys. 14 , 115012 ( 2012 ). NJOPFM 1367-2630 10.1088/1367-2630/14/11/115012 [4] 4 H. Honda , J. Theor. Biol. 72 , 523 ( 1978 ). JTBIAP 0022-5193 10.1016/0022-5193(78)90315-6 [5] 5 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 [6] 6 L. Hufnagel , A. A. Teleman , H. Rouault , S. M. Cohen , and B. I. Shraiman , Proc. Natl. Acad. Sci. U.S.A. 104 , 3835 ( 2007 ). PNASA6 0027-8424 10.1073/pnas.0607134104 [7] 7 T. Nagai and H. Honda , Philos. Mag. B 81 , 699 ( 2001 ). PMABDJ 1364-2812 10.1080/13642810108205772 [8] 8 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 [9] 9 D. Bi , X. Yang , M. C. Marchetti , and M. L. Manning , Phys. Rev. X 6 , 021011 ( 2016 ). PRXHAE 2160-3308 10.1103/PhysRevX.6.021011 [10] 10 T. Su and G. Lan , arXiv:1610.04254 . [11] 11 D. Bi , J. H. Lopez , J. Schwarz , and M. L. Manning , Soft Matter 10 , 1885 ( 2014 ). SMOABF 1744-683X 10.1039/c3sm52893f [12] 12 D. Bi , J. Lopez , J. Schwarz , and M. L. Manning , Nat. Phys. 11 , 1074 ( 2015 ). NPAHAX 1745-2473 10.1038/nphys3471 [13] 13 D. Barton , S. Henkes , C. Weijer , and R. Sknepnek , PLoS Comput. Biol. 13 , e1005569 ( 2017 ). PCBLBG 1553-7358 10.1371/journal.pcbi.1005569 [14] 14 S. Banerjee , K. J. C. Utuje , and M. C. Marchetti , Phys. Rev. Lett. 114 , 228101 ( 2015 ). PRLTAO 0031-9007 10.1103/PhysRevLett.114.228101 [15] 15 M. H. Köpf and L. M. Pismen , Soft Matter 9 , 3727 ( 2013 ). SMOABF 1744-683X 10.1039/c3sm26955h [16] 16 J. Ranft , M. Basan , J. Elgeti , J.-F. Joanny , J. Prost , and F. Jülicher , Proc. Natl. Acad. Sci. U.S.A. 107 , 20863 ( 2010 ). PNASA6 0027-8424 10.1073/pnas.1011086107 [17] From here on, we use the term VMs to refer to discrete models with the tissue energy given in Eq. (1) . [18] 18 S. Eran , M. Marder , and H. L. Swinney , Am. Sci. 92 , 254 ( 2004 ). AMSCAC 0003-0996 10.1511/2004.47.932 [19] 19 Y. Klein , E. Efrati , and E. Sharon , Science 315 , 1116 ( 2007 ). SCIEAS 0036-8075 10.1126/science.1135994 [20] 20 S. Armon , E. Efrati , R. Kupferman , and E. Sharon , Science 333 , 1726 ( 2011 ). SCIEAS 0036-8075 10.1126/science.1203874 [21] 21 R. Kupferman , M. Moshe , and J. P. Solomon , Arch. Ration. Mech. Anal. 216 , 1009 ( 2015 ). AVRMAW 0003-9527 10.1007/s00205-014-0825-y [22] 22 H. G. E. Hentschel , M. Moshe , I. Procaccia , and K. Samwer , Philos. Mag. 96 , 1399 ( 2016 ). PMHABF 1478-6435 10.1080/14786435.2016.1163433 [23] 23 V. Blåsjö , Am. Math. Mon. 112 , 526 ( 2005 ). AMMYAE 0002-9890 10.1080/00029890.2005.11920227 [24] 24 K. E. Kasza , A. C. Rowat , J. Liu , T. E. Angelini , C. P. Brangwynne , G. H. Koenderink , and D. A. Weitz , Curr. Opin. Cell Biol. 19 , 101 ( 2007 ). COCBE3 0955-0674 10.1016/j.ceb.2006.12.002 [25] 25 I. Levental , P. C. Georges , and P. A. Janmey , Soft Matter 3 , 299 ( 2007 ). SMOABF 1744-683X 10.1039/B610522J [26] 26 P. Fernández , P. A. Pullarkat , and A. Ott , Biophys. J. 90 , 3796 ( 2006 ). BIOJAU 0006-3495 10.1529/biophysj.105.072215 [27] 27 M. Merkel and M. L. Manning , New J. Phys. 20 , 022002 ( 2018 ). NJOPFM 1367-2630 10.1088/1367-2630/aaaa13 [28] 28 S. Garcia , E. Hannezo , J. Elgeti , J.-F. Joanny , P. Silberzan , and N. S. Gov , Proc. Natl. Acad. Sci. U.S.A. 112 , 15314 ( 2015 ). PNASA6 0027-8424 10.1073/pnas.1510973112 [29] 29 S. M. Zehnder , M. Suaris , M. M. Bellaire , and T. E. Angelini , Biophys. J. 108 , 247 ( 2015 ). BIOJAU 0006-3495 10.1016/j.bpj.2014.11.1856 [30] 30 F. Giavazzi , C. Malinverno , S. Corallino , F. Ginelli , G. Scita , and R. Cerbino , J. Phys. D 50 , 384003 ( 2017 ). JPAPBE 0022-3727 10.1088/1361-6463/aa7f8e [31] 31 D. M. Sussman and M. Merkel , Soft Matter 14 , 3397 ( 2018 ). SMOABF 1744-6848 10.1039/C7SM02127E [32] 32 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 [33] 33 B. Audoly and Y. Pomeau , Elasticity and Geometry: From Hair Curls to the Non-linear Response of Shells ( Oxford University , New York, 2010 ). [34] 34 W. T. Koiter , Proc. K. Ned. Akad. Wet., Ser. B: Phys. Sci. 69 , 1 ( 1966 ). [35] 35 E. Efrati , E. Sharon , and R. Kupferman , J. Mech. Phys. Solids 57 , 762 ( 2009 ). JMPSA8 0022-5096 10.1016/j.jmps.2008.12.004 [36] For simple elastic solids, g ¯ is Euclidean and, in Cartesian coordinates, can be written as g ¯ 11 = g ¯ 22 = 1 and g ¯ 12 = g ¯ 21 = 0 . [37] 37 E. Kröner , in Physics of Defects , Les Houches Summer School, Session XXXV , edited by R. Balian , M. Kleman , and J. P. Poirier ( North-Holland , Amsterdam, 1981 ). [38] 38 M. Moshe , I. Levin , H. Aharoni , R. Kupferman , and E. Sharon , Proc. Natl. Acad. Sci. U.S.A. 112 , 10873 ( 2015 ). PNASA6 0027-8424 10.1073/pnas.1506531112 [39] 39 E. Sharon and E. Efrati , Soft Matter 6 , 5693 ( 2010 ). SMOABF 1744-683X 10.1039/c0sm00479k [40] 40 See Supplemental Material at http://link.aps.org/supplemental/10.1103/PhysRevLett.120.268105 for detailed derivation of the model, and for a Mathematica code recovering the results presented in the paper. [41] 41 H. S. Seung and D. R. Nelson , Phys. Rev. A 38 , 1005 ( 1988 ). PLRAAN 0556-2791 10.1103/PhysRevA.38.1005 [42] 42 R. Dasgupta , H. G. E. Hentschel , and I. Procaccia , Phys. Rev. Lett. 109 , 255502 ( 2012 ). PRLTAO 0031-9007 10.1103/PhysRevLett.109.255502 [43] 43 R. Dasgupta , H. G. E. Hentschel , and I. Procaccia , Phys. Rev. E 87 , 022810 ( 2013 ). PRESCM 1539-3755 10.1103/PhysRevE.87.022810 [44] 44 E. Efrati , E. Sharon , and R. Kupferman , Soft Matter 9 , 8187 ( 2013 ). SMOABF 1744-683X 10.1039/c3sm50660f

PY - 2018/6/29

Y1 - 2018/6/29

N2 - We study the mechanical behavior of two-dimensional cellular tissues by formulating the continuum limit of discrete vertex models based on an energy that penalizes departures from a target area A0 and a target perimeter P0 for the component cells of the tissue. As the dimensionless target shape index s0=(P0/A0) is varied, we find a transition from a soft elastic regime for a compatible target perimeter and area to a stiffer nonlinear elastic regime frustrated by geometric incompatibility. We show that the ground state in the soft regime has a family of degenerate solutions associated with zero modes for the target area and perimeter. The onset of geometric incompatibility at a critical s0c lifts this degeneracy. The resultant energy gap leads to a nonlinear elastic response distinct from that obtained in classical elasticity models. We draw an analogy between cellular tissues and anelastic deformations in solids.

AB - We study the mechanical behavior of two-dimensional cellular tissues by formulating the continuum limit of discrete vertex models based on an energy that penalizes departures from a target area A0 and a target perimeter P0 for the component cells of the tissue. As the dimensionless target shape index s0=(P0/A0) is varied, we find a transition from a soft elastic regime for a compatible target perimeter and area to a stiffer nonlinear elastic regime frustrated by geometric incompatibility. We show that the ground state in the soft regime has a family of degenerate solutions associated with zero modes for the target area and perimeter. The onset of geometric incompatibility at a critical s0c lifts this degeneracy. The resultant energy gap leads to a nonlinear elastic response distinct from that obtained in classical elasticity models. We draw an analogy between cellular tissues and anelastic deformations in solids.

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

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

U2 - 10.1103/PhysRevLett.120.268105

DO - 10.1103/PhysRevLett.120.268105

M3 - Article

C2 - 30004729

AN - SCOPUS:85049369940

VL - 120

JO - Physical Review Letters

JF - Physical Review Letters

SN - 0031-9007

IS - 26

M1 - 268105

ER -