Sphere packing in a cylinder

From HandWiki
Short description: Three-dimensional packing problem
SpherePacking
Illustration of a columnar structure assembled by golf balls.

Sphere packing in a cylinder is a three-dimensional packing problem with the objective of packing a given number of identical spheres inside a cylinder of specified diameter and length. For cylinders with diameters on the same order of magnitude as the spheres, such packings result in what are called columnar structures.

These problems are studied extensively in the context of biology, nanoscience, materials science, and so forth due to the analogous assembly of small particles (like cells and atoms) into cylindrical crystalline structures.

The book "Columnar Structures of Spheres: Fundamentals and Applications"[1] serves as a notable contributions to this field of study. Authored by Winkelmann and Chan, the book reviews theoretical foundations and practical applications of densely packed spheres within cylindrical confinements.

Appearance in science

Columnar structures appear in various research fields on a broad range of length scales from metres down to the nanoscale. On the largest scale, such structures can be found in botany where seeds of a plant assemble around the stem. On a smaller scale bubbles of equal size crystallise to columnar foam structures when confined in a glass tube. In nanoscience such structures can be found in man-made objects which are on length scales from a micron to the nanoscale.

arum maculatum
The berries of the Arum maculatum form a columnar structure (Bushy Park).

Botany

Columnar structures were first studied in botany due to their diverse appearances in plants.[2] D'Arcy Thompson analysed such arrangement of plant parts around the stem in his book "On Growth and Form" (1917). But they are also of interest in other biological areas, including bacteria,[3] viruses,[4] microtubules,[5] and the notochord of the zebra fish.[6]

One of the largest flowers where the berries arrange in a regular cylindrical form is the titan arum. This flower can be up to 3m in height and is natively solely found in western Sumatra and western Java.

On smaller length scales, the berries of the Arum maculatum form a columnar structure in autumn. Its berries are similar to that of the corpse flower, since the titan arum is its larger relative. However, the cuckoo-pint is much smaller in height (height ≈ 20 cm). The berry arrangement varies with the stem to berry size.

Another plant that can be found in many gardens of residential areas is the Australian bottlebrush. It assembles its seed capsules around a branch of the plant. The structure depends on the seed capsule size to branch size.

Foams

Foam
Spherical soap bubbles confined in a cylindrical glass tube.

A further occurrence of ordered columnar arrangement on the macroscale are foam structures confined inside a glass tube. They can be realised experimentally with equal-sized soap bubbles inside a glass tube, produced by blowing air of constant gas flow through a needle dipped in a surfactant solution.[7] By putting the resulting foam column under forced drainage (feeding it with surfactant solution from the top), the foam can be adjusted to either a dry (bubbles shaped as polyhedrons) or wet (spherical bubbles) structure.[8]

Due to this simple experimental set-up, many columnar structures have been discovered and investigated in the context of foams with experiments as well as simulation. Many simulations have been carried out using the Surface Evolver to investigate dry structure or the hard sphere model for the wet limit where the bubbles are spherical.

In the zigzag structure the bubbles are stacked on top of each other in a continuous w-shape. For this particular structure a moving interface with increasing liquid fraction was reported by Hutzler et al. in 1997.[9] This included an unexpected 180° twist interface, whose explanation is still lacking.

The first experimental observation of a line-slip structure was discovered by Winkelmann et al. in a system of bubbles.[10]

Further discovered structures include complex structures with internal spheres/foam cells. Some dry foam structures with interior cells were found to consist of a chain of pentagonal dodecahedra or Kelvin cells in the centre of the tube.[11] For many more arrangements of this type, it was observed that the outside bubble layer is ordered, with each internal layer resembling a different, simpler columnar structure by using X-ray tomography.[7]

Nanoscience

Columnar structures have also been studied intensively in the context of nanotubes. Their physical or chemical properties can be altered by trapping identical particles inside them.[12][13][14] These are usually done by self-assembling fullerenes such as C60, C70, or C78 into carbon nanotubes,[12] but also boron nitride nanotubes[15]

