The Calculus of Explanations

Putting the fun in functionals


The longest arc: Part II

Last time we saw that the problem of finding the circle with the longest arc-length wasn’t possible analytically, and we had to resort to numerical or graphical methods.

Now we extend the problem to find the sphere centered at the origin with the largest surface area inside a unit sphere shifted vertically by one unit.

Which sphere has the biggest surface area? See a demonstration of the problem in Desmos.

Here we define a sphere (S_1) centered at \left(0,0,1\right):

\displaystyle S_1\::\:x^2 + y^2 + \left(z - 1\right)^2 = 1^2

And the ‘expanding’ sphere (S_2) centered at the origin with radius \rho:

\displaystyle S_2\::\:x^2 + y^2 + z^2 = \rho^2

Where \rho \in \mathbb{R}^+. Here we seek to solve for \rho that optimises the surface area of the intersection of S_1 and S_2. We first solve for the intersection of S_1 and S_2:

\begin{aligned}S_1\::\: \:  x^2 + y^2 &= 1 - \left(z - 1\right)^2, \\ S_2\::\: \: x^2 + y^2 &= \rho^2 - z^2.\end{aligned}

Setting the right hand sides equal gives us that z=\rho^2/2 for the points at which the two spheres intersect.

To identify the region this corresponds to in x and y we substitute z = \rho^2/2 into either S_1 or S_2. Here, we choose S_2 for simplicity:

\begin{aligned} x^2 + y^2 + z^2 &= \rho^2 \\ x^2 + y^2 + \left( \frac{\rho^2}{2} \right)^2 &= \rho^2 \\ x^2 + y^2 &= \rho^2 - \frac{\rho^4}{4} \\ x^2 + y^2 &= \frac{1}{4}  \rho^2 \left(4 -  \rho^2\right) \end{aligned}

Here we observe that S_1 and S_2 intersect at z = \rho^2/2 over a circle centered at the origin with radius \frac{1}{2} \rho \sqrt{4 - \rho^2}. Call this region R.

We first observe that the radius of the expanding sphere is \rho < 2. From there, we can define the surface area A(\rho) of S_2 over R as the double integral:

A\left(\rho\right) = \displaystyle \iint_R \sqrt{1 + \left(\frac{\partial z}{\partial x}\right)^2 + \left(\frac{\partial z}{\partial y}\right)^2}\:\mathrm{d}R

Recall for S_2 that

\begin{aligned}x^2 + y^2 + z^2 &= \rho^2, \\ \frac{\partial}{\partial x}\left( x^2 + y^2 + z^2\right) &= \frac{\partial}{\partial x}\left( \rho^2\right), \\ 2x + 2z \displaystyle\frac{\partial z}{\partial x} &= 0, \\ \frac{\partial z}{\partial x} &= -\frac{x}{z}.\end{aligned}

Similarly, we have \displaystyle \frac{\partial z}{\partial y} = -\frac{y}{z}. Thus, we have:

\begin{aligned}\displaystyle A\left(\rho\right) &= \iint_R \sqrt{1 + \left(-\frac{x}{z}\right)^2 + \left(-\frac{y}{z}\right)^2}\:\mathrm{d}R \\ &= \iint_R \frac{\sqrt{x^2 + y^2 + z^2}}{z} \:\mathrm{d}R \\&= \iint_R \frac{\rho}{\sqrt{\rho^2 - \left(x^2 + y^2\right)}} \:\mathrm{d}R \end{aligned}

This integral is crying out for a substitution, into the polar coordinates given by x / \rho=r\cos \theta and y / \rho = r\sin \theta. Under this transformation the integral becomes

\begin{aligned}\displaystyle A\left(\rho\right) &= \displaystyle \rho^2\int_{-\pi}^{\pi}\:\mathrm{d}\theta \int_{0}^{\frac{\sqrt{4 - \rho^2}}{2}} \frac{r}{\sqrt{1 - r^2}}\:\mathrm{d}r \\ &=\rho^2\left[ \theta\right]_{-\pi}^{\pi} \left[-\sqrt{1 - r^2} \right]_0^{\frac{\sqrt{4 - \rho^2}}{2}} \\ &= \rho^2 \cdot 2\pi \left(1 - \sqrt{1 - \left(\frac{\sqrt{4 - \rho^2}}{2}\right)^2}\right)\\ &= 2\pi\rho^2 \left(1 - \frac{\rho}{2}\right) \\ A\left(\rho\right)  & = 2\pi\left(\rho^2 - \frac{\rho^3}{2}\right) \end{aligned}

A remarkably simple formula for the surface area of the expanding sphere!

To find the one with the maximum surface area, we take the derivative and set it equal to zero,

2\rho - \displaystyle \frac{3\rho^2}{2} = 0

This equation has two solutions; \rho = 0 (clearly a minimum), and \rho = \displaystyle \frac{4}{3}, which is the radius of the sphere with maximum surface area. If you like, show the second derivative is negative to prove this is a local maximum!

If we put this value back into our surface area formula we get that the maximum area is

A\left(\displaystyle \frac{4}{3}\right) = \displaystyle 2\pi\left(\left(\frac{4}{3}\right)^2 - \displaystyle \frac{1}{2}\left(\frac{4}{3}\right)^3\right) = \displaystyle \frac{32}{27}\pi.

As we mentioned in the previous post, this seems to be a problem that (in some sense) becomes easier in one higher dimension.

However you feel about double integrals aside, it is interesting that in the case of the circles, the equation to solve for the optimal radius was analytically intractable, whereas for spheres, it is a simple quadratic.

This allows us to achieve a closed form solution for the maximum area in 3d, something that wasn’t possible in 2d.

Actually, we think that n=3 dimensions is the only time this happens, and in all other n it becomes necessary to solve A'(\rho) = 0 numerically.

Interestingly, this seems to be true for different reasons depending on whether the dimension n is odd or even.

But that’s a story for next time.

Stay curious.



One response to “The longest arc: Part II”

  1. […] in for the next post where we’ll be looking at the same problem but with spheres instead of circles!In the […]

    Like

Leave a comment

Hi there! I’m a creator with a background in maths education and a passion for making complicated ideas simple, with gorgeous animations and visual explanations.

Whether you’re just entering high school or you’re a post-grad student, or just have an interest in maths in general, I aim to make content that can appeal to everyone, in the form of videos and blog posts.

You can find all my previous posts here.

Check out my YouTube channel for more maths education content.