arXiv:0810.5724v1 [cond-mat.soft] 31 Oct 2008 iesos[7 n hr r orglrtseain of tesselations regular no are two above there R generalized and indeed, law [17] – growth dimensions core the its was at recently geometrical useful only understanding three-dimensions, our and In cuts coherent frustration a foams. provide of [16] description Neumann- law von the growth and Mullins [15] theorem decoration the there foams. projec- Plateau-law-abiding The all are make temperature. insight land- we finite offers that the at tions bubble and evolution 4-sided [14], foam symmetric into foams the random near of scape understand- statistics toward step the new weak ing a Our is a approach tetrahedron. of general symmetric more exception the around the maximum with local bubbles, tetrahedral and A area the both determine We A following. the [12] in that polyhedra” polyhedra” bubbles Plateau of “average “isotropic class or of large notion a the of extend properties the study we n )ec dei h odro he ae,adtedi- 2 the is and faces, faces three these of between border angle the is hedral edge each 2) and ol aif w ue:1 ahvre emntsfour terminates vertex each angle 1) an each edges, rules: that two Recall foams satisfy dry three-dimensional would that polyhedra. 11] [10, regular posited Plateau space-filling with Plateau’s of laws incompatibility with the by along promoted light, disorder with the borders borders these topo- bubble of interaction foam the the of and is network the issue ex- of At and complexity [8] logical [9]. probe treatment experimental theoretical have interior haustive direct bulk both 7]. foam’s evaded [6, a long cosmology of and granular details space-time of the to picture Unfortunately, even Edwards’ and the [5], to matter and 4], copolymers [3, from colloids been con- systems: have fuzzy other they applied many as 2] model of [1, to grown array used has broad foams in a Interest mouth-filling in texts. space-and arise are which foams structures souffles, to cream ing 3 2 3 h tr ntodmnin sqaiaieydifferent: qualitatively is dimensions two in story The rmbe obed east ose n shav- and mousse, to metals bread, to beer From n volume and /V hc aif lta’ as apl,i swl known well is it Happily, laws. Plateau’s satisfy which snal osatfrtecaso dodecahedral of class the for constant nearly is eateto hsc n srnm,Uiest fPennsyl of University Astronomy, and Physics of Department 2rearrangements. T2 rjcinto projection rudcraue esuyacaso eua om on foams regular of class a study We curvature. ground oee,teerglrtseain ontebdi Euclide in embed not do tesselations regular these however, elhsregular has cell n hp itiuin u prahas losu ocomme to re us of explain allows maximum also and local approach extend weak Our and distribution. bubbles, shape foam and bulk of classes broad esltoswt symbols with tesselations N V eua esltoso pc r hrceie hog th through characterized are space of tesselations Regular hda 1]wihw ilcl “average call will we which [13] -hedra” fteeclsadfidta pcfi area, specific that find and cells these of δ cos = R 3 − p en reua el ossetwt lta’ as ean We laws. Plateau’s with consistent cells irregular define 1 gnlsides, -gonal ( − 3 1 pr rmteother, the from apart ) A 2 3 { alD oe n adl .Kamien D. Randall and Modes D. Carl π/ peia om nFa Space Flat in Foams Spherical /V p, .I hsletter this In 3. 3 , ttemxmlysmercttaerlbbl htpartici that bubble tetrahedral symmetric maximally the at q 3 } etn tec etx and vertex, each at meeting l aif lta’ asfreulbimfas o gener For foams. equilibrium for laws Plateau’s satisfy all Dtd uy1,2018) 10, July (Dated: orfcsme ta de h eurmn fdihedral of requirement The 2 edge. of as an angles rules at Plateau’s meet satisfy faces not four does honeycomb cubic The ruda de n ec h Schl¨afli is the symbol hence and edge, an around aetresc de n thus must cell and each edges that such follows three it have edges, four terminate must ubro hs ae etn tec etxof vertex each at meeting faces these of number oe nfltsaei h ui oecm.Snethe Since honeycomb. al- cubic tesselation the cubes, such is are only space cells the flat example, in of lowed way By edge. stenme fegso ahfc ftecell, the of face each on edges of number the is and eua esltoso pc 2] ocasf uhtes- such classify Schl¨afli symbols, To the of use [20]. we language space selations geometric the of to the tesselations connection is simple regular language laws a That Plateau’s for allows consider curvature. that the to in language expedient lines a is these in it along foams, progress to of order case In sys- the 19]. lifted to curvature [18, be background tem may a frustration of of introduction type the this by cases, some in that, coordination. natur foam-like dodecahedral their the and Note cells the further. lines lighter origin, I.1 ipeflt(huhnncnoml rjcinof projection non-conformal) (though into flat 120-cell simple the A 1: FIG. S r 3 stenme fclso ae etn rudeach around meeting faces or cells of number the is hc,truhcnoml stereographic conformal, through which, i Schl¨afli symbols eir ai,Piaepi,P 90-36 USA 19104-6396, PA Philadelphia, vania, nsae u eur nfr back- uniform a require but space, an to omsaiiyb dniyn a identifying by stability foam on nt etosrain nfa structure foam on observations cent π/ r etn nec edge. each on meeting mle that implies 3 R p 3 and 4 = akrlnsaepoetdcoe othe to closer projected are lines Darker . ltclycaatrz a characterize alytically { ,q r q, p, q ;fu ue ttogether fit cubes four 3; = r } .Sneec vertex each Since 3. = q hr each where , .Tu esethat see we Thus 3. = ae in pates Regular { al ,q r q, p, p , } where , { one q 4 , sthe is 3 cell, , of e 4 } p . 2 any regular space tesselation of the form p, 3, 3 neces- sarily represents a valid foam. If we divide{ each}p-gonal face into p isoceles triangles, each will have one angle of δ 2π/p and two angles 2 . On a uniformly curved space of constant sectional curvature K, this triangle, ∆, lies in a two-manifold of curvature K. The Gauss-Bonnet theorem relates the excess angle of the triangle to the integrated curvature:

2π 1 1 1 +2 cos− π = KdA (1) p × 2 −3 − Z∆

In the flat limit K = 0 and thus p 5.1043; real foams ≈ FIG. 2: (color online). The foam-like tetrahedral covering are known to have 5.1 sides per face on average [8, 21] of the sphere and its Reuleaux Triangle image under stereo- ∼ as do the average polyhedra [12, 13]. Of course, p must be graphic projection from a vertex. The full family of possible a positive integer and the theory of polytopes [20] limits conformal distortions of this triangle is accessible by direct us to p 5 (we shun the frumious horocycle – 6, 3, 3 ). geometric calculation. It follows≤ that the regular polytopes can be embedded{ } in a three-space with K > 0, i.e. the three-sphere, S3. These “polychora” are the analogs of the two-dimensional that can fit in a two-dimensional foam. If, for instance, tilings of S2 which correspond to the platonic solids. We we project through (0, 0, 1), then the top vertex is pro- choose to study two tesselations of S3, the so-called 5- jected to infinity and the remaining three vertices become √3 1 √3 1 1 cell, 3, 3, 3 and the 120-cell, 5, 3, 3 (Figure 1). The ′ = ( , − ), ( − , − ), (0, ) . Since great cir- { } { } T { 2√2 2√2 2√2 2√2 √2 } 120-cell is of interest because this regular polytope re- cles map to circles or lines, the edges of the tetrahedron quires the least amount of background curvature to real- become arcs of circles and, by symmetry, we have the fa- ize among the family in question. Meanwhile, the 5-cell mous Reuleaux Triangle [25], familiar in two-dimensional is worthwhile because it features tetrahedral bubbles like foams as the ideal corner bubble during T2 topological those seen at “T2” topological rearrangements, where a rearrangements. Were we to project from another point, foam bubble shrinks until it disappears, leaving behind or equivalently, rotate the sphere and project through the a four-fold vertex [9]. Related constructions have been pole, we would find asymmetric foams. Moreover, it is used to show the uniqueness of k-bubble clusters in d not just the central bubble that satisfies Plateau’s laws; dimensions [22]. the entire projection is guaranteed to have all the correct But how can we make contact with physical foams re- angles for an equilibrium foam. The central cell on which siding in actual, flat space? We observe that the con- we focus could exist in any equilibrium R3 foam. straints of Plateau’s laws that forced the form p, 3, 3 R2 { } Working in , we can calculate both the perime- are all constraints on angles, hence mappings from the ter and area of a general triangular bubble in a two- R3 3-sphere to that preserve angles will also preserve dimensional foam. Here Plateau’s rules are greatly sim- Plateau’s laws. We choose here to work with the con- plified: edges must be arcs of circles and must meet at formal stereographic projection, long familiar to cartog- an angle of 2π/3. In order to parameterize an arbitrary raphers [23] where X = x/(1 w), Y = y/(1 w), and − − triangle, we set two of the vertices at ( 1, 0) and (1, 0), Z = z/(1 w) where X,Y,Z are cartesian coordinates and allow the third to float as (x, y), where− x and y are − R3 in the target , and x,y,z,w are the coordinates in both non-negative. It is tedious but straightforward to S3 R4 with x2 + y2 + z2 + w2 = 1. ⊂ calculate both the perimeter and the area as a function In two dimensions, this projection is the standard of x and y. The exact formulae are complicated and the S2 R2 stereographic projection of to and we may project results are better characterized graphically in Figure 3. the p, 3 Platonic solids (the tetrahedron, cube, and The machinery carries through in exactly the same dodecahedron){ } onto the plane. For concreteness, we way in the three dimensional case as well, allowing us will consider the 3, 3 tetrahedron shown in Fig. 2; to produce a continuous family of space foam bubbles. the conformal projection{ } preserves the 2π/3 angle be- In principle, this technique can produce bubble families tween edges and so we will find a two-dimensional cell from each of the 3, 3, 3 , the 4, 3, 3 , and the 5, 3, 3 that satisfies Plateau’s rules [24]. With one vertex with 4-faced, 6-faced,{ and} 12-faced{ bubbles,} respectively.{ } at the pole, (0, 0, 1), the remaining set of vertices is The class of transformations that belong to the confor- 2 √2 1 2 √2 1 2√2 1 3 = ( 3 , 3 , 3 ), ( 3 , 3 , 3 ), (0, 3 , 3 ) mal group of R is restricted to solid body rotations T { − − − − − − } and theyq are all connectedq by arcs of great circles. We and translations, reflections, dilations, and spherical in- have the freedom to choose from which point we project versions; the degrees of freedom available in the stere- S2, and each projection will give us two-dimensional cells ographic projection precisely account for each of these 3

cide with the highest symmetry target bubble, the av- erage polyhedra [12, 13], whose vertices coincide with the vertices of a regular dodecahedron and whose faces are all equivalent pentagonal spherical caps. This high- symmetry target bubble may be obtained by orienting 1 the pre-image 120-cell so that its vertices lie at 2√2 of the permutations of [20] ( 2, 2, 0, 0), ( √5, 1, 1, 1), 2 ± ± 2 1 ± 1± ±1 ± ( φ, φ, φ, φ− ), and ( φ , φ− , φ− , φ− ) and ±1 ± ± ± ± ± ± 2 ± 2 2√2 of the even permutations of ( φ , 1, φ− , 0), 1 ± ±1 ± ( √5, φ, φ− , 0), and ( 2, φ, 1, φ− ) where φ = (1± + √±5)/2± is the Golden mean.± ± ± ± With this characterization of the reference 120-cell in hand, we may begin to apply elements of SO(4) and project down to flat space. In analogy with the example 3 FIG. 3: (color online). Contour plot of A 2 /V as a function calculation shown above for the projection of the tetra- of the triangle coordinate (x,y). The Reuleaux triangle is a hedron from S2 we obtain the equations of the spheres local maximum, surrounded by three deep valleys. in R3 defining the faces of the image dodecahedral bub- ble. Qualitatively, rotating the pre-image 120-cell around one of the w planes results in a bubble that has been possibilities. Generically, parallel translation of the hy- stretched, with a vertex or edge pulled farther away from persphere being projected relative to the target hyper- the others. This stretching also controls the sign of the plane yields translations. Normal translation gives rise face curvature for the bubbles; we consider relatively to dilations. Rotations and reflections of the hypersphere small deformations in order to preserve the positivity of about hyperplanes that lie in the target hyperplane give the curvature for all twelve faces. the solid body transformations. The rest of the rotations Taking the appropriate face and cell intersections, we give the spherical inversions. Since the shape of a bub- numerically determined the volume and surface area of ble or bubble cluster is sensitive only to inversions from the projected bubble. Perturbing away from the maxi- among these possible transformation types, we can ig- mally symmetric bubble by rotating around the w-z plane nore the translational degrees of freedom and focus only by an angle τ, we calculate the non-dimensionalized com- 3 on the out-of-hyperplane rotations. bination A 2 /V . One would expect this quantity to dom- Guided by experiment [8], we choose first to examine inate the bubble’s free energy and, through the Boltz- the most prevalent from among these bubble types, those mann factor, control the likelihood of finding a bubble with twelve faces. As in the two-dimensional example, we of this shape in an arbitrary dry foam. Interestingly, will have the freedom to choose the coordinate orienta- however, we find that to within the errors of our calcula- 3 3 tion of the 120-cell inside S and the ability to choose the tion A 2 /V is essentially constant with respect to τ (see base point of the stereographic projection; either one of Figure 4). This result extends the observations of Hilgen- these choices completely exhausts the degrees of freedom feldt et al. who found that for their idealized Isotropic available, leaving the other to parameterize the resulting Plateau Polyhedra the non-dimensional surface area is family of bubbles. We choose to pin the base point of the largely insensitive to the number of faces [12]. In ad- projection to the north pole of S3, (0, 0, 0, 1), and will ro- dition, the demonstration that even “stretching” foam tate the 120-cell inside S3 to generate different bubbles. bubbles does not change their reduced surface area-to- In so doing it becomes transparent to identify the param- volume ratio serves to support and strengthen the recent eterization of the foam bubble family with the solid body results of Kraynik et al., who found in simulations using rotation taking place on the 120-cell pre-image. Said an- Surface Evolver that the average reduced surface area- other way, we may unambiguously index our final bulk to-volume ratio of native bulk foam bubbles is constant foam bubbles by elements of SO(4). Note that, due to [14]. The implication, then, is that the bubble shape dis- the high symmetry of the pre-image 120-cell, it will suf- tribution in a bulk foam is dominated by the topology of fice to consider the action of the projection on only one the network of Plateau borders, with little contribution cell, since one cell may be carried into the position of any from bubble geometry. other by some element of SO(4). In analogy to the two- Pursuing the tetrahedral target bubbles generated by 3 dimensional case, planar regions in S map to spherical 3, 3, 3 is qualitatively similar to what was done with 3 caps in R . the{ 120-cell.} The unrotated, reference state vertices S3 √15 1 If we view our parameter space as a higher dimen- for the 5-cell in are : (0, 0, 0, 1), 0, 0, 4 , 4 ,  −  sional rotation space, we will need to specify an unro- 5 1 5 1 1 5 1 5 1 5 1 tated, reference state. We choose this state to coin- 0, 6 , 4 3 , 4 , 2 2 , 2 6 , 4 3 , 4 , and  q − q −   q − q − q −  4

along the lines of the average polyhedra [12, 13]. It is a pleasure to acknowledge discussions with B.G. Chen, D.J. Durian, O.L. Halt, J. M. Kikkawa, R.B. Kus- ner, E. Matsumoto, and V. Vitelli. This work was sup- ported by NSF Grant DMR05-47320, a gift from L. J. Bernstein, and a gift from H. H. Coburn.

[1] K. Feitosa, O. L. Halt, R. D. Kamien, and D. J. Durian, Europhys. Lett. 76, 683 (2006). [2] B. Dollet and F. Graner, J. Fluid Mech. 585, 181 (2007). [3] G.M. Grason, B.A. DiDonna, and R.D. Kamien, Phys. FIG. 4: (color online). Dimensionless surface area to volume Rev. Lett. 91, 058304 (2003); G.M. Grason, Phys. Rep. ratios for dodecahedral (lower curve) and tetrahedral bub- 433, 1 (2006). bles (upper curve) distorted by rotations in the w-z plane [4] P. Ziherl and R.D. Kamien, Phys. Rev. Lett. 85, 3528 by τ followed by sterographic projection. Two different ro- (2000). tational starting points are used for the tetrahedra, with [5] S. F. Edwards and R. B. S. Oakeshott, Physica A, 3 157 no significant difference in A 2 /V . Note the near indepen- (3), 1080 (1989). dence of both bubble types’ ratios on the mapping parame- [6] J.R .Gott, III, M. Dickinson, and A.L. Melott, Ap. J. 306 ter, τ, with the exception of the maximum and valley near , 341 (1986); C. Criado and N. Alamo, Class. Quan- tum Grav. 21 2241 (2004). τ = 0 for the tetrahedron. Likewise, the distorted shapes 2 are still highly constrained by Plateau’s laws. For a sphere, [7] J. A. Wheeler, Ann. Phys. (6), 604 (1957). 3 [8] E. B. Matzke, American Journal of Botany, 33, 58 A 2 /V = 6√π 10.635 attains its minimum. ≈ (1945). [9] D. Weaire and S. Hutzler, The Physics of Foams, Claren- don Press (Oxford, 1999). 1 5 1 5 1 5 1 [10] J.A.F. Plateau, Statique Experimentale et Theorique des 2 2 , 2 6 , 4 3 , 4 . Proceeding as before, we − q − q − q −  Liquides Soumis aux Seules Forces Moleculaires (Gau- perturbed the maximally symmetric tetrahedral bubble thier Villiard, Paris, 1873). by rotating an angle τ about the w-z plane. As is the [11] J.E. Taylor, Ann. of Math. 103, 489 (1976). case for the dodecahedral bubbles, we find a largely flat [12] S. Hilgenfeldt, A.M. Kraynik, D.A. Reinelt, and J.M. 3 dependence of A 2 /V on τ (Figure 4). Furthermore, in Sullivan, Europhys. Lett. 67, 484 (2004). agreement with past work [14], we find that the values of [13] M.E. Glicksman, Phil. Mag. 85, 3 (2005). the reduced surface area-to-volume ratio for the tetrahe- [14] A. M. Kraynik, D. A. Reinelt, and F. van Swol, Phys. Rev. Lett. 93, 208301 (2004). dral bubbles are essentially indistinguishable from those [15] F. Bolton and D. Weaire, Phil. Mag. B 63, 795 (1991). of the dodecahedral bubbles, differing by less than one [16] J. von Neumann, in Metal Interfaces, ed. by C. Her- percent on average. It is noteworthy, though, that the 3 ring, pp. 108110 (American Society for Metals, Cleve- tetrahedra are uniformly higher in A 2 /V , consistent with land, 1952); W.W. Mullins, in Metal Surfaces: Struc- the expectation that the closer a cell is to achieving 13.7 ture, Energetics, and Kinetics ed. by W.D. Robertson faces, the more efficiently it will enclose its volume. and N.A. Gjostein, pp 17-66 (American Society for Met- als, Metals Park, Ohio, 1963). Our numerical results in three dimensions also suggest 3 [17] R.D. MacPherson and D.J. Srolovitz, Nature 446, 1053 there is a local maximum of A 2 /V at the maximally sym- (2007). metric Reuleaux tetrahedron. Whether the existence of [18] J. P. Sethna, D. C. Wright, and N. D. Mermin, Phys. this local maximum is an impediment to T2 transitions at Rev. Lett. 51, 467 (1983); J. P. Sethna, Phys. Rev. Lett. finite temperature is a natural question to pursue. Since 51, 2198 (1983). this feature is preserved in two dimensions, where P 2/A [19] C. D. Modes and R. D. Kamien, Phys. Rev. Lett. 99, 77 as shown in Figure 3, it may be possible to do both ex- 235701 (2007); Phys. Rev. E , 041125 (2008). [20] H. S. M. Coxeter, Regular Polytopes, Methuen & Co. Ltd. periment and simulation to address this issue. (London, 1948). We have generated a continuous class of dodecahe- [21] R. Kusner, Proc. R. Soc. Lond. 439, 683 (1992). dral and tetrahedral bubbles via conformal projection of [22] A.M. Amilibia, Asian J. Math. 5, 23 (2001). known polychora and have shown that the specific area [23] Ptolemy, Planisphaerium (ca. 150); see P. Kunitzsch and of these bubbles is roughly constant over our parameter R. Lorch, Maslama’s Notes on Ptolemy’s Planisphaerium space. Our methods have also allowed us to uncover a and Relatex Texts (M¨unchen, Verlag der Bayerischen 3 Akademie der Wissenschaften, 1994). 2 possible local maximum of A /V at the T2-event bubbles [24] A different but similar approach was studied in A. in both two and three dimensions. It remains an inter- Mughal and D. Weaire, arXiv:0801.2474v1 (2008). esting open question as to whether this approach can be [25] F. Reuleaux, The Kinematics of Machinery: Outlines of generalized to study p, 3, 3 polytopes for non-integer p a Theory of Machines, Macmillan & Co. (London 1876). { }