Such structures also assemble when particles are coated on the surface of a spherocylinder as in the context of pharmaceutical research. Lazáro et al. examined the morphologies of virus capsid proteins self-assembled around metal nanorods.[16] Drug particles were coated as densely as possible on a spherocylinder to provide the best medical treatment.

Wu et al. built rods of the size of several microns. These microrods are created by densely packing silica colloidal particles inside cylindrical pores. By solidifying the assembled structures the microrods were imaged and examined using scanning electron microscopy (SEM).[17]

Columnar arrangements are also investigated as a possible candidate of optical metamaterials (i.e. materials with a negative refractive index) which find applications in super lenses[18] or optical cloaking.[19] Tanjeem et al. are constructing such a resonator by self-assembling nanospheres on the surface of the cylinder.[20][21] The nanospheres are suspended in an SDS solution together with a cylinder of diameter [math]\displaystyle{ D }[/math], much larger than the diameter of the nanospheres [math]\displaystyle{ d }[/math] ([math]\displaystyle{ D/d \approx 3 \text{ to } 5 }[/math]). The nanospheres then stick to the surface of the cylinders by a depletion force.

Classification using phyllotactic notation

The most common way of classifying ordered columnar structures uses the phyllotactic notation, adopted from botany. It is used to describe arrangements of leaves of a plant, pine cones, or pineapples, but also planar patterns of florets in a sunflower head. While the arrangement in the former are cylindrical, the spirals in the latter are arranged on a disk. For columnar structures phyllotaxis in the context of cylindrical structures is adopted.

The phyllotactic notation describes such structures by a triplet of positive integers [math]\displaystyle{ (l = m+n, m, n) }[/math] with [math]\displaystyle{ l \geq m \geq n }[/math]. Each number [math]\displaystyle{ l }[/math], [math]\displaystyle{ m }[/math], and [math]\displaystyle{ n }[/math] describes a family of spirals in the 3-dimensional packing. They count the number of spirals in each direction until the spiral repeats. This notation, however, only applies to triangular lattices and is therefore restricted to the ordered structures without internal spheres.

Types of ordered columnar structures without internal spheres

Ordered columnar structures without internal spheres are categorised into two separate classes: uniform and line-slip structures. For each structure that can be identified with the triplet [math]\displaystyle{ (l, m, n) }[/math], there exist a uniform structure and at least one line slip.

Uniform structure

An example uniform structure and its corresponding rolled-out contact network. The identical vicinity of each sphere defines a uniform structure.

A uniform structure is identified by each sphere having the same number of contacting neighbours.[22][1] This gives each sphere an identical neighbourhood. In the example image on the side each sphere has six neighbouring contacts.

The number of contacts is best visualised in the rolled-out contact network. It is created by rolling out the contact network into a plane of height [math]\displaystyle{ z }[/math] and azimuthal angle [math]\displaystyle{ \theta }[/math] of each sphere. For a uniform structure such as the one in the example image, this leads to a regular hexagonal lattice. Each dot in this pattern represents a sphere of the packing and each line a contact between adjacent spheres.

For all uniform structures above a diameter ratio of [math]\displaystyle{ D/d \gt 2.0 }[/math], the regular hexagonal lattice is its characterising feature since this lattice type has the maximum number of contacts.[22][1] For different uniform structures [math]\displaystyle{ (l, m, n) }[/math] the rolled-out contact pattern only varies by a rotation in the [math]\displaystyle{ z \text{-}\theta }[/math] plane. Each uniform structure is thus distinguished by its periodicity vector [math]\displaystyle{ V }[/math], which is defined by the phyllotactic triplet [math]\displaystyle{ (l, m, n) }[/math].

Line-slip structure

An example line-slip structure and its corresponding rolled-out contact network. A line slip is identified by the loss of contacts.

For each uniform structure, there also exists a related but different structure, called a line-slip arrangement.[22][1]

The differences between uniform and line-slip structures are marginal and difficult to spot from images of the sphere packings. However, by comparing their rolled-out contact networks, one can spot that certain lines (which represent contacts) are missing.

