Wolfram Language

Compute Lattice Packing Densities

The "Lattice" entity domain contains useful information about classical named lattices.

Copy to clipboard.
In[1]:=1
Out[1]=1
Copy to clipboard.
In[2]:=2
Out[2]=2

For example, consider the body-centered cubic (BCC) lattice, which can be immediately visualized using the "Image" property.

Copy to clipboard.
In[3]:=3
Out[3]=3

An important property of a lattice packing is the fraction of space occupied by spheres packed with the relevant configuration. For a BCC lattice, this packing density is given by the following.

Copy to clipboard.
In[4]:=4
Out[4]=4

In other words, spheres fill about 68% of the total space in body-centered cubic packing.

Copy to clipboard.
In[5]:=5
Out[5]=5

The minimal vectors generating this lattice are available as the "MinimalVectors" property.

Copy to clipboard.
In[6]:=6
Out[6]=6

Given these, you can make a space-filling visualization from which this density can be directly computed. Begin by constructing a list giving the centers of spheres appearing in the unit cell from the minimal vectors.

Copy to clipboard.
In[7]:=7
Copy to clipboard.
In[8]:=8
Out[8]=8

Now restrict the plot range to the region of the unit cell.

Copy to clipboard.
In[9]:=9
Out[9]=9

The visualization can be made slightly more elaborate by using filled balls instead of hollow ones.

Copy to clipboard.
In[10]:=10
Copy to clipboard.
In[11]:=11
Out[11]=11

As can be seen from this diagram, there are eight eighth-spheres and one full sphere in the unit cell, contributing full spheres. Taking as the radius of the spheres, by the Pythagorean theorem applied to a space diagonal, the overall cube has edge length , meaning the packing density of spheres is given by the expression already seen.

Copy to clipboard.
In[12]:=12
Out[12]=12

Now consider the face-centered cubic (FCC) lattice.

Copy to clipboard.
In[13]:=13
Out[13]=13

This lattice packing has a higher packing density.

Copy to clipboard.
In[14]:=14
Out[14]=14

That is, spheres therefore fill about 74% of the total space in face-centered cubic packing.

Copy to clipboard.
In[15]:=15
Out[15]=15

As before, make a space-filling visualization by first finding the centers of spheres in an FCC packing.

Copy to clipboard.
In[16]:=16
Out[16]=16
Copy to clipboard.
In[17]:=17
Out[17]=17
Copy to clipboard.
In[18]:=18
Out[18]=18

As can be seen from this diagram, there are eight eighth-spheres and six half-spheres in the unit cell, contributing full spheres. Taking as the radius of the spheres, by the Pythagorean theorem applied to a face diagonal, the overall cube has edge length , meaning the packing density of spheres is given by the expression already retrieved.

Copy to clipboard.
In[19]:=19
Out[19]=19

Hexagonal close packing (HCP) is another closely related packing.

Copy to clipboard.
In[20]:=20
Out[20]=20

It has the same packing density as face-centered cubic.

Copy to clipboard.
In[21]:=21
Out[21]=21

As in face-centered cubic close packing, each sphere is surrounded by 12 other spheres in HCP. Taking the Delaunay mesh of the minimal vectors for this packing is equivalent to connecting the centers of the external 12 spheres and allowing the edges to be the faces of an enclosing polyhedron.

Copy to clipboard.
In[22]:=22
Out[22]=22

This results in the polyhedron known as the triangular orthobicupola.

Copy to clipboard.
In[23]:=23
Out[23]=23

Related Examples