Archimedes' Revenge Arxiv:2108.05195V1 [Math.HO] 6

Mark B. Villarino and Joseph C. Várilly Escuela de Matemática, Universidad de Costa Rica, 10101 San José, Costa Rica

1 Introduction

Some 2300 years ago the great Archimedes proved that the volume of the orthogonal intersection of two circular cylinders of equal radius, today called a bicylinder, is two thirds of the volume of its circumscribed cube [1]. (It seems have remained unnoticed that if the two cylinders intersect at an angle, then the volume of the solid is two thirds that of the circumscribed box. This follows at once by making an aﬃne transformation of ℝ3, which preserves ratios of volumes.) Later, modern mathematicians found the volume of the intersection of three such cylinders, a so-called tricylinder or “Steinmetz solid”, and its computation has caused countless desperation headaches in generations of calculus students. In 2017, Oliver Knill proposed, under the title “Archimedes’ Revenge,” that one prove that the volume of the intersection 푹 of the three (solid) hyperboloids

푥2 + 푦2 − 푧2 6 1; 푦2 + 푧2 − 푥2 6 1; 푧2 + 푥2 − 푦2 6 1 (1)

is equal to log 256. He proposed it as a challenge for the Harvard Maths 21A summer school in 2017. One of Knill’s students, Runze Li, oﬀered an integral [2] which computes the volume (see below); this motivated the solution we give here. arXiv:2108.05195v1 [math.HO] 6 Aug 2021

Figure 1: The solid contained by the three hyperboloids

1 2 The volume in the ﬁrst octant

By symmetry, the intersection 푹 of the three hyperboloids is the union of eight congruent solids, one in each octant. We will compute the volume of the component solid 푹1 in the ﬁrst octant. The intricate internal complexity of the solid 푹1 is shown by a brute-force triple integral for its volume 푉1, namely:

√ √ √ √ √︃ √ ( 2 2 2 2 2 1 2 2 ) ∫ 1/ 2 ∫ 푦 ∫ 푥 +푦 ∫ 1/ 2 ∫ 푥 +푦 ∫ 푦 + 2 ∫ 1−푥 +푦 푉1 = + √ + √ √ 1푑푦 푑푥 푑푧 0 0 0 푦 푦2−푥2 1/ 2 푥2−푦2 √ √ √ √ √ ∫ 3/2(∫ 1/ 2 ∫ 1+푥2−푦2 ∫ 푦 ∫ 1+푥2−푦2 ∫ 1 ∫ 1−푥2+푦2 ) 푑푦 푑푥 푑푧 + √ √︃ √ + √ √ + √ 1 2 1 2 2 2 2 2 2 1/ 2 푦 − 2 푦 −푥 1/ 2 푥 +푦 −1 푦 푥 +푦 −1 √ √ √ √ ∫ 1 (∫ 1/ 2 ∫ 1+푥2−푦2 ∫ 푦 ∫ 1+푥2−푦2 ∫ 1 ∫ 1−푥2+푦2 ) 푑푦 푑푥 푑푧. + √ √︃ √ + √ √ + √ 1 2 1 2 2 2 2 2 2 3/2 푦 − 2 푦 −푥 1/ 2 푥 +푦 −1 푦 푥 +푦 −1

This is by no means the end of the story. The ﬁrst two integrals in the second and third lines cannot be computed directly, but one must change the order of integration to calculate them. Thus a complete computation of the volume푉1 by this direct approach is quite daunting and tedious.

3 Symmetry

We will show how symmetry considerations reduce the computation of푉1 to a single integral! The solid 푹, whose total volume must be determined, is shown in Figure 1 – that appears on the webpage [2] of Oliver Knill, who kindly contributed the image. This Mathematica graphic shows its main features: a highly symmetrical solid whose boundary is a framework with several line segments supporting negatively curved surface patches taken from the three hyperboloids. (These line segments also form the edges of Kepler’s stella octangula.) The origin of ℝ3 is located at the (hidden) center of the solid. It evidently has the reﬂection symmetries 푥 ↔ −푥, 푦 ↔ −푦, 푧 ↔ −푧; and the cyclic symmetry (푥,푦, 푧) ↦→ (푦, 푧, 푥) of the 120◦ rotation around the diagonal line 푥 = 푦 = 푧. The key to understanding the solid in Figure 1 is that a one-sheeted hyperboloid is a ruled surface (indeed, a doubly ruled surface: it can be generated by either of two families of nonintersecting straight lines). Consider the intersection of any two of the hyperboloids of Eq. (1): 푦2 + 푧2 − 푥2 = 1 푧2 = 1 푧 = ±1 ⇐⇒ ⇐⇒ . 푧2 + 푥2 − 푦2 = 1 푥2 = 푦2 푥 ± 푦 = 0 The two lines 푧 = 1, 푥 ± 푦 = 0 are the horizontal lines at the upper boundary in Figure 1; the other two form the lower boundary at the bottom. The solid is constrained to lie inside the third hyperboloid 푥2 + 푦2 − 푧2 6 1, which cuts oﬀ these four lines at 푥2 + 푦2 6 2, yielding four line segments with eight endpoints (±1, ±1, ±1), all signs being allowed.

2 It should now be clear that the portion of the solid in the ﬁrst octant, 푹1, contains one vertex (1, 1, 1) of Kepler’s star, and three cross-points of the boundary segments, namely, the standard basis vectors of ℝ3. The solid 푹1 is composed of ﬁve pieces:

two back-to-back tetrahedra Π1 and Π2 with a common base, namely the equilateral triangle with vertices (1, 0, 0), (0, 1, 0), (0, 0, 1); and opposite vertices (0, 0, 0) and (1, 1, 1), respectively. Their union Π1 ∪ Π2 is a triangular dipyramid (Figure 2).

three congruent curved pieces 푆1, 푆2, 푆3, each of which is bounded by an outer face of the larger tetrahedron Π2 and one of the hyperboloids.

푧 (0, 0, 1) • (1, 1, 1) •

Π1 Π2

(0, 0, 0) • • 푦 (0, 1, 0)

• (1, 0, 0) 푥

Figure 2: The two adjacent tetrahedra: Π1 (hidden) and Π2 (in the foreground)

1 1 1 The volumes of the tetrahedra are: Vol(Π ) = and Vol(Π ) = , giving a total of . 1 6 2 3 2

4 Volume of one curved piece

We shall compute the volume of the solid piece 푆2 bounded by: the face of Π2 with vertices (1, 0, 0), (0, 0, 1), (1, 1, 1); a segment of the unit circle in the 푥푧-plane; and the hyperboloid

•

Figure 3: The curved piece 푆2 in the ﬁrst octant

3 2 2 2 푧 +푥 −푦 = 1. The face of Π2 lies in the plane 푥 −푦+푧 = 1. The projection of this tetrahedral face onto the 푥푦-plane is the triangle with vertices (0, 0), (1, 0), (1, 1). See Figure 3. The portion of the hyperboloid that forms the roof of 푆2 is shown in Figure 3 with both rulings by line segments. Its curved boundary (a quarter circle) can be parametrized by either (sin휃, 0, cos휃) or (cos 휙, 0, sin 휙) with the angles in the interval [0, 휋/2]. The rulings proceed from there to either of the straight sides, as follows:

푡 ↦→ (sin휃 + 푡 cos휃, 0, cos휃 − 푡 sin휃), 0 6 푡 6 sec휃 − tan휃; 푠 ↦→ (cos 휙 − 푠 sin 휙, 0, sin 휙 + 푠 cos 휙), 0 6 푠 6 sec 휙 − tan 휙.

The triple integral for the volume of the curved piece 푆2 is thus: √ ∫ 1 ∫ 1 ∫ 1+푦2−푥2 퐼 := 1푑푧 푑푥 푑푦. 0 푦 1−푥+푦

The inner integration is immediate: ∫ 1 ∫ 1 √︃ 퐼 = 1 + 푦2 − 푥2 − (1 − 푥 + 푦) 푑푥 푑푦. 0 푦 We write ∫ 1 ∫ 1 √︃ ∫ 1 ∫ 1 2 2 퐼1 := 1 + 푦 − 푥 푑푥 푑푦 and 퐼2 := (1 − 푥 + 푦) 푑푥 푑푦. 0 푦 0 푦 1 A routine calculation gives 퐼2 = . 3 √ Using the indeﬁnite integral of 푎2 − 푥2 with 푎 = √︁1 + 푦2 and then integrating by parts twice, we obtain

∫ 1 푥 √︃ 1 + 푦2 푥 푥=1 퐼 + 푦2 − 푥2 + 푑푦 1 = 1 arcsin √︁ 0 2 2 1 + 푦2 푥=푦 ∫ 11 + 푦2 1 1 + 푦2 푦 − 푑푦 = arcsin √︁ arcsin √︁ 0 2 1 + 푦2 2 1 + 푦2 휋 1 ∫ 1 푦 + 1푦3 휋 1 ∫ 1 푦 + 1푦3 = + 3 푑푦 − − 3 푑푦 2 2 6 2 0 1 + 푦 6 2 0 1 + 푦 ∫ 1 푦 + 1푦3 1 ∫ 1 2푦 log 2 1 = 3 푑푦 = + 푦 푑푦 = + . 2 2 0 1 + 푦 3 0 1 + 푦 3 6

Thus the volume of the curved piece 푆2 is log 2 1 1 log 2 1 퐼 = 퐼 − 퐼 = + − = − . 1 2 3 6 3 3 6

4 5 Total Volume

Therefore the volume of the solid 푹1 in the ﬁrst quadrant is 1 log 2 1 1 푉 = 3퐼 + = 3 − + = log 2; 1 2 3 6 2 and the total volume of the solid 푹 in Archimedes’ revenge is:

Vol(푹) = 8푉1 = 8 log 2 = log 256. qed!

6 Comment

The idea to exploit the symmetry of a component in one octant was already suggested 1 by Li [2]. His solution claims that the solid 푹1 is composed of one tetrahedron of volume 2 (a minor oversight) and the three congruent curved pieces. He computes the volume of a curved piece using the following integral whose origin is somewhat mysterious, namely:

1 ∫ 1 휋 퐼 := (푧2 + 1) − 2 arctan 푧 + 푧2 − 1 푑푧 2 0 2 and indeed our attempt to decipher it led us to the solution oﬀered here. (We would welcome any reader’s commentary and explanation for this integrand.)

Acknowledgements

We are grateful to Oliver Knill for sharing the image in Figure 1 that appears on his webpage [2] (which also includes Runze Li’s suggested solution), and for helpful comments. We also acknowledge support from the Vicerrectoría de Investigación of the Universidad de Costa Rica.

References

[1] Archimedes, The Works of Archimedes: The Method, T. L. Heath, ed., Dover Publica- tions, New York, 1960.

[2] From the webpages of Oliver Knill, 2017: http://people.math.harvard.edu/ ~knill/teaching/summer2017/exhibits/revenge/