All spheres in a uniform structure have the same number of contacts, but the number of contacts for spheres in a line slip may differ from sphere to sphere. For the example line slip in the image on the right side, some spheres count five and others six contacts. Thus a line slip structure is characterised by these gaps or loss of contacts.

Such a structure is termed line slip because the losses of contacts occur along a line in the rolled-out contact network. It was first identified by Picket et al., but not termed line slip.[23]

The direction, in which the loss of contacts occur can be denoted in the phyllotactic notation [math]\displaystyle{ (l, m, n) }[/math], since each number represents one of the lattice vectors in the hexagonal lattice.[22][1] This is usually indicated by a bold number.

By shearing the row of spheres below the loss of contact against a row above the loss of contact, one can regenerate two uniform structures related to this line slip. Thus, each line slip is related to two adjacent uniform structures, one at a higher and one at a lower diameter ratio [math]\displaystyle{ D/d }[/math].[22][1][24]

Winkelmann et al. were the first to experimentally realise such a structure using soap bubbles in a system of deformable spheres.[10]

Dense sphere packings in cylinders

File:PackingRatio.pdf Columnar structures arise naturally in the context of dense hard sphere packings inside a cylinder. Mughal et al. studied such packings using simulated annealing up to the diameter ratio of [math]\displaystyle{ D/d = 2.873 }[/math] for cylinder diameter [math]\displaystyle{ D }[/math] to sphere diameter [math]\displaystyle{ d }[/math].[24] This includes some structures with internal spheres that are not in contact with the cylinder wall.

They calculated the packing fraction for all these structures as a function of the diameter ratio. At the peaks of this curve lie the uniform structures. In-between these discrete diameter ratios are the line slips at a lower packing density. Their packing fraction is significantly smaller than that of an unconfined lattice packing such as fcc, bcc, or hcp due to the free volume left by the cylindrical confinement.

The rich variety of such ordered structures can also be obtained by sequential depositioning the spheres into the cylinder.[25] Chan reproduced all dense sphere packings up to [math]\displaystyle{ D/d \lt 2.7013 }[/math] using an algorithm, in which the spheres are placed sequentially dropped inside the cylinder.

Mughal et al. also discovered that such structures can be related to disk packings on a surface of a cylinder.[24] The contact network of both packings are identical. For both packing types, it was found that different uniform structures are connected with each other by line slips.[24]

Fu et al. extended this work to higher diameter ratios [math]\displaystyle{ D/d \lt 4.0 }[/math] using linear programming and discovered 17 new dense structures with internal spheres that are not in contact with the cylinder wall.[26]

A similar variety of dense crystalline structures have also been discovered for columnar packings of spheroids through Monte Carlo simulations.[27] Such packings include achiral structures with specific spheroid orientations and chiral helical structures with rotating spheroid orientations.

Columnar structures created by rapid rotations

Columnar structures are assembled by using rapid rotations around a central axis to drive the spheres towards this axis.

A further dynamic method to assemble such structures was introduced by Lee et al.[28] Here, polymeric beads are placed together with a fluid of higher density inside a rotating lathe.

When the lathe is static, the beads float on top of the liquid. With increasing rotational speed, the centripetal force then pushes the fluid outwards and the beads toward the central axis. Hence, the beads are essentially confined by a potential given by the rotational energy[math]\displaystyle{ E_{\text{rot}} = \frac{1}{2} m R^2 \omega^2, }[/math]where [math]\displaystyle{ m }[/math] is the mass of the beads, [math]\displaystyle{ R }[/math] the distance from the central axis, and [math]\displaystyle{ \omega }[/math] the rotational speed. Due to the [math]\displaystyle{ R^2 }[/math] proportionality, the confining potential resembles that of a cylindrical harmonic oscillator.

Depending on number of spheres and rotational speed, a variety of ordered structures that are comparable to the dense sphere packings were discovered.

