First an easy non dense packing. Suppose we have ball centers at <i, j, k> for arbitrary integral values of i, j and k. They touch if r = 1/2. A unit cell has volume 1 and holds just one ball. The density is thus 4/3 π r^{3} = π/6 = 0.524. Note that each ball touches 6 others.
Suppose that we start with a 2D square array of balls in a plane: <i, j, 0>. We can pile another array closer to the first if we offset the new array thus: <i+1/2, j+1/2, z>. If we set z = 1/√2 then the distance between neighboring centers in the original plane is the same as between interplanar neighbors. r = 1/2 and the cell is [0,1] × [0,1] × [0,√2] whose volume is √2. The cell holds two balls. The density is (π/3)/√2 = 0.740. Each ball touches 12 others.
Lets try two populations of centers, one at <i, j, k> and the other at <i+1/2, j+1/2, k+1/2>. The distance between neighboring centers is now √(3/4). There are two balls per cell; the cell size is 1 and the density is thus 0.680. Each ball touches just 8 others. This seems to be body centered.
Now we try “face centered” as follows:
four populations of centers:
Above are the obvious arrangements starting with perpendicular vectors between centers. What of configurations that begin with hexagonal arrangements? There are two kinds of coordinate systems, the familiar orthogonal coordinate systems that we use above, and oblique coordinates which are especially convenient here. We shall begin with oblique coordinates where one must declare the dot product of each pair of coordinate bases. If e_{j} is the jth basis vector then g_{ij} = e_{i} ⋅ e_{j} which is the metric tensor for this coordinate system. The distance between points given in oblique coordinates is √(∑_{i}∑_{j} g_{ij}dx^{i}dx^{j}).
We choose basis vectors that have dot product 1 if they are the same and 1/2 otherwise. Or:
g_{ij} = 

In each of the configurations so far, translation of the entire array by the vector connecting centers of two neighbors, leaves the array unchanged. Here is a packing without that property. When we placed the second layer in the preceding example we had a choice. We placed the second layer so that a ball at <0, 0, 1> in that layer touched each of <0, 0, 0>, <1, 0, 0> and <0, 1, 0> in the first layer. Alternatively we could have placed the second layer at <i+1/3, j+1/3, 1>. That choice places a ball at <1/3, 1/3, 1> which touches each of <1, 1, 0>, <1, 0, 0> and <0, 1, 0> in the first layer. The choice is fresh for each layer. Chemistry often likes to alternate choices which gives hexagonal closest packing. Regardless of these choices the density tends to .740 . As we shall see, a consistent choice is equivalent to cubic closest packing!
Density and number of balls that touch a given ball are properties of packings and do not depend on coordinates. The face centered cubic packing and hexagonal packing with uniform layering choice share these properties. Are they the same packing?
Yes, except for ball radius. In the coordinates used for face centered cubic closest packing we select the following vectors as a basis for a new coordinate system: <1/2, 1/2, 0>, <1/2, 0, 1/2>, <0, 1/2, 1/2>. Each of these vectors corresponds to a vector connecting centers of neighboring balls in the packing. The metric tensor for the coordinates system with this basis is
g_{ij} = 

The other packings can be discriminated either on the basis of density of number of neighbors.
See diamond’s structure for an even less dense packing.
See some sphere packing in higher dimensions.
If n+1 unit radius balls touch each other in n dimensions, then their centers determine a regular simplex whose volume is given by the CayleyMenger determinant. This simplex, and n! others congruent to it, compose and partition a unit cell of all those points each of whose coordinates belong to [0,1]—a parallelepiped. The content of this cell is √(determinant of g).
For the array of balls with radius = 1/2 and centered on points whose oblique coordinates are all integers.
The metric tensor is
g_{ij} = (δ_{ij} + 1)/2 and its determinant seems to be (n+1)2^{−n}.
The density of this packing ρ = (volume of sphere)/(volume of cell)
= (π^{n/2}(1/2)^{n}/(n/2)!)/√((n+1)2^{−n})
= π^{n/2}2^{−n/2}/((n/2)!√(n+1)).
Here are some <n, ρ_{n}> values
0 1.000000e+00  8 8.455650e02  16 2.229534e04  24 9.421750e08 
1 1.000000e+00  9 4.609806e02  17 9.178457e05  
2 9.068997e01  10 2.402823e02  18 3.680773e05  
3 7.404805e01  11 1.201847e02  19 1.439750e05  
4 5.517277e01  12 5.786488e03  20 5.499537e06  
5 3.798813e01  13 2.688948e03  21 2.053627e06  
6 2.441515e01  14 1.208824e03  22 7.504107e07  
7 1.476489e01  15 5.267991e04  23 2.685649e07 (see)  
// Compute the distance from the origin to each // lattice point in the (2b+1)^n box. #include <stdio.h> #include <math.h> #define n 10 #define b 2 #define do(x) for(ndx[x] = b; ndx[x]<=b; ++ndx[x]) static double g[n][n]; static double const eps = 1.e10; int ndx[n]; void pc(char * m){printf(m); for(int j=0; j<n; ++j) printf(" %d", ndx[j]); printf("\n");} int main(){int nc=0; for(int j=0; j<n; ++j) for(int k=0; k<n; ++k) g[j][k] = j==k?1.:0.5; do(0) do(1) do(2) do(3) do(4) do(5) do(6) do(7) do(8) do(9) {double d = 0; for(int i=0; i<n; ++i) for(int j=0; j<n; ++j) d += g[i][j]*ndx[i]*ndx[j]; if(d < 1eps) pc("Ouch "); if(fabs(d  1.)< eps) {pc("Bing "); ++nc;}} printf("%d neighbors\n", nc); return 0;}
There are more at distance 1 from origin at least according to M. Leech. What is the group of this lattice while holding the origin fixed? Any set of n neighbors that touch each other by pairs serve as a coordinate system with the same metric tensor. The transformation matrix should consist solely of integers.
This relates to String theory?