Geometry Under Compression–Decompression: A Soft-Matter Atlas

Preface This document grew out of a line of work in which the familiar Einstein field equation, Gµν = 8πTµν, was replaced by a simpler rule designed to describe how shapes in soft and biological matter emerge from mechanical constraints rather than from spacetime curvature. In this alternative viewpoint, geometry is governed not by the Einstein tensor, but by a balance between local compression and decompression forces. This was expressed through a Compression–Decompression Equation of the form ∇·C= Λ, where Crepresents the net compression–decompression field acting on a material interface, and Λ encodes the constraints imposed by tension, bending, pressure differences, flow, or phase separation. Changing the terms inside Λ changes the resulting geometry. From this starting point, a simple question guided the exploration: If different physical “settings” modify the balance encoded in ∇·C= Λ, can we classify the shapes that naturally arise from each setting? Working through examples showed that the answer is yes. A wide variety of geometries—spheres, cylinders, sheets, helices, folds, branching networks, toroidal shapes, vesicles, multi-phase domains, minimal surfaces, foams, filament networks, and many others— appear as solutions corresponding to particular choices of tension, bending rigidity, pressure, adhesion, anisotropy, or flow. This document organizes twenty-five such geometry families into a coherent atlas. For each class, the sections describe: •the physical setting that selects the geometry, •the governing equation or variational principle, •examples where the shape appears in nature or materials, •and the conditions under which the shape is stable or unstable. The goal is not to propose a universal theory, but to provide a clear, practical mapping of how different shapes follow from changes in basic physical constraints. The compression–decompression idea serves as a unifying thread: change the constraints, change the geometry. This atlas summarizes what emerged from that exploration: a concise way to understand why so many forms in soft and biological matter are variations of the same underlying principles. 1 Geometry Under Compression–Decompression: A Soft-Matter View Ricardo Miguel Machado Fernandes Contents 1 Setup: Energy, Compression, and the Spherical Bias 7 1.1 First variation and the Young–Laplace equation ............... 7 1.2 Spherical bias: constant mean curvature implies sphere ........... 8 2 Breaking Spherical Symmetry: Adding Geometric “Settings” 8 2.1 Preferred curvature as a symmetry-breaking mechanism ........... 8 2.2 Shape equation with preferred curvature ................... 8 2.3 Example: Cylindrical tubes from C0= 0 ................... 9 2.4 Interpretation .................................. 9 3 Anisotropic Tension: Direction-Dependent “Settings” 9 3.1 Surface tension as a tensor ........................... 9 3.2 Example: Formation of a planar sheet .................... 10 3.3 Example: Cylindrical shapes from anisotropy ................. 10 3.4 Cross-like tension patterns ........................... 10 3.5 Interpretation .................................. 11 4 One-Dimensional Geometry: Rods, Twist, and Helices 11 4.1 Kirchhoff rod energy .............................. 11 4.2 Straight rod vs. helical rod ........................... 11 4.3 Geometry of the equilibrium helix ....................... 12 4.4 Interpretation .................................. 12 5 Crystals: Faceted Shapes from Anisotropic Surface Energy 12 5.1 Surface energy with orientation dependence ................. 12 5.2 The Wulff shape ................................. 12 5.3 Example: Cubic symmetry ........................... 13 5.4 Interpretation .................................. 13 6 Branching and Fractals: Minimizing Transport Cost under Constraints 13 6.1 The transport cost functional ......................... 14 6.2 Murray’s law .................................. 14 6.3 Why branching reduces energy ......................... 14 6.4 Fractal scaling .................................. 14 6.5 Interpretation .................................. 15 2 7 Toroidal Geometry: Curvature Balance and Topological Constraints 15 7.1 Geometry of a standard torus ......................... 15 7.2 Isotropic surface tension cannot produce a torus ............... 16 7.3 Curvature elasticity enables toroidal shapes ................. 16 7.4 Toroidal vesicles ................................. 16 7.5 Interpretation .................................. 16 8 Wrinkling and Folding: Instabilities of Compressed Sheets 17 8.1 Elastic energy of a thin sheet ......................... 17 8.2 Onset of wrinkling ............................... 17 8.3 Wrinkle wavelength ............................... 17 8.4 From wrinkles to folds ............................. 18 8.5 Interpretation .................................. 18 9 Ellipsoids and External Fields: Breaking Symmetry from Outside 18 9.1 Droplet under gravity .............................. 18 9.2 Bond number and shape deformation ..................... 19 9.3 Electric field deformation ............................ 19 9.4 Rotating droplets ................................ 19 9.5 Interpretation .................................. 19 10 Biconcave and Multi-Lobed Shapes: Area–Volume Constrained Minimizers 20 10.1 Fixed area and fixed volume .......................... 20 10.2 Reduced volume as a shape parameter .................... 20 10.3 Shape equation for closed vesicles ....................... 20 10.4 Why the discocyte minimizes energy ..................... 21 10.5 Stomatocytes and multi-lobed shapes ..................... 21 10.6 Interpretation .................................. 21 11 Multi-Phase Geometry: Janus Droplets and Compartmentalized Shapes 21 11.1 Two-phase droplet with interfacial energy .................. 22 11.2 Contact angle condition ............................ 22 11.3 Geometric regimes ............................... 22 11.4 Shape equation for the Janus droplet ..................... 22 11.5 Examples of multi-phase equilibrium shapes ................. 23 11.6 Interpretation .................................. 23 12 Reaction–Diffusion Patterns: Turing Mechanisms Interacting with Curvature 23 12.1 Reaction–diffusion equations on curved surfaces ............... 24 12.2 Curvature coupling and geometric feedback .................. 24 12.3 Pattern selection on curved surfaces ...................... 24 12.4 Example: stripe and spot formation on vesicles ............... 25 12.5 Interpretation .................................. 25 3 13 Minimal Surfaces: Catenoids, Helicoids, and Soap Films 26 13.1 The minimal-surface equation ......................... 26 13.2 Catenoid: minimal surface between two rings ................. 26 13.3 Helicoid: minimal surface with screw symmetry ............... 26 13.4 Minimal surfaces in biological and physical systems ............. 27 13.5 Minimal surfaces as limiting cases of the compressive framework . . . . . . 27 14 Foams and Space-Filling Geometries: Plateau Rules and Kelvin Cells 27 14.1 Plateau’s laws for foam geometry ....................... 28 14.2 Volume constraints and bubble shapes .................... 28 14.3 Kelvin cells and optimal space-filling structures ............... 28 14.4 Foams in biological contexts .......................... 28 14.5 Interpretation .................................. 29 15 Filament Networks and Tension-Web Geometries: Actin, ECM, and Cytoskeletal Meshes 29 15.1 Force balance on a filament node ....................... 29 15.2 Continuum limit: stress tensor of a filament network ............ 30 15.3 Geometry from network anisotropy ...................... 30 15.4 Energetics of filament bending and stretching ................ 30 15.5 Biological and physical examples ....................... 30 15.6 Interpretation .................................. 31 16 Polyhedral and Geodesic Geometries: Viral Capsids, Fullerenes, and Tensegrity Shells 31 16.1 Icosahedral symmetry as a fundamental organizing principle ........ 31 16.2 Discrete elastic energy for triangulated surfaces ............... 32 16.3 Viral capsid geometry ............................. 32 16.4 Fullerenes and carbon polyhedra ........................ 32 16.5 Tensegrity shells ................................. 32 16.6 Interpretation .................................. 33 17 Morphogenetic Flow Geometries: Shear, Vortices, and Rotational Patterning 33 17.1 Flow-induced deformation of membranes ................... 33 17.2 Rotational flows and vortex geometries .................... 34 17.3 Shear-alignment and elongation ........................ 34 17.4 Flow-induced instabilities and pattern formation ............... 34 17.5 Active morphogenetic flows .......................... 35 17.6 Interpretation .................................. 35 18 Shear-Induced Geometries: Elongation, Banding, and Teardrop Morphologies 36 18.1 Shear deformation of a droplet ......................... 36 18.2 Small-deformation analysis: Taylor deformation parameter ......... 36 18.3 Teardrop geometry under high shear ..................... 37 18.4 Shear banding and layered geometries ..................... 37 18.5 Vesicle inclination and tank-treading ..................... 37 18.6 Interpretation .................................. 38 4 19 Wetting and Adhesion Geometries: Sessile Droplets, Contact Angles, and Capillary Bridges 38 19.1 Sessile droplets and the contact angle condition ............... 38 19.2 Geometry of a spherical cap .......................... 39 19.3 Wetting ridges on soft substrates ....................... 39 19.4 Capillary bridge geometries .......................... 39 19.5 Pinned contact lines and asymmetric shapes ................. 39 19.6 Interpretation .................................. 40 20 Shell Buckling and Core–Shell Deformation: Indented Spheres, Dimples, and Pressurized Shells 40 20.1 Elastic energy of a thin shell .......................... 40 20.2 Buckling of a pressurized spherical shell .................... 41 20.3 Geometry of a single dimple .......................... 41 20.4 Multi-dimple and ridge formation ....................... 41 20.5 Core–shell systems and asymmetric shapes .................. 42 20.6 Indentation by localized forces ......................... 42 20.7 Interpretation .................................. 42 21 Tubular Instabilities: Rayleigh–Plateau Breakup, Pearling, and Unduloid Geometries 43 21.1 Rayleigh–Plateau instability .......................... 43 21.2 Pearling of lipid tubes ............................. 43 21.3 Unduloids as equilibrium shapes ........................ 44 21.4 Axial tension and geometric bifurcations ................... 44 21.5 Viscoelastic and active tubular instabilities .................. 44 21.6 Interpretation .................................. 45 22 Triply Periodic Minimal Surfaces: Gyroid, Schwarz P, and Schwarz D Geometries 45 22.1 Minimal-surface condition ........................... 45 22.2 Schwarz P surface (primitive surface) ..................... 46 22.3 Schwarz D surface (diamond surface) ..................... 46 22.4 Gyroid surface: chiral TPMS .......................... 46 22.5 Energetic origin of TPMS geometries ..................... 47 22.6 Comparison with foams and membranes ................... 47 22.7 Interpretation .................................. 47 23 Porous and Sponge Geometries: Percolation Networks, Trabecular Bone, and Reticulated Structures 48 23.1 Percolation networks and critical geometry .................. 48 23.2 Trabecular bone as an optimal porous structure ............... 48 23.3 Reticulated polymers and stochastic foams .................. 49 23.4 Mechanical optimality and geometric scaling ................. 49 23.5 Sponge-like biological tissues .......................... 49 23.6 Interpretation .................................. 50 5 24 Polygonal Cell Tilings: Voronoi Geometry, Vertex Models, and Epithelial Packing 50 24.1 Voronoi geometry as a first approximation .................. 51 24.2 Vertex models for epithelial tissues ...................... 51 24.3 Topological constraints: Euler relation .................... 51 24.4 Mechanical equilibrium at vertices ....................... 52 24.5 Disordered and anisotropic tessellations .................... 52 24.6 Geometric transitions and rigidity ....................... 52 24.7 Interpretation .................................. 52 25 Jamming and Marginal Rigidity: Force Chains, Granular Packing, and Critical Geometries 53 25.1 Contact network and isostaticity ........................ 53 25.2 Force chains and emergent anisotropy ..................... 54 25.3 Packing geometry and Voronoi cells ...................... 54 25.4 Critical scaling at the jamming transition ................... 54 25.5 Granular architecture and anisotropic loading ................ 55 25.6 Biological jamming ............................... 55 25.7 Interpretation .................................. 55 6 1 Setup: Energy, Compression, and the Spherical Bias We consider a soft, deformable body (e.g. a liquid droplet or fluid membrane) with surface ∂Ω enclosing a volume V= Vol(Ω). The simplest macroscopic energy functional is E[Ω] = σ A[Ω] + p V [Ω],(1) where •σ > 0 is the (isotropic) surface tension, •A[Ω] is the surface area of ∂Ω, •pis a Lagrange multiplier enforcing fixed volume (physically: a pressure difference), •V[Ω] is the enclosed volume. We are interested in shapes Ω that minimize E[Ω] under a fixed volume constraint. 1.1 First variation and the Young–Laplace equation Let the surface ∂Ω be deformed by a small normal displacement δf(x)n(x), where nis the outward unit normal and δf is a scalar field on the surface. The first variations of area and volume are (standard results from differential geometry) δA =−2Z∂Ω H δf dA, (2) δV =Z∂Ω δf dA, (3) where His the mean curvature of the surface. The first variation of the energy (1) is δE =σ δA +p δV =Z∂Ω−2σH +pδf dA. (4) For Ω to be an equilibrium shape, we require δE = 0 for all variations δf. This implies that the integrand must vanish: −2σH +p= 0 =⇒2σH =p. (5) Equation (5) is the Young–Laplace relation: the mean curvature Hof the surface is everywhere constant and proportional to the pressure jump p. 7 1.2 Spherical bias: constant mean curvature implies sphere A classical result in differential geometry (Alexandrov’s theorem) states that any compact, embedded surface in R3with constant mean curvature must be a sphere. Combined with the Young–Laplace relation (5), this implies: Proposition (Spherical bias). For a body with isotropic surface tension σand fixed volume, any equilibrium shape that minimizes the energy (1)must be a sphere. In physical terms: under isotropic compression–decompression mediated by surface tension, the preferred (energy-minimizing) shape is spherical. Non-spherical shapes cannot be stable minima without additional structure or anisotropy. 2 Breaking Spherical Symmetry: Adding Geometric “Settings” The spherical bias described in Section 1 arises from the fact that the energy functional (1) contains only isotropic surface tension. To obtain non-spherical equilibrium shapes, the energy must contain additional geometric terms—“settings”—which introduce anisotropy or preferred curvature. In this section we introduce the simplest such modification and derive its consequences. This shows explicitly how non-spherical geometries can arise from soft-matter energy minimization. 2.1 Preferred curvature as a symmetry-breaking mechanism Many soft biological and physical interfaces possess an internal microstructure (lipids, proteins, polymers, anisotropic crystals) that imposes a preferred curvature on the surface. A standard model capturing this effect is the Helfrich energy: EH[Ω] = Z∂Ωhσ+κ 22H−C02idA +p V [Ω],(6) where •κis the bending rigidity, •C0is the spontaneous curvature (a geometric “setting”), •His the mean curvature, •the term (2H−C0)2penalizes deviation from preferred curvature. 2.2 Shape equation with preferred curvature Varying the energy (6) with respect to normal displacements again yields a curvature-force balance. The resulting Euler–Lagrange equation is 2σH −κ(2H−C0)2H2−2K+C0H=p, (7) where Kis the Gaussian curvature. 8 2.3 Example: Cylindrical tubes from C0= 0 For a perfect cylinder of radius R, the two principal curvatures are κ1=1 R, κ2= 0, so the mean and Gaussian curvatures are H=1 2R, K = 0. Substituting these into the Helfrich energy (6) and minimizing with respect to Ryields Rtube =1 C0 .(8) Thus a tube of radius 1/C0is energetically preferred whenever C0= 0. 2.4 Interpretation The key consequence is that introducing a preferred curvature C0changes the equilibrium geometry from a sphere (constant H) to a family of possible shapes, including cylinders and planar lamellae. Such non-spherical geometries require the presence of an internal setting (here C0) that breaks isotropy. In physical terms: the system must possess a soft, reconfigurable medium in which curvature can be redistributed. Only then can the membrane explore shape-space and settle into a non-spherical configuration. 3 Anisotropic Tension: Direction-Dependent “Settings” In Sections 1–2 the surface tension was taken to be isotropic: σ= constant. This produces spherical equilibrium shapes. To obtain non-spherical shapes even in the absence of preferred curvature, the tension must depend on direction. This is the geometric meaning of anisotropy: the interface resists deformation differently along different surface directions. 3.1 Surface tension as a tensor Instead of a scalar σ, we now treat surface tension as a rank-2 tensor σacting on tangent directions of the surface. In principal curvature coordinates this reduces to two principal tensions: T1, T2. The generalized Laplace law becomes ∆p=T1 R1 +T2 R2 ,(9) where R−1 1=κ1and R−1 2=κ2are the principal curvatures. 9 7.2 Isotropic surface tension cannot produce a torus If surface tension were isotropic (Section 1), the equilibrium shape would require constant mean curvature: 2σH =p. Because H(θ) for a torus is not constant (cf. (16)), a torus cannot satisfy this equation for any σand p. Thus: No torus is an equilibrium shape under isotropic tension. 7.3 Curvature elasticity enables toroidal shapes If the membrane has bending elasticity (Helfrich-type energy), the total energy becomes E=Z∂Ωhσ+κ 2(2H−C0)2idA +pV. (18) The additional curvature term allows non-spherical shapes to be equilibria. For suitable values of κand C0, the torus becomes a local or global energy minimum. 7.4 Toroidal vesicles In lipid vesicles, the interplay between bending energy and fixed area-to-volume ratio creates equilibrium shapes of non-zero genus. When the reduced volume v=3V 4πR3 0 is sufficiently small and the bending modulus κis large, the Helfrich energy (18) admits toroidal minima. This phenomenon has been observed experimentally and predicted theoretically: vesicles with excess area can “wrap” into toroidal geometries to reduce bending energy. 7.5 Interpretation The torus illustrates a general principle: •Isotropic compression–decompression enforces spherical geometry. •Toroidal geometry requires additional “settings”: curvature elasticity, area–volume constraints, or topological restrictions. •The transition from spherical to toroidal shapes requires a soft or fluid-like phase in which curvature can be redistributed. •Once the toroidal configuration is attained, it may be stabilized into a more rigid structure. Thus the torus fits naturally into the geometry-formation framework: it is a nonspherical shape that becomes accessible only when internal or external settings break the spherical bias. 16 8 Wrinkling and Folding: Instabilities of Compressed Sheets When a thin elastic sheet is subjected to in-plane compression, it cannot decrease its area uniformly without undergoing large strains. Instead, the sheet escapes the compressive load by deflecting out of the plane. This leads to the formation of wrinkles or folds, creating a new geometry that inherits the compression–decompression constraints. 8.1 Elastic energy of a thin sheet The standard model for an elastic sheet with thickness hcombines bending and stretching contributions. For small out-of-plane displacement w(x, y), the energy is approximately E[w] = ZD 2(∇2w)2+Y h 2ε2 ∥dx dy, (19) where •D=Y h3 12(1 −ν2)is the bending modulus, •Yis Young’s modulus, •νis Poisson’s ratio, •ε∥is the in-plane compressive strain. 8.2 Onset of wrinkling A flat configuration becomes unstable when the compressive stress σ∥satisfies σ∥> σcrit =√D Y h−1.(20) Beyond this threshold, the sheet lowers its energy by forming a periodic out-of-plane pattern: w(x, y)=Acos(kx), where kis the wrinkle wavenumber selected by minimizing (19). 8.3 Wrinkle wavelength Minimizing the energy with respect to the wavenumber kyields the characteristic wrinkle wavelength λ= 2π2D σ∥1/4 .(21) This predicts that wrinkles are finer when bending rigidity Dis small and compressive stress σ∥is large. 17 8.4 From wrinkles to folds As compression increases beyond the wrinkling regime, the sheet enters a nonlinear state. Wrinkles localize into sharp folds, producing piecewise-flat geometry separated by ridges or creases. Mathematically, this corresponds to energy concentration in narrow regions where curvature becomes large: κ∼1 ℓ, with ℓthe fold width. 8.5 Interpretation Wrinkling and folding arise because a two-dimensional sheet cannot resolve compressive loads by uniform shrinkage. Instead, the sheet explores the third dimension and adopts a geometry of oscillatory or localized curvature to lower its energy. This phenomenon fits naturally into the compressive framework: •isotropic compression favors spherical shapes; •planar geometries resist uniform contraction; •the sheet escapes by forming wrinkles or folds; •the resulting pattern is governed by material “settings” (bending modulus D, inplane stress σ∥). Such patterns can only develop when the sheet is in a soft or fluid-like state that allows reconfiguration. Once formed, they may be stabilized by subsequent stiffening. 9 Ellipsoids and External Fields: Breaking Symmetry from Outside Spherical geometry minimizes surface energy under isotropic conditions. However, external fields such as gravity, electric fields, adhesion, or shear flow introduce directional forces that break this isotropy. The resulting equilibrium shapes are no longer spherical but become ellipsoids, teardrops, or other axisymmetric configurations. 9.1 Droplet under gravity For a droplet in a gravitational field g, the pressure difference across the interface becomes height dependent: p(z)=p0+ρgz, where ρis the density and zis the vertical coordinate. The Laplace condition becomes 2σH = ∆p(z)=ρgz. (22) Because the right-hand side depends on position, the curvature Hcannot be constant. The surface therefore elongates or flattens depending on droplet size and Bond number. 18 9.2 Bond number and shape deformation The competition between gravity and surface tension is governed by the Bond number Bo = ρgR2 σ,(23) where Ris the characteristic radius. •For Bo ≪1, surface tension dominates and the shape remains nearly spherical. •For Bo ≳1, gravity dominates and significant deformation occurs, producing oblate or prolate ellipsoids. 9.3 Electric field deformation If the droplet is placed in an external electric field E, the Maxwell stress introduces an anisotropic normal pressure: pEM =ϵ 2(E2 ⊥−E2 ∥), where E⊥and E∥are the normal and tangential components of E. The shape equation becomes 2σH =pcapillary +pEM. Strong electric fields stretch the droplet along the field direction, ultimately forming a cone with semi-angle ≈49.3◦(Taylor’s cone angle). 9.4 Rotating droplets A droplet rotating with angular velocity ωexperiences a centrifugal potential that lowers pressure at the equator. The modified pressure is p(r, θ) = p0+1 2ρω2r2sin2θ. The equilibrium shapes are •oblate spheroids for moderate rotation, •dumbbell-like or multi-lobed shapes for high rotation. Rotation thus breaks spherical symmetry and produces axisymmetric non-spherical geometries. 9.5 Interpretation Ellipsoids and related shapes arise when external fields impose directional stresses that break isotropy. The droplet can no longer maintain constant mean curvature, and the resulting equilibrium geometry reflects the balance between surface tension and external forcing. Key points: 19 •gravity, electric fields, shear, and rotation break spherical symmetry; •deformation magnitude depends on nondimensional parameters such as the Bond number; •the droplet must be in a fluid-like state to reconfigure; •once symmetry is broken, new geometries (ellipsoids, cones, dumbbells) emerge as equilibrium solutions. This is another validation of the compression–decompression framework: the sphere is the default, but additional “settings” supplied by external fields lead to alternative stable geometries. 10 Biconcave and Multi-Lobed Shapes: Area–Volume Constrained Minimizers Red blood cells and many lipid vesicles adopt a characteristic biconcave “discocyte” geometry. This shape is neither spherical nor toroidal, yet it arises naturally from the minimization of bending energy under fixed area and fixed volume constraints. The same energy that produced tubes (Section 2) also produces these higher-order shapes when global constraints dominate. 10.1 Fixed area and fixed volume Let the membrane enclose volume Vand have fixed surface area A. These constraints are represented by Lagrange multipliers λAand λV. The total energy is E[Ω] = Z∂Ωhκ 2(2H−C0)2idA +λA(A−A0)+λV(V−V0).(24) 10.2 Reduced volume as a shape parameter The family of equilibrium shapes is parameterized by the reduced volume v=3V 4πR3 0 , R0=pA/4π. (25) •v= 1 corresponds to a perfect sphere. •v≈0.6 corresponds to a biconcave discocyte. •v < 0.5 produces stomatocytes (cup-shaped vesicles). 10.3 Shape equation for closed vesicles Variation of (24) yields the closed-vesicle shape equation: κ2(2H−C0)(2H2−2K+C0H)−∇2(2H−C0)=λAH+λV.(26) 20 10.4 Why the discocyte minimizes energy For reduced volume v≈0.6, a sphere cannot satisfy both constraints (fixed area A0 and fixed volume V0) without incurring large bending energy. The biconcave shape redistributes curvature such that regions of positive and negative curvature balance more efficiently. Energetically: Ediscocyte < Esphere for v≈0.6. 10.5 Stomatocytes and multi-lobed shapes For lower reduced volume (v < 0.5) or nonzero spontaneous curvature C0, the shape equation (26) predicts multi-lobed solutions: •stomatocytes (cup shapes), •pear-shaped vesicles, •multi-lobed ellipsoids, •axisymmetric “pearling” instabilities. These shapes arise purely from curvature elasticity and global constraints, with no need for external fields. 10.6 Interpretation The biconcave shape of the red blood cell is a canonical example of geometry emerging from the competition between bending energy, area constraint, and volume constraint. This mechanism embodies the compression–decompression framework: •A sphere is optimal only for v= 1. •When v < 1, spherical symmetry is incompatible with the constraints. •The membrane seeks a new geometry that satisfies constraints while minimizing curvature energy. •The resulting shapes—discocytes, stomatocytes, pearlings—are accessible only in a soft or fluid-like state where curvature can redistribute. Once established, the shape may be stabilized by subsequent structural reinforcement. 11 Multi-Phase Geometry: Janus Droplets and Compartmentalized Shapes Many biological and physical systems contain multiple coexisting phases with distinct surface energies. Examples include phase-separated organelles, double-emulsion droplets, patchy membranes, and amphiphilic aggregates. These systems give rise to new equilibrium geometries that do not occur in single-phase materials. The key idea is that each interface carries its own tension, producing a competition that defines the final geometry. 21 11.1 Two-phase droplet with interfacial energy Consider two immiscible phases Aand Bseparated by an interface Γ, all enclosed by a surface ∂Ω. The total energy is E=σAAA+σBAB+σABAAB +p V, (27) where •σA: surface tension of region Awith the environment, •σB: surface tension of region Bwith the environment, •σAB: interfacial tension between Aand B, •AA, AB, AAB: respective areas, •V: enclosed volume. 11.2 Contact angle condition At the line where the three interfaces meet, force balance gives Young’s law: σB−σA=σAB cos θ, (28) where θis the angle between phases Aand Bat the three-phase boundary. 11.3 Geometric regimes Depending on the relative magnitudes of σA,σBand σAB, three regimes occur: •σAB ≪ |σA−σB|: complete engulfing (core–shell geometry), •σAB ≫ |σA−σB|: separated domains (two separate droplets), •intermediate case: a “Janus” configuration, where the droplet divides into two lobes sharing a curved interface. 11.4 Shape equation for the Janus droplet For an axisymmetric interface under radial coordinate r(s) and height z(s), the Young– Laplace relation becomes 2σAHA=pA,2σBHB=pB,2σABHAB =pA−pB,(29) where HA,HB, and HAB are the mean curvatures of the respective surfaces. Equations (29) determine the geometry of the two lobes and the internal interface. 22 11.5 Examples of multi-phase equilibrium shapes Depending on the tensions and volume ratios, the Janus system can produce: •spherical core–shell structures, •two-lobed droplets, •dumbbell-like droplets, •flattened or bowl-shaped compartments, •multi-domain vesicles, •patchy membranes and surface patterns. 11.6 Interpretation Multi-phase systems introduce new geometric settings: each interface carries its own tension, and their competition determines curvature and topology. The geometry emerges from balancing the three surface energies and satisfying the contact-angle condition. Key points: •a sphere is optimal only when all tensions are equal, •unequal tensions produce multi-lobed or compartmentalized shapes, •the interface must be fluid or soft to reconfigure, •once formed, the geometry may be “frozen” by solidification. Thus multi-phase droplets provide another class of non-spherical shapes that arise naturally from the compression–decompression framework when additional settings are present. 12 Reaction–Diffusion Patterns: Turing Mechanisms Interacting with Curvature Beyond purely mechanical or energetic constraints, geometry can emerge from chemical pattern formation. A fundamental model for such patterning is the reaction–diffusion (RD) system introduced by Turing, in which two or more chemical species diffuse and react to produce spatially heterogeneous concentration patterns. When these patterns form on curved surfaces, such as membranes or vesicles, the resulting chemical fields influence the mechanical energy of the surface and thereby shape the geometry. Conversely, the curvature of the surface influences the RD dynamics. This feedback creates a coupled chemo-mechanical pattern-forming system that can generate spots, stripes, labyrinths, and curvature-dependent domains. 23 12.1 Reaction–diffusion equations on curved surfaces Let u(x, t) and v(x, t) be two reacting and diffusing species defined on a surface ∂Ω with Laplace– Beltrami operator ∆s. A general RD system is ∂u ∂t =Du∆su+f(u, v),(30) ∂v ∂t =Dv∆sv+g(u, v),(31) where Duand Dvare diffusion coefficients, and f, g describe local reaction kinetics. Patterns form when Du=Dv, so that differential diffusion destabilizes the homogeneous steady state. On a curved surface, the Laplace–Beltrami operator couples pattern formation directly to surface geometry: ∆su=∇s·∇su, where ∇sis the tangential gradient. 12.2 Curvature coupling and geometric feedback Reaction–diffusion patterns modify curvature whenever the chemical fields affect local material properties such as: •spontaneous curvature C0(u, v), •surface tension σ(u, v), •bending rigidity κ(u, v), •line tension between chemically distinct domains. Incorporating chemical dependence into the Helfrich energy yields E=Z∂Ωκ(u, v) 22H−C0(u, v)2+σ(u, v)dA +pV. (32) Thus chemical patterns act as “settings” that modulate curvature and tension. Regions with high umay prefer higher curvature, while regions with high vmay prefer flatter or differently curved geometry. The membrane deforms until mechanical forces and chemical gradients reach equilibrium. 12.3 Pattern selection on curved surfaces Curvature affects RD patterns in several ways: •The eigenmodes of ∆sdepend on geometry, so the pattern wavelength and symmetry are curvature-dependent. •High-curvature regions (large |H|) amplify or suppress patterns depending on the reaction kinetics. 24 •Coupling C0(u, v) to RD fields produces localized buds, dimples, and curvaturedriven segregation. For a spherical surface, pattern modes correspond to spherical harmonics Yℓm. Lower ℓmodes dominate as curvature increases. As the membrane deforms, the spectrum of the Laplace– Beltrami operator changes, shifting the RD pattern in real time. 12.4 Example: stripe and spot formation on vesicles If fand gare chosen such that the Turing instability yields striped patterns, the membrane may develop alternating bands of high and low curvature if C0or κdepend on u. In contrast, spot-like patterns produce isolated buds on the membrane. Mechanically: RD pattern =⇒C0(u, v) =⇒curvature modulation =⇒geometric deformation. This mechanism can create: •spotted vesicles, •stripe-wrapped vesicles, •labyrinthine domains, •chemically segregated compartments, •curvature-induced budding. 12.5 Interpretation Reaction–diffusion systems introduce chemical “settings” that modify the mechanical response of the membrane. In the presence of curvature elasticity, these patterns generate non-spherical and spatially modulated geometries. The curvature affects diffusion, and diffusion affects curvature, producing a tightly coupled chemo-mechanical feedback loop. Key principles: •Without curvature coupling, RD patterns form but do not deform geometry. •Without RD patterns, curvature remains homogeneous and tends toward spherical or tubular geometries. •With coupling, the system explores new geometric classes (spots, stripes, labyrinths) that serve as equilibrium or metastable states under the combined energy (32). Thus reaction–diffusion processes provide another mechanism for breaking spherical symmetry within the compression–decompression framework, generating rich patterned geometries in soft and fluid-like interfaces. 25 16.2 Discrete elastic energy for triangulated surfaces A polyhedral or triangulated surface can be modeled using the discrete Helfrich–like bending energy E=X edges e κe(1 −cos θe),(39) where θeis the dihedral angle across edge e, and κeis the bending modulus associated with that edge. Flat faces have θe= 0 and contribute no energy; curvature locks into edges and vertices. 16.3 Viral capsid geometry Viral capsids are typically composed of Tidentical triangular subunits arranged in a quasiequivalent triangulation of the sphere. Caspar and Klug showed that such triangulations are classified by T=h2+hk +k2, with integers h, k ≥0. The capsid geometry consists of: •12 pentamers at the vertices carrying positive Gaussian curvature, •(10T−12) hexamers forming nearly flat facets, •edges supporting concentrated bending energy as given by (39). This arrangement minimizes total elastic energy for a shell composed of identical subunits. 16.4 Fullerenes and carbon polyhedra Buckminsterfullerene (C60) forms a truncated icosahedron composed of 12 pentagons and 20 hexagons. The geometry is determined by: •sp2hybridization enforcing planar hexagonal tiling, •the need to close the surface requiring 12 pentagons, •curvature localized at pentagonal sites. The result is a nearly spherical polyhedron with high structural stability, deriving from geometric and energetic constraints. 16.5 Tensegrity shells Tensegrity structures combine tension-bearing cables and compression struts to form lightweight yet rigid polyhedral frameworks. Their geometry is determined by force balance: X j Fij = 0 for each node i. Biological examples include: 32 •the cytoskeleton modeled as a tensegrity network, •extracellular matrices forming geodesic-like scaffolds, •protein complexes that adopt polyhedral assemblies. 16.6 Interpretation Polyhedral and geodesic geometries arise when: •discrete building blocks enforce faceted surfaces, •curvature is concentrated at a finite set of vertices, •bending energy is minimized by flattening faces, •symmetry reduces the total energy of assembling identical units. Within the compression–decompression framework: •the sphere is the default smooth solution, •discrete or modular constraints force curvature localization, •energy minimization yields faceted or geodesic structures, •formation requires a soft or mobile phase before solidification. Thus polyhedral capsids, fullerenes, and tensegrity shells represent a geometry class where discrete symmetry, minimal bending energy, and modular assembly combine to break spherical symmetry in a predictable, mathematically structured way. 17 Morphogenetic Flow Geometries: Shear, Vortices, and Rotational Patterning Beyond static energy minimization, many biological and physical systems generate geometry through flow. When tissues, membranes, or fluid-like layers move under shear, rotation, or internal circulation, the resulting velocity fields deform the material into characteristic shapes. These morphogenetic flows represent a dynamic geometry class where shapes arise from the interplay between stresses, viscosity, and curvature in a soft medium. 17.1 Flow-induced deformation of membranes Let v(x, t) be the velocity field tangential to a membrane. The local strain rate is given by the symmetrized gradient εij =1 2∇ivj+∇jvi.(40) Shear flow produces anisotropic stretching, generating preferred directions of curvature. If bending rigidity is small compared to viscous stresses, the membrane deforms primarily according to the flow field. The membrane shape evolves according to dX dt =v∥+vnn, where vnis the normal velocity from curvature forces. 33 17.2 Rotational flows and vortex geometries Rotational or swirling flows are characterized by nonzero vorticity ω=∇×v. Regions of high vorticity induce spiral or vortex-like deformations: •swirling patterns on cell sheets, •vortex-driven rearrangements in epithelial monolayers, •ring-like deformations in fluid membranes, •spiral morphologies in active matter. In thin films, vorticity couples to curvature through the “geometric torque” τgeo ∼κ∇H, causing membranes to buckle or twist in rotating patterns. 17.3 Shear-alignment and elongation Under uniform shear, droplets, cells, and vesicles elongate along the principal strain direction. The shape tensor Qij evolves as ˙ Qij =εij −1 τQij,(41) where τis a relaxation time. Steady-state solutions of (41) produce: •ellipsoidal vesicles aligned with flow, •teardrop-shaped cells, •elongated droplets in emulsions, •anisotropic tissue domains. Shear thus acts as an external “setting” that biases geometry toward elongated forms. 17.4 Flow-induced instabilities and pattern formation Flow can destabilize uniform geometry, producing: •Kelvin–Helmholtz instabilities (wavelike interfaces), •Tollmien–Schlichting waves (shear bands), •swirling cell collectives, •ring-like vortical deformations. 34 These patterns occur when shear overcomes surface tension or elastic restoring forces. For a membrane with tension σunder shear rate ˙γ, the instability threshold is approximately ˙γ > rσ ρR3,(42) where ρis density and Ra characteristic radius. Above this threshold, the membrane cannot maintain isotropic shape and develops flow-driven geometry. 17.5 Active morphogenetic flows Biological tissues often exhibit active flows generated by motor proteins or actomyosin contraction. The active stress tensor has the form σactive ij =ζ qij, where qij is a nematic order parameter and ζa contractile coefficient. Such flows generate: •topological defects (spirals, asters), •convergent–extension flows, •rotational vortices in epithelial tissues, •swirl patterns in stem-cell colonies. These geometries arise from force generation rather than external constraints, illustrating that shape can emerge from internal flow activity within a soft medium. 17.6 Interpretation Morphogenetic flow geometries extend the compression–decompression framework into the dynamic regime. Here, shapes are not determined solely by static energy minimization but by the interplay of: •viscous stresses, •shear and vorticity, •curvature elasticity, •active internal force generation. Key principles: •Shear and rotation break spherical symmetry dynamically. •Flow alignment leads to elongated, teardrop, or ellipsoidal shapes. •Vorticity induces spiral and vortex-driven geometries. •Active flows in biological tissues generate patterns not accessible through passive minimization alone. 35 •All such geometries require a fluid or soft state to form. Thus morphogenetic flows represent a geometry class governed by dynamic stress fields, extending the range of shapes accessible under the compression–decompression paradigm. 18 Shear-Induced Geometries: Elongation, Banding, and Teardrop Morphologies Shear forces represent a directional form of compression–decompression. When a soft or fluid interface is subjected to shear, the deformation is anisotropic: material lines aligned with the flow stretch, while those perpendicular to it compress. This imbalance produces a distinct class of geometries characterized by elongation, asymmetry, and directional curvature. Shear-induced shapes arise in droplets, vesicles, cells, emulsions, tissues, and soft materials across scales. 18.1 Shear deformation of a droplet A simple model for droplet deformation under steady shear flow v= (˙γy, 0,0) gives the velocity gradient tensor ∇v=  0 ˙γ0 0 0 0 0 0 0 . Surface tension σresists deformation, while shear rate ˙γpromotes elongation. The dimensionless capillary number Ca = µ˙γR σ(43) governs the shape response, where µis viscosity and Ris the droplet radius. •Ca ≪1: droplet remains nearly spherical. •Ca ≳1: significant elongation occurs. •Ca ≫1: droplet forms a teardrop or becomes unstable. 18.2 Small-deformation analysis: Taylor deformation parameter For mild shear, the droplet deforms into an ellipsoid with deformation parameter D=L−B L+B,(44) where Land Bare major and minor axes. Taylor’s law gives D=19λ+ 16 16λ+ 16 Ca,(45) with viscosity ratio λ=µin/µout. As deformation increases, the droplet becomes progressively elongated and asymmetric. 36 18.3 Teardrop geometry under high shear At high capillary number, curvature concentrates near the rear end of the droplet. The front becomes blunt, while the tail stretches into a sharp tip. The shape approximates a “teardrop”: κtail ≫κfront. This asymmetry arises because shear continuously pulls fluid from the rear, while tension cannot restore symmetry fast enough. Teardrop shapes appear in: •vesicles in microfluidic shear, •red blood cells under fast flow, •immiscible droplets in emulsions, •soft biological tissues experiencing directional stress. 18.4 Shear banding and layered geometries Shear can induce instabilities leading to banded morphologies. In a non-Newtonian or viscoelastic medium, stress is non-monotonic in the shear rate: σ(˙γ) decreases for some ranges of ˙γ. This produces: •shear bands with different strain rates, •striped or layered geometries, •alternating thick and thin regions of a film, •localized deformation channels. Mathematically, banding corresponds to coexistence of two solutions ˙γ1and ˙γ2to the constitutive relation, linked by a Maxwell-like equal-area construction. 18.5 Vesicle inclination and tank-treading Vesicles and red blood cells in shear flow exhibit periodic rotation (tank-treading) or swinging. The inclination angle θsatisfies tan 2θ=2˙γ ω0−˙γ,(46) where ω0is a characteristic membrane relaxation frequency. These dynamics produce continuously varying shapes: •tilted ellipsoids, •oscillatory deformations, •rotating teardrop geometries. 37 18.6 Interpretation Shear-induced geometries represent a dynamic pathway for breaking spherical or symmetric shapes: •Shear stretches the material in one direction and compresses it in another. •The resulting anisotropy drives elongation or teardrop shapes. •At high shear, curvature localization produces asymmetric forms. •Shear banding generates layered or striped geometry. •Vesicle inclination and tank-treading create continuously changing anisotropic shapes. All such deformations require the material to be soft or fluid, enabling it to explore shape space before solidification or stabilization. Shear-added-deformation thus fits consistently within the compression–decompression framework by introducing directional stress as a “setting” that selects new geometric families. 19 Wetting and Adhesion Geometries: Sessile Droplets, Contact Angles, and Capillary Bridges When a soft or fluid interface interacts with a solid substrate, surface adhesion and wetting forces introduce new geometric constraints. These constraints break spherical symmetry and generate characteristic shapes such as sessile droplets, capillary bridges, wetting ridges, and pinned interfaces. The geometry is governed by the balance of surface tensions at the contact line, as expressed by Young’s and related equations. 19.1 Sessile droplets and the contact angle condition A fluid droplet resting on a solid surface forms a sessile droplet shape. Let σSL,σSG, and σLG be the surface tensions of the solid–liquid, solid–gas, and liquid–gas interfaces. At equilibrium, the contact angle θsatisfies Young’s equation: σSG −σSL =σLG cos θ. (47) This relation determines the geometry at the contact line, which in turn shapes the entire droplet: •θ < 90◦: droplet spreads (wetting). •θ > 90◦: droplet beads (non-wetting). •θ= 90◦: neutral wetting. The sessile droplet shape satisfies the modified Young–Laplace equation 2σLGH= ∆p(z), with boundary condition (47) enforced at the solid surface. 38 19.2 Geometry of a spherical cap For weak gravity (Bond number ≪1), the sessile droplet approximates a spherical cap of radius Rand height h. The relation between height and contact angle is h=R(1 −cos θ).(48) The footprint radius ais a=Rsin θ. Thus the entire geometry is determined by a single parameter θ, which acts as a geometric “setting” imposed by surface energies. 19.3 Wetting ridges on soft substrates For soft solids, the substrate deforms under the droplet due to capillary forces. The wetting ridge height hwsatisfies hw∼σLG E, where Eis the Young’s modulus of the substrate. Wetting ridges illustrate a three-way geometric coupling: •liquid surface tension shapes the droplet, •solid elasticity reshapes the substrate, •adhesion determines the contact angle and footprint. 19.4 Capillary bridge geometries A droplet confined between two plates forms a capillary bridge whose shape satisfies 2σH = ∆p+σadhκcontact,(49) where σadh is an adhesion term at the plates and κcontact is the curvature at the boundary. Common geometries include: •axisymmetric catenoid-like bridges (for zero pressure difference), •bulging or necked bridges with varying curvature, •transitions between convex and concave geometries as plate separation changes. 19.5 Pinned contact lines and asymmetric shapes If the contact line is pinned (cannot move), the droplet may develop: •asymmetric shapes under gravity, •faceted or multi-curvature boundaries, •hysteresis between advancing and receding angles, •“stretched” droplets when the substrate is tilted. Pinned contact lines impose geometric constraints similar to imposed boundary conditions in elastic membranes. 39 19.6 Interpretation Wetting and adhesion geometries reveal how surface interactions generate new shapes by modifying curvature boundary conditions. Key principles are: •Young’s law sets the contact angle, acting as a geometric constraint or “setting”. •Sessile droplets break spherical symmetry through interaction with a substrate. •Capillary bridges represent minimal-area surfaces with imposed boundaries. •Wetting ridges show that geometry can arise from a competition between surface tension and solid elasticity. •All wetting and adhesion geometries require a fluid or soft interface to adapt the shape to the imposed constraints. Thus wetting and adhesion provide another geometry class in the compression–decompression framework, generating shapes where curvature and surface tension interact with solid boundaries to produce non-spherical forms. 20 Shell Buckling and Core–Shell Deformation: Indented Spheres, Dimples, and Pressurized Shells Thin shells under pressure or indentation exhibit a rich variety of geometric behaviors. Unlike membranes (which resist bending weakly) or solid bodies (which resist stretching strongly), shells combine stretching and bending in a highly nonlinear way. This leads to characteristic shapes such as dimples, ridges, collapsed shells, egg-like forms, and multidimple patterns. Shell buckling represents another geometry class where spherical symmetry breaks when external loading or internal pressure violates stability conditions. 20.1 Elastic energy of a thin shell For a shell of thickness hand mid-surface with curvature tensor Cij, the Koiter elastic energy is E=Y h 2(1 −ν2)Z(εijεij) + h2 12 (Cij −C0 ij)2dA, (50) where: •Yis Young’s modulus, •νis Poisson’s ratio, •εij is the in-plane strain, •C0 ij is the preferred curvature (zero for a sphere). Stretching energy dominates when curvature changes significantly; bending energy dominates when curvature deviates locally. 40 20.2 Buckling of a pressurized spherical shell For a spherical shell of radius Rand thickness h, subjected to external pressure p, the critical buckling pressure is: pcrit =2Y p3(1 −ν2)h R2 .(51) When p>pcrit, the sphere undergoes symmetry-breaking and forms a dimple. This is the same dimple mechanism observed in: •pollen grains, •virus capsids under osmotic stress, •microcapsules, •aerosphere collapse, •egg-shell deformation. 20.3 Geometry of a single dimple After buckling, the surface exhibits a localized region of negative Gaussian curvature where stretching is minimized at the cost of bending. The indentation depth δsatisfies approximately δ∼Rp pcrit −11/2 . The dimple resembles the biconcave indentation of vesicles but arises from elastic-shell mechanics rather than membrane bending alone. 20.4 Multi-dimple and ridge formation For stronger compression, multiple buckling sites appear. Interactions between dimples produce: •multi-dimple patterns, •ridged shapes, •faceted shell geometries, •“crumpled” but energetically ordered deformations. Ridges form where bending energy concentrates. Their energy scales as Eridge ∼Y h5/2R−1/2. 41 23 Porous and Sponge Geometries: Percolation Networks, Trabecular Bone, and Reticulated Structures Porous, sponge-like geometries fill volumes with interconnected cavities and channels. These structures arise when competing forces drive the formation of a continuous solid matrix interwoven with fluid or air domains. Such geometries appear in natural sponges, foams, bone trabeculae, porous rocks, polymer networks, filtration membranes, and cellular architectures of tissues. Porous geometries represent a regime where the interface between phases becomes highly convoluted, balancing mechanical support with transport efficiency. 23.1 Percolation networks and critical geometry Percolation theory models porous structures as random occupation of sites or bonds in a lattice. A cluster becomes system-spanning when the occupation probability pexceeds the percolation threshold pc. Near the critical point p=pc, the cluster has a fractal dimension Df≈(91/48 ≈1.896 (2D), 2.52 (3D).(56) These fractal networks resemble natural porous geometries where: •cavities form irregular labyrinths, •surfaces have rough, multiscale curvature, •connectivity emerges spontaneously at critical density. Such porous networks represent an emergent geometry resulting from density-driven compression and the avoidance of collapse. 23.2 Trabecular bone as an optimal porous structure Trabecular bone organizes into a sponge-like network of rods and plates aligned with stress fields. According to Wolff’s law, bone remodels to minimize strain energy: δEstrain = 0,(57) leading to alignment of trabeculae with principal compressive and tensile directions. Trabecular architecture emerges from: •mechanical loading, •anisotropic remodeling, •constrained volume occupation, •metabolic cost minimization. The resulting geometry is a mechanically optimized porous network. 48 23.3 Reticulated polymers and stochastic foams Polymer foams and hydrogels exhibit reticulated geometries formed through: •gas expansion during polymerization, •phase separation, •removal of one phase (solvent extraction), •surface-tension-driven coarsening. The resulting pores are typically: •irregularly shaped, •connected, •multiscale, •organized by energy minimization under isotropic or anisotropic compression. 23.4 Mechanical optimality and geometric scaling Porous structures achieve a balance between stiffness and mass. For a network with volume fraction ϕ, Eeff ∼ϕn, where ndepends on whether struts deform by bending (n≈2) or stretching (n≈1). This scaling determines the optimum density and geometric arrangement of the porous structure. Furthermore, the interface between solid and pore space often obeys surface minimization principles similar to foams: 2σH =psolid −ppore, where the pressures reflect mechanical loading or fluid infiltration. 23.5 Sponge-like biological tissues Many tissues exhibit porous or labyrinthine architectures: •lymphatic sinuses, •cancellous (trabecular) bone, •porous plant parenchyma, •spongy mesophyll in leaves, •reticular connective tissue. They combine: •mechanical support, 49 •fluid or gas transport, •metabolic exchange, •shock absorption, •low material cost. Their geometry results from distributed compression with localized stiffness variations and tension-bearing fibers. 23.6 Interpretation Porous and sponge geometries form through the cooperative action of: •phase separation, •mechanical loading, •surface tension, •percolation phenomena, •anisotropic remodeling. Within the compression–decompression framework: •high compression tends to collapse pores, •surface tension tends to round interfaces, •structural reinforcement preserves voids, •percolation thresholds define connectivity, •transport demands sustain labyrinthine patterns. Thus porous, sponge-like geometries form a distinct class of shapes dictated by the competition between mechanical integrity, minimal material usage, and the need for efficient internal transport. 24 Polygonal Cell Tilings: Voronoi Geometry, Vertex Models, and Epithelial Packing Many biological tissues form two-dimensional tilings composed of polygonal cells. These patterns arise from the balance of intracellular pressure, junctional tension, cell–cell adhesion, and constraints on area and perimeter. The result is a geometry class governed by piecewise-linear boundaries forming polygons whose shapes depend on local mechanical interactions. Polygonal tilings also appear in foams, granular materials, surface networks, and optimization problems, highlighting their generality. 50 24.1 Voronoi geometry as a first approximation A Voronoi tessellation partitions the plane into regions closest to a set of points. If cell centers are at positions {xi}, the region for cell iis Vi={x:|x−xi|<|x−xj|,∀j=i}. Although real tissues are not exact Voronoi diagrams, this model captures: •polygonal shape emergence, •adjacency relations, •average coordination number (≈6), •isotropic packing tendencies. 24.2 Vertex models for epithelial tissues A more accurate representation uses vertex models where each cell is a polygon with vertices {vk}. The energy for the tissue is typically E=X iKA 2(Ai−A0i)2+KP 2(Pi−P0i)2+X edges e ΓeLe,(58) where: •Aiis the cell area, •A0iis its preferred area (pressure constraint), •Piis the cell perimeter, •P0iis preferred perimeter (adhesion–tension balance), •Leis edge length, •Γeis line tension for each edge. Minimization of (58) generates realistic epithelial geometries. 24.3 Topological constraints: Euler relation For a planar tiling with Ffaces, Eedges, and Vvertices, Euler’s formula gives V−E+F= 1. If most vertices are 3-fold coordinated (as in epithelial tissues), one can derive the average number of edges per face: ⟨n⟩= 6. This explains the ubiquity of hexagonal packing in many biological tissues and foams. 51 24.4 Mechanical equilibrium at vertices At each vertex, force balance requires X e Γeˆ te= 0, where ˆ teare tangents to edges meeting at the vertex. For equal line tensions, edges meet at 120◦angles, producing hexagonal tilings—similar to Plateau’s rules for foams. 24.5 Disordered and anisotropic tessellations If tensions or preferred perimeters vary across the tissue, the tiling becomes anisotropic or disordered, producing: •elongated polygons, •regions of high and low coordination number, •anisotropic stress patterns, •neighbor exchanges (T1 transitions), •rosette formations (multiple cells meeting at a point). These phenomena appear during: •embryonic morphogenesis, •tissue stretching or compression, •wound healing, •collective cell migration. 24.6 Geometric transitions and rigidity The vertex model exhibits a rigidity transition controlled by the shape index p=P √A. For p < pc≈3.81, the tissue behaves like a solid; for p > pc, it becomes fluid-like. This provides a geometric criterion for fluid–solid transitions in epithelial sheets. 24.7 Interpretation Polygonal cell tilings show how geometry emerges from the balance of: •pressure (setting the area), •tension (setting the perimeter), •adhesion (setting edge energies), 52 •topological constraints (Euler relation), •mechanical equilibrium at vertices. Key insights: •Hexagonal packing is the “default” when tension is uniform. •Anisotropic tension generates elongated polygons and oriented structures. •Changes in pressure or preferred perimeter drive transitions between shapes. •The tiling becomes rigid or fluid depending on geometric parameters. •Soft, reconfigurable states are required to reach the optimal packing before stabilization. Thus polygonal tilings form another geometry class under the compression–decompression framework, governed by energy minimization and topological constraints in two dimensions. 25 Jamming and Marginal Rigidity: Force Chains, Granular Packing, and Critical Geometries Jammed systems arise when a collection of discrete elements (grains, cells, droplets, or particles) are compressed to a point where they cannot rearrange without deforming. At this threshold, known as the jamming transition, the system becomes mechanically rigid, yet remains geometrically disordered. Jamming introduces a geometry class where shape emerges not from a smooth surface but from the arrangement of discrete bodies constrained by contact forces and global compression. 25.1 Contact network and isostaticity A jammed packing of Nfrictionless particles is isostatic when the number of contacts Z satisfies Z= 2d, (59) where dis the dimensionality (e.g. Z= 6 in 3D, Z= 4 in 2D). At this threshold: •the system is marginally rigid, •small perturbations produce large rearrangements, •force chains develop spontaneously. This condition marks the geometric boundary between fluid-like and solid-like behavior. 53 25.2 Force chains and emergent anisotropy In a jammed system under compression, forces do not distribute uniformly. Instead, they organize into filament-like “force chains” satisfying local balance: X j∈contacts Fij =0. These chains: •carry most of the load, •form anisotropic structures, •define the effective mechanical geometry, •produce patterns reminiscent of tension-web networks (cf. Section 15). 25.3 Packing geometry and Voronoi cells Each grain occupies a Voronoi region whose shape is determined by contact geometry. As compression increases, Voronoi cells: •become more isotropic near jamming, •elongate under shear, •develop angular facets under strong confinement, •form disordered polygonal tilings in 2D and polyhedral packings in 3D. Thus the geometry of jamming interacts with the polygonal tilings of Section 24. 25.4 Critical scaling at the jamming transition Near the jamming density ϕc, key quantities obey scaling laws: Z−Zc∼(ϕ−ϕc)1/2,(60) G∼(ϕ−ϕc)1/2,(61) δ∼(ϕ−ϕc),(62) where: •Zc= 2dis the isostatic coordination, •Gis the shear modulus, •δis particle overlap. These relations show how mechanical rigidity emerges from a purely geometric constraint: the network becomes just constrained enough to resist deformation. 54 25.5 Granular architecture and anisotropic loading When granular systems are subject to anisotropic compression, they develop: •shear bands, •layering, •arches and vaults, •localized regions of high curvature in the contact network. Arches form when force chains percolate into closed loops that support weight without global rearrangement. Such structures illustrate geometry emerging from frictional and contact forces. 25.6 Biological jamming Cells can jam in crowded tissues, forming near-rigid yet disordered architectures. This occurs in: •embryonic epithelia, •tumors, •wound-healing fronts, •confluent monolayers. The cell packing geometry obeys similar principles to granular matter: •contact networks define mechanical rigidity, •transitions between fluid-like and solid-like states correspond to geometric rearrangements, •force chains form through adherens junctions and cytoskeletal tension. 25.7 Interpretation Jamming and marginal rigidity represent a geometry class governed by: •discrete particle contacts, •local force balance, •isostatic constraints, •anisotropic stress redistribution, •emergent force-chain networks. Within the compression–decompression framework: •spherical packing is unstable under crowding, 55 •contact forces break symmetry as density increases, •force chains generate anisotropic geometry, •marginal rigidity corresponds to a critical geometric state, •the system becomes shape-bearing only when sufficiently compressed. Thus jamming adds a final geometry family to the atlas, capturing how shape emerges from compression in systems composed of discrete units rather than continuous surfaces or networks. 56 References [1] W. Helfrich, “Elastic Properties of Lipid Bilayers: Theory and Possible Experiments,” Zeitschrift f¨ur Naturforschung C, vol. 28, pp. 693–703, 1973. [2] P. Canham, “The Minimum Energy of Bending as a Possible Explanation of the Biconcave Shape of the Human Red Blood Cell,” Journal of Theoretical Biology, vol. 26, pp. 61–81, 1970. [3] E. Evans, “Minimum Energy Analysis of Membrane Deformation Applied to Red Blood Cells,” Biophysical Journal, vol. 30, pp. 265–284, 1980. [4] C. D. Murray, “The Physiological Principle of Minimum Work,” PNAS, vol. 12, pp. 207–214, 1926. [5] A. M. Turing, “The Chemical Basis of Morphogenesis,” Philosophical Transactions of the Royal Society B, vol. 237, pp. 37–72, 1952. [6] J. Plateau, Statique Exp´erimentale et Th´eorique des Liquides Soumis aux Seules Forces Mol´eculaires. Paris: Gauthier-Villars, 1873. [7] W. T. Koiter, “On the Stability of Elastic Equilibrium,” NASA Technical Translation, 1966. [8] D. Weaire and R. Phelan, “A Structure for Foam of Minimal Surface Area,” Philosophical Magazine Letters, vol. 69, pp. 107–110, 1994. [9] H. A. Schwarz, “Gesammelte Mathematische Abhandlungen,” vol. 1, 1867 (Schwarz P and D surfaces). [10] S. Gandy et al., “Biological Gyroids: Structural and Optical Properties,” Proceedings of the Royal Society A, 2001. [11] A. J. Liu and S. R. Nagel, “Jamming is Not Just Cool Any More,” Nature, vol. 396, pp. 21–22, 1998. 57