# Three-dimensional Apollonian spirals

Given mutually (externally) tangent spheres $S_1$, $S_2$, $S_3$, $S_4$, let $S_n$ be the unique sphere externally tangent to $S_{n-1}$, $S_{n-2}$, $S_{n-3}$, and $S_{n-4}$ for $n \geq 5$. Let $P_{\infty}$ be the point toward which the $S_n$'s tend, let $P_n$ be the center of $S_n$, and let $Q_n$ be the point of tangency between $S_{n}$ and $S_{n+1}$. What can be said about the asymptotic distribution of the unit vectors pointing from $P_{\infty}$ to $P_n$? What about the unit vectors pointing from $P_{\infty}$ to $Q_n$?

This was part of what I had in mind when I posted Spirals in Apollonian circle-packings. My hope is that this process (or something like it) gives a way to generate uniformly-distributed points on the sphere; my fear is that (as in the case of a rotating ellipsoid or rectangular parallelepiped) there's a kind of instability at play, so that asymptotically one gets a uniform distribution on some great circle rather than the whole sphere.