Mean-field elastic moduli of a three-dimensional, cell-based vertex model

Abstract

The mechanics of a foam depends on bubble shape, bubble network topology, and the material at hand, be it metallic or polymeric, for example. While the shapes of bubbles are the consequence of minimizing surface area for a given bubble volume in a space-filling packing, if one were to consider biological tissue as a foam-like material, the zoology of observed shapes of cells perhaps motivates different energetic contributions. Building on earlier two-dimensional results, here, we focus on a mean field approach to obtain the elastic moduli for an ordered, three-dimensional vertex model. We use the space-filling shape of a truncated octahedron and an energy functional containing a restoring surface area spring and a restoring volume spring. The tuning of the three-dimensional shape index exhibits a rigidity transition via a compatible-incompatible transition. Specifically, for smaller shape indices, both the target surface area and volume cannot be achieved, while beyond some critical value of the three-dimensional shape index, they can be, resulting in a zero-energy state. In addition to analytically determining the location of the transition in mean field, we find that the rigidity transition and the elastic moduli depend on the parameterization of the cell shape. This parameterization effect is more pronounced in three dimensions than in two dimensions given the zoology of shapes that a polyhedron can take on (as compared to a polygon). We also uncover nontrivial dependence of the elastic moduli on the deformation protocol in which some deformations result in affine motion of the vertices, while others result in nonaffine motion. Such dependencies on the shape parameterization and deformation protocol give rise to a nontrivial shape landscape and, therefore, nontrivial mechanical response even in the absence of topology changes.