# Packing regular tetrahedra of edge length 1 with a vertex at the origin in a unit sphere

I'm currently venturing through Paul Sally's Fundamentals of Mathematical Analysis. This is an unusual textbook in terms of the difficulty of exercises. I've already been stunned by the very first one:

How many regular tetrahedra of edge length 1 can be packed inside a unit sphere with each one has a vertex located at the origin?

My instinct told me this is highly nontrivial, as I've struggled with it for several weeks. Could somebody give a hint on how to tackle it? Thanks in advance.

Edit: As jmerry and Erik Parkinson shows, the solution has been squeezed between $$20$$ and $$22$$. Can further progress be made?

• Clearly the answer is at least $8$, as you can let one tetrahedron lie comfortably within each octant. More than that I have no idea though. Can you get two tetrahedra in each octant? I don't remember enough about the angles of regular tetrahedra to answer that. Jan 4, 2019 at 7:45
• You can surely put six in a hemisphere, so at least twelve in a sphere. Jan 4, 2019 at 7:51
• What section of the book is this exercise in? Jan 5, 2019 at 2:50
• @Cheerful Parsnip The very first section...And I noticed this is classified as challenge problem rather than exercise. I guess it is set to stump students. Jan 5, 2019 at 3:09
• I have an idea. If you take a regular tetrahedron with edge length 1, and one vertex at origin, you have an triangular pyramid with solid angle about $0.551$ steradians. If we use the largest right circular cone that fits within the pyramid (intersection with the outer face of the tetrahedron is the inscribed circle of the triangle), we might be able to show 20 is the optimum because you cannot fit 21 or 22 points with at least the cone apex angle minimum separation between them. I think. I'm not absolutely sure, though; I'm not good with proofs. Jan 8, 2019 at 4:17

To complement Erik Parkinson's construction and lower bound, here's a rather tight upper bound: $$22$$.

Suppose we place one of these tetrahedrons in the sphere with one vertex at the center and three on the sphere. What is the area of the spherical triangle generated from the three vertices on the sphere and great-circle edges?

The area of a spherical polygon (in the unit sphere) is equal to the excess of the sum of interior angles over the standard for a plane polygon. For example, a triangle with one vertex at the north pole and two vertices $$90^\circ$$ apart on the equator (one face of a regular octahedron) has all three angles equal to $$90^\circ$$ (or equivalently $$\frac{\pi}{2}$$) for a total of $$270^\circ$$. That's an excess of $$\frac{\pi}{2}$$, exactly one eight of the total area of $$4\pi=720^\circ$$.

For the triangles we care about, I'll use coordinates. The first vertex $$O$$ goes at the center $$(0,0,0)$$. Place the second and third on the equator, at latitudes $$\pm 30^\circ$$: $$A=(\frac{\sqrt{3}}{2},\frac12,0)$$ and $$B=(\frac{\sqrt{3}}{2},-\frac12,0)$$.
The third vertex $$C$$ is then directly above the center of $$OAB$$, for $$x$$ coordinate $$\frac{\sqrt{3}}{3}$$ and $$y$$ coordinate zero. Since it's on the sphere, we solve the sphere's equation $$x^2+y^2+z^2=1$$ to get $$C=(\frac{\sqrt{3}}{3},0,\frac{\sqrt{6}}{3})$$.

Now, we need to calculate some angles on the sphere. For that we'll need to calculate some tangent vectors, and equations for the planes they lie in.
The equation for the tangent plane to the sphere at a point $$(u,v,w)$$ is simply $$ux+vy+wz=1$$ - for a smooth surface defined as a level surface $$f(x,y,z)=k$$, the gradient of $$f$$ is perpendicular to the tangent plane. At $$A$$, this equation becomes $$\frac{\sqrt{3}}{2}x +\frac12 y=1$$.
The great circle $$AB$$ lies in the plane $$OAB$$, which has equation $$z=0$$.
The great circle $$AC$$ lies in the plane $$OAC$$. That plane's equation must have the form $$ux+vy+wz=0$$ for some $$u,v,w$$. Since it includes $$A$$ and $$C$$, we get $$\frac{\sqrt{3}}{2}u+\frac12 v = 0 = \frac{\sqrt{3}}{3}u + \frac{\sqrt{6}}{3}w$$. Solving this (and normalizing), we get $$\frac{\sqrt{2}}{3}x-\frac{\sqrt{6}}{3}y-\frac13z=0$$.

The tangent to the great circle $$AB$$ at $$A$$ lies in both the tangent plane $$\frac{\sqrt{3}}{2}x +\frac12 y=1$$ and the great circle's plane $$z=0$$. Therefore, it's perpendicular to both normals and we can find it with a cross product: $$\left(\frac{\sqrt{3}}{2},\frac12,0\right)\times (0,0,1) = \left(\frac12,-\frac{\sqrt{3}}{2},0\right)$$.
Similarly, we have two planes for the tangent to the great circle $$AC$$ at $$A$$, and we get it with the cross product $$\left(\frac{\sqrt{2}}{3},-\frac{\sqrt{6}}{3},-\frac13\right)\times \left(\frac{\sqrt{3}}{2},\frac12,0\right) = \left(\frac16,-\frac{\sqrt{3}}{6},\frac{4\sqrt{2}}{6}\right)$$.

Finally, the angle on the sphere is the angle between these two tangent vectors. They're both normalized to length $$1$$ already, so we just take their dot product $$\frac1{12}+\frac{3}{12}+\frac{0}{12}=\frac13$$ to find the angle $$\arccos\left(\frac13\right)$$.