A comprehensive theory to this experiment was developed by Winkelmann et al.[29] It is based on analytic energy calculations using a generic sphere model and predicts peritectoid structure transitions.

See also

References

  1. 1.0 1.1 1.2 1.3 1.4 1.5 Chan, Jens Winkelmann, Ho-Kei (2023-03-31). Columnar Structures of Spheres: Fundamentals and Applications. New York: Jenny Stanford Publishing. doi:10.1201/9780429092114. ISBN 978-0-429-09211-4. https://www.taylorfrancis.com/books/edit/10.1201/9780429092114/columnar-structures-spheres-ho-kei-chan-jens-winkelmann. 
  2. Erickson, R. O. (1973-08-24). "Tubular Packing of Spheres in Biological Fine Structure" (in en). Science 181 (4101): 705–716. doi:10.1126/science.181.4101.705. ISSN 0036-8075. PMID 4579682. Bibcode1973Sci...181..705E. https://www.science.org/doi/10.1126/science.181.4101.705. 
  3. Amir, A.; Nelson, D. R. (2012-06-19). "Dislocation-mediated growth of bacterial cell walls" (in en). Proceedings of the National Academy of Sciences 109 (25): 9833–9838. doi:10.1073/pnas.1207105109. ISSN 0027-8424. PMID 22660931. Bibcode2012PNAS..109.9833A. 
  4. Hull, Roger (1976-01-01), "The Structure of Tubular Viruses", in Lauffer, Max A.; Bang, Frederik B.; Maramorosch, Karl et al. (in en), Advances in Virus Research Volume 20, Advances in Virus Research, 20, Academic Press, pp. 1–32, doi:10.1016/s0065-3527(08)60500-x, ISBN 9780120398201, PMID 775945, http://www.sciencedirect.com/science/article/pii/S006535270860500X, retrieved 2020-04-13 
  5. Bryan, Joseph (1974-12-01). "Microtubules" (in en). BioScience 24 (12): 701–711. doi:10.2307/1297089. ISSN 0006-3568. https://academic.oup.com/bioscience/article/24/12/701/241436. 
  6. Norman, James; Sorrell, Emma L.; Hu, Yi; Siripurapu, Vaishnavi; Garcia, Jamie; Bagwell, Jennifer; Charbonneau, Patrick; Lubkin, Sharon R. et al. (2018-11-05). "Tissue self-organization underlies morphogenesis of the notochord". Philosophical Transactions of the Royal Society B: Biological Sciences 373 (1759): 20170320. doi:10.1098/rstb.2017.0320. PMID 30249771. 
  7. 7.0 7.1 Meagher, A. J.; García-Moreno, F.; Banhart, J.; Mughal, A.; Hutzler, S. (2015-05-20). "An experimental study of columnar crystals using monodisperse microbubbles" (in en). Colloids and Surfaces A: Physicochemical and Engineering Aspects. A Collection of Papers Presented at the 10th Eufoam Conference, Thessaloniki, Greece,7-10 July 2014 473: 55–59. doi:10.1016/j.colsurfa.2014.12.020. ISSN 0927-7757. http://www.sciencedirect.com/science/article/pii/S092777571400939X. 
  8. Weaire, D.; Hutzler, S.; Verbist, G.; Peters, E. (2007), "A Review of Foam Drainage" (in en), Advances in Chemical Physics (John Wiley & Sons, Ltd): pp. 315–374, doi:10.1002/9780470141618.ch5, ISBN 978-0-470-14161-8, https://onlinelibrary.wiley.com/doi/abs/10.1002/9780470141618.ch5, retrieved 2020-04-16 
  9. Hutzler, S.; Weaire, D.; Crawford, R. (1997-06-01). "Moving boundaries in ordered cylindrical foam structures". Philosophical Magazine B 75 (6): 845–857. doi:10.1080/13642819708205711. ISSN 1364-2812. Bibcode1997PMagB..75..845H. https://doi.org/10.1080/13642819708205711. 
  10. 10.0 10.1 Winkelmann, J.; Haffner, B.; Weaire, D.; Mughal, A.; Hutzler, S. (2017-07-31). "Corrected Article: Simulation and observation of line-slip structures in columnar structures of soft spheres [Phys. Rev. E 96, 012610 (2017)]". Physical Review E 97 (5): 059902. doi:10.1103/PhysRevE.97.059902. PMID 29906839. Bibcode2017PhRvE..97e9902W. 
  11. Saadatfar, M.; Barry, J.; Weaire, D.; Hutzler, S. (2008-09-01). "Ordered cylindrical foam structures with internal bubbles". Philosophical Magazine Letters 88 (9–10): 661–668. doi:10.1080/09500830802307658. ISSN 0950-0839. Bibcode2008PMagL..88..661S. https://doi.org/10.1080/09500830802307658. 
  12. 12.0 12.1 Troche, Karla S.; Coluci, Vitor R.; Braga, Scheila F.; Chinellato, David D.; Sato, Fernando; Legoas, Sergio B.; Rurali, Riccardo; Galvão, Douglas S. (2005-02-01). "Prediction of Ordered Phases of Encapsulated C60, C70, and C78 Inside Carbon Nanotubes". Nano Letters 5 (2): 349–355. doi:10.1021/nl047930r. ISSN 1530-6984. PMID 15794624. Bibcode2005NanoL...5..349T. https://doi.org/10.1021/nl047930r. 
  13. Sanwaria, Sunita; Horechyy, Andriy; Wolf, Daniel; Chu, Che-Yi; Chen, Hsin-Lung; Formanek, Petr; Stamm, Manfred; Srivastava, Rajiv et al. (2014). "Helical Packing of Nanoparticles Confined in Cylindrical Domains of a Self-Assembled Block Copolymer Structure" (in en). Angewandte Chemie International Edition 53 (34): 9090–9093. doi:10.1002/anie.201403565. ISSN 1521-3773. PMID 24989683. https://onlinelibrary.wiley.com/doi/abs/10.1002/anie.201403565. 
  14. Yamazaki, T; Kuramochi, K; Takagi, D; Homma, Y; Nishimura, F; Hori, N; Watanabe, K; Suzuki, S et al. (2008-01-30). "Ordered fullerene nanocylinders in large-diameter carbon nanotubes". Nanotechnology 19 (4): 045702. doi:10.1088/0957-4484/19/04/045702. ISSN 0957-4484. PMID 21817519. Bibcode2008Nanot..19d5702Y. https://iopscience.iop.org/article/10.1088/0957-4484/19/04/045702. 
  15. Chopra, Nasreen G.; Luyken, R. J.; Cherrey, K.; Crespi, Vincent H.; Cohen, Marvin L.; Louie, Steven G.; Zettl, A. (1995-08-18). "Boron Nitride Nanotubes" (in en). Science 269 (5226): 966–967. doi:10.1126/science.269.5226.966. ISSN 0036-8075. PMID 17807732. Bibcode1995Sci...269..966C. https://www.science.org/doi/10.1126/science.269.5226.966. 
  16. Lázaro, Guillermo R.; Dragnea, Bogdan; Hagan, Michael F. (2018-07-18). "Self-assembly of convex particles on spherocylindrical surfaces" (in en). Soft Matter 14 (28): 5728–5740. doi:10.1039/C8SM00129D. ISSN 1744-6848. PMID 29796568. Bibcode2018SMat...14.5728L. 
  17. Wu, Gaoxiang; Cho, Hyesung; Wood, Derek A.; Dinsmore, Anthony D.; Yang, Shu (2017-04-12). "Confined Assemblies of Colloidal Particles with Soft Repulsive Interactions". Journal of the American Chemical Society 139 (14): 5095–5101. doi:10.1021/jacs.6b12975. ISSN 0002-7863. PMID 28367624. https://doi.org/10.1021/jacs.6b12975. 
  18. Ozbay, Ekmel (2008-11-01). "The Magical World of Photonic Metamaterials" (in EN). Optics and Photonics News 19 (11): 22–27. doi:10.1364/OPN.19.11.000022. ISSN 1541-3721. https://www.osapublishing.org/opn/abstract.cfm?uri=opn-19-11-22. 
  19. Cai, Wenshan; Chettiar, Uday K.; Kildishev, Alexander V.; Shalaev, Vladimir M. (April 2007). "Optical cloaking with metamaterials" (in en). Nature Photonics 1 (4): 224–227. doi:10.1038/nphoton.2007.28. ISSN 1749-4893. Bibcode2007NaPho...1..224C. https://www.nature.com/articles/nphoton.2007.28. 
  20. Manoharan, Vinothan. "Self-assembled metamaterial". https://manoharan.seas.harvard.edu/self-assembled-metamaterials. 
  21. Tanjeem, Nabila. "Tanjeem's website". https://scholar.harvard.edu/nabilatanjeem/research. 
  22. 22.0 22.1 22.2 22.3 22.4 Winkelmann, Jens (2020). Structures of columnar packings with soft and hard spheres (Thesis thesis). Trinity College Dublin. School of Physics. Discipline of Physics. hdl:2262/91733.
  23. Pickett, Galen T.; Gross, Mark; Okuyama, Hiroko (2000-10-23). "Spontaneous Chirality in Simple Systems". Physical Review Letters 85 (17): 3652–3655. doi:10.1103/PhysRevLett.85.3652. PMID 11030973. Bibcode2000PhRvL..85.3652P. https://link.aps.org/doi/10.1103/PhysRevLett.85.3652. 
  24. 24.0 24.1 24.2 24.3 Mughal, A.; Chan, H. K.; Weaire, D.; Hutzler, S. (2012-05-11). "Dense packings of spheres in cylinders: Simulations". Physical Review E 85 (5): 051305. doi:10.1103/PhysRevE.85.051305. PMID 23004748. Bibcode2012PhRvE..85e1305M. https://link.aps.org/doi/10.1103/PhysRevE.85.051305. 
  25. Chan, Ho-Kei (2011-11-14). "Densest columnar structures of hard spheres from sequential deposition". Physical Review E 84 (5): 050302. doi:10.1103/PhysRevE.84.050302. PMID 22181358. Bibcode2011PhRvE..84e0302C. https://link.aps.org/doi/10.1103/PhysRevE.84.050302. 
  26. Fu, Lin; Steinhardt, William; Zhao, Hao; Socolar, Joshua E. S.; Charbonneau, Patrick (2016-02-23). "Hard sphere packings within cylinders" (in en). Soft Matter 12 (9): 2505–2514. doi:10.1039/C5SM02875B. ISSN 1744-6848. PMID 26843132. Bibcode2016SMat...12.2505F. https://pubs.rsc.org/en/content/articlelanding/2016/sm/c5sm02875b. 
  27. Jin, Weiwei; Chan, Ho-Kei; Zhong, Zheng (2020-06-16). "Shape-Anisotropy-Induced Ordered Packings in Cylindrical Confinement". Physical Review Letters 124 (24): 248002. doi:10.1103/PhysRevLett.124.248002. PMID 32639829. Bibcode2020PhRvL.124x8002J. 
  28. Lee, Taehoon; Gizynski, Konrad; Grzybowski, Bartosz A. (2017). "Non-Equilibrium Self-Assembly of Monocomponent and Multicomponent Tubular Structures in Rotating Fluids" (in en). Advanced Materials 29 (47): 1704274. doi:10.1002/adma.201704274. ISSN 1521-4095. PMID 29112327. https://onlinelibrary.wiley.com/doi/abs/10.1002/adma.201704274. 
  29. Winkelmann, J.; Mughal, A.; Williams, D. B.; Weaire, D.; Hutzler, S. (2019-02-25). "Theory of rotational columnar structures of soft spheres". Physical Review E 99 (2): 020602. doi:10.1103/PhysRevE.99.020602. PMID 30934268. Bibcode2019PhRvE..99b0602W. https://link.aps.org/doi/10.1103/PhysRevE.99.020602. 

External links