Numerically, $$\arccos\left(\frac13\right)$$ is approximately $$1.23$$, or $$70.5^\circ$$. By symmetry, the other two angles of the spherical triangle are the same, so the excess is $$3\arccos\left(\frac13\right)-\pi\approx 0.551$$, or about $$31.6^\circ$$. Dividing that into the total area $$4\pi$$ or $$720^\circ$$, we get about $$22.8$$. Since the number of tetrahedra must be an integer, we round down to get that at most $$22$$ tetrahedra of this type can possibly fit in the sphere. It took a lot of space to write down, but all of this is very routine.

And now, that leaves determining whether $$21$$ or $$22$$ tetrahedra are possible. I expect that part to be harder than everything that has come before.

• You mean this is actually an open problem? Jan 4, 2019 at 9:40
• @YuiTo Cheng: I don't know whether this is an open problem. I suspect someone, somewhere, has solved it. There's certainly no full solution in the answers here so far. Jan 4, 2019 at 9:42
• @jmerry: Just to verify, I integrated the solid angle distended by a face when viewed from the opposite vertex in a regular tetrahedron with edge length 1. (To ensure a different approach, I used Monte Carlo integration with rays with uniform direction distribution, and about $10^{10}$ rays, about $10^9$ hitting the face.) The solid angle is indeed about $0.551$ steradians (full solid angle being $4\pi \approx 12.566$ steradians), meaning 22.8 such tetrahedra cover the full solid angle. In short, numerical physical/non-mathematical simulation completely agrees: 22 is the strict upper limit. Jan 8, 2019 at 4:06
• When I spoke to Tom Hales (who proved the optimality of the best sphere packing in 3-dimensions) about this problem 5 years ago, he believed that this was open. Apr 27 at 5:28

Just to illustrate how the packing might work, here are five regular tetrahedra with a common edge (and therefore a common vertex). There is a very narrow angular gap between the two tetrahedra that don't have two adjacent neighbors; the dihedral angle between the faces on either side of that gap is about $$7.36$$ degrees. When you put another $$15$$ tetrahedra around one vertex in this fashion, you get a lot more of these gaps and some larger gaps. The question then is whether there's enough "wiggle room" in the arrangement so that you can make room between the tetrahedra to fit one more in. It doesn't look likely, but proving that is another matter.

• I added an "extended comment" to describe the same via comparing the icosahedral net (with all icosahedron vertices on the unit sphere) and twenty nets of equilateral triangles with unit edge length (corresponding to the outer faces of the tetrahedra). I was trying to find a good solid angle argument, but couldn't find any... Jan 8, 2019 at 23:44

This is not an answer, but an extended comment.

The edge length of an unit icosahedron (all vertices at unit distance from origin) is $$a = \frac{1}{\sin\frac{2 \pi}{5}} \approx 1.0514622$$ according to the Wikipedia Regular icosahedron article. If we compare the net of such an icosahedron to twenty equilateral triangles with edge length 1 (corresponding to the faces of the tetrahedra opposite to the vertex at origin), it is rather obvious that there is not enough room to squeeze in another tetrahedron in there. (This corresponds to actually making 20 unit tetrahedra, and trying to squeeze in an additional tetrahedron while keeping them all within an unit sphere. Remember, the above are nets, not a solid angle diagram.)

My math-fu is not advanced enough to express that as a convincing mathematical argument, though.

It might be interesting to note that R. H. Hardin, N. J. A. Sloane and W. D. Smith list the putative covering radii of points on sphere as (excerpt) $$\begin{array}{c|c} \text{points} & \text{degrees} \\ \hline 19 & 30.3686748 \\ 20 & 29.6230958 \\ 21 & 28.8244768 \\ 22 & 27.8100588 \\ 23 & 27.4818687 \\ \end{array}$$ This means that if we take right circular cones with apex angle say $$2 \cdot 29.6230958 = 59.2461916$$ degrees, up to twenty of them can be connected by their apexes; if the angle is say $$2 \cdot 30.3686748 = 60.7373496$$ degrees, only 19.

Unfortunately, the incircle and circumcircle of the outer face of the tetrahedron differ too much to allow a stricter bound (than 20 to 22 tetrahedra).

• The approach of visualizing via a net is interesting, but you ought to center the new triangles within each cell, and then also keep in mind that the new vertices for the regular tetrahedra are going to be closer to the origin, so they should be a bit smaller even. Given that the inscribed sphere radius for icosahedron (distance from center to midface) is 0.75576... edge lengths, and the height of a tetrahedron is 0.81650... edge lengths, this means if you want length-1 edges for your tetrahedron then the icosahedron edges should have length of 1.08046... Jan 6, 2021 at 5:54

The optimal is 20, forming an icosahedron. This is clearly the max based on the fact that it pacts perfectly around the origin. Not a very satisfying answer mathematically though.

• No, it doesn't pack perfectly around the origin - or rather, those aren't regular tetrahedra. If we take a face of the regular icosahedron and connect it to the center, the face angles at the central vertex are each $72^\circ$ rather than $60^\circ$. We can shrink them to regular tetrahedra and fit 20 that way, but there will definitely be gaps. Jan 4, 2019 at 8:16
• The origin angles are not $72°$ either. In a regular icosahedron draw a fivefold axis. Let A be either vertex through which is passes, and B, C, D, E, F be the adjacent vertices to A in rotational order about the axis. The axis intersects plane BCD (which of course also contains E and F) at point G. It is at this point, not the center, where angle BGC measures $72°$. But yes, the angle at the center is still greater than $60°$. Jan 4, 2019 at 10:44
• To the nearest tenth the central angle actually measures $63.4°$. Its cosine is the square root of $1/5$. Jan 4, 2019 at 10:52