1. Motivation and statement of results
A classical theorem states that 
$\mathbb{R}^3$ can be decomposed as a union of disjoint unit circles. This statement may be folklore, but it seems to have been popularized by Conway and Croft, who gave a non-constructive proof using transfinite induction [Reference Conway and Croft10]. Since then, a number of fascinating questions, regarding the possibility of decomposing a fixed space as the union of disjoint copies of another space, have been studied from set theoretic, topological and geometric perspectives. A collection of decomposition questions appear in an expository article by Gardner [Reference Gardner11], whereas a collection of geometric constructions using transfinite induction can be found in [Reference Komjáth and Pach19]. Similar questions appear every so often on MathOverflow; see [20, 22–24]. It seems to be unknown whether there exists a Borel decomposition of 
$\mathbb{R}^3$ into unit circles; see [21].
With different topological or geometric constraints, versions of this question have appeared in a wide range of articles [Reference Bankston and McGovern2–Reference Cobb9, Reference Jonsson and Wästlund18, Reference Moussong and Simányi26, Reference Wilker31], and several explicit constructions have appeared. For example, Szulkin constructed an explicit covering of 
$\mathbb{R}^3$ using geometric (round) circles of varying radii, which fails to be a foliation only on the union of a countable discrete set of concentric circles [Reference Szulkin29], whereas foliations of 
$\mathbb{R}^3$ by topological circles were thoroughly studied by Vogt [Reference Vogt30]. Bankston and Fox [Reference Bankston and Fox1], generalizing a construction of Kakutani, gave a decomposition of 
$\mathbb{R}^{n+2}$ by tamely embedded, unknotted, pairwise unlinked (topological) n-spheres, but it seems unknown whether there exists any decomposition of 
$\mathbb{R}^{n+2}$ by unit n-spheres for n > 1.
Here, we address the existence question for continuous unit sphere decompositions of open subsets E in Euclidean space 
$\mathbb{R}^d$. A decomposition of E by unit n-spheres is called continuous if the map taking 
$p \in E$ to its containing sphere is continuous; that is, if the sphere centres and their normal spaces vary continuously with p. In this case, we say that E is fibered by unit n-spheres. Global fibrations of 
$\mathbb{R}^d$ by unit spheres cannot exist for topological reasons, as can be seen from the long exact sequence associated to the fibration. However, the existence question for local fibrations, has, to the authors’ surprise, turned out to be non-trivial.
A simple example of a unit circle fibration can be described as follows. Consider an open torus 
$\textstyle T\subset\mathbb{R}^3$ with major and minor radii both equal to 1 (the boundary of T has a singularity in the centre). The interior of this torus may be foliated by a collection of tori Tr, each with major radius 1, but with minor radii r ranging from 0 to 1. Now each individual torus admits two foliations (a left-handed one and a right-handed one) by Villarceau circles. The radius of a Villarceau circle is equal to the major radius of the torus, and hence all such circles have radius 1. Choosing the right-handed foliation on each torus yields a collection of pairwise linked unit circles which fill the interior of T. See Figure 1.
Fibration of a toroidal region in 
$\mathbb{R}^3$ by linked unit circles.
Remark 1. We note that the stereographic projection of the Hopf fibration of S 3 by unit circles produces an image which looks very similar to that of Figure 1. While the circles there are also Villarceau circles on tori, they are not unit circles.
In 
$\S$3, we construct explicit unit n-sphere fibrations in 
$\mathbb{R}^{2n+1}$, for 
$n \in \left\{1, 3, 7 \right\}$. When n = 1 this construction produces the fibered torus described above. The dimensions of our construction are sharp, as shown by our main result:
Theorem 1. Suppose that there exists a unit n-sphere fibration of an open subset 
$E \subset \mathbb{R}^d$. Then 
$d \geq 2n+1$, and if 
$d=2n+1$, then 
$n \in \left\{0, 1, 3, 7 \right\}$.
If there exists a unit n-sphere fibration of an open region 
$E \subset \mathbb{R}^d$, then there exists a unit n-sphere fibration of 
$E \times \mathbb{R} \subset \mathbb{R}^{d+1}$, by ‘stacking’ unit n-sphere fibrations of E. Thus, it would be interesting to determine:
Question 1. Given 
$n \in \mathbb{N}$, do there exist d and an open subset 
$E \subset \mathbb{R}^d$ such that E is fibered by unit n-spheres? If so, what is the minimum such dimension d?
The answer to this question is known for two other types of geometric fibrations, and it would not be surprising if the answer to Question 1 matches one of these described below:
A great sphere fibration is a sphere bundle with total space Sd and fibres which are great n-spheres. The standard example is the Hopf fibration of S 3 by unit circles, which arises by choosing an orthogonal complex structure on 
$\mathbb{R}^4$ and intersecting S 3 with all of the complex lines. Similar constructions with quaternions or octonions yield Hopf fibrations with fibres S 3 or S 7. Algebraic topology imposes strong restrictions on the possible dimensions of great sphere fibrations; in particular the fibres must have dimension 0, 1, 3 or 7. With this in mind, it would not be surprising if only these n could serve as the fibre dimension of a unit sphere fibration. See [Reference Gluck and Warner12–Reference Gluck, Warner and Ziller14, Reference McKay25] for more information about great sphere fibrations.
A skew fibration is a vector bundle with total space 
$\mathbb{R}^d$ and fibres which are pairwise skew affine n-planes. In [Reference Ovsienko and Tabachnikov27], Ovsienko and Tabachnikov showed that a skew fibration of 
$\mathbb{R}^d$ by n-planes exists if and only if 
$n \leq \rho(d-n) - 1$, where ρ is the Hurwitz–Radon function, defined as follows: Decompose 
$q \in \mathbb{N}$ as the product of an odd number and 
$2^{4a+b}$, where 
$0 \leq b \leq 3$, then 
$\rho(q) = 2^b + 8a$. It follows from unboundedness of ρ that for every n there exists d such that 
$\mathbb{R}^d$ may be fibered by skew affine n-planes. With this in mind, it would not be surprising if every n could serve as the fibre dimension of a unit sphere fibration. See [Reference Harrison15–Reference Harrison17, Reference Ovsienko and Tabachnikov28] for more information about skew fibrations.
Remark 2. In the proof of Theorem 1, we will demonstrate a correspondence between linked unit sphere fibrations and skew fibrations. In particular, we will use the fact that skew fibrations of 
$\mathbb{R}^{2n+1}$ by n-planes exist if and only if 
$n + 1 = \rho(n+1)$, which occurs if and only if 
$n \in \left\{0, 1, 3, 7 \right\}$.
Although we have found some relationships among unit sphere fibrations, great sphere fibrations and skew fibrations, we have found that the techniques used to study unit sphere fibrations are somewhat different from those used to study the other types of fibrations. For example, the collection of n-spheres in a great sphere fibration of Sd corresponds to a submanifold of 
$\operatorname{Gr}_{n+1}(d+1)$, and an important component of studying great sphere fibrations is understanding the topology of the Grassmann manifold. Similarly, the study of skew fibrations requires understanding the topology of the affine Grassmann and spaces of non-singular bilinear maps. By contrast, geometry plays a large role in the study of unit sphere fibrations, and it seems unlikely that their study can be completely reduced to questions of topology or linear algebra.
2. Proof of the main result
The proof of the main result uses the notion of linkedness for unit n-spheres. Recall that two disjoint topological n-spheres S 1 and S 2 in 
$\mathbb{R}^d$ are linked if S 1 is non-trivial as an element of 
$\pi_n(\mathbb{R}^d - S_2)$. Linking of topological n-spheres can potentially occur in dimensions 
$n+2 \leq d \leq 2n+1$, but linking for unit n-spheres is much more restrictive.
Lemma 1. Let S 1 and S 2 be two disjoint unit n-spheres in 
$\mathbb{R}^d$, 
$n+2 \leq d \leq 2n+1$, and let Pi be the affine 
$(n+1)$-plane containing Si. The spheres S 1 and S 2 are linked if and only if S 2 intersects P 1 exactly once in each connected component of 
$P - S_1$. In particular, two unit n-spheres can only be linked in dimension 
$d = 2n+1$.
Proof. If S 2 does not intersect P 1 inside S 1, then S 2 lives in the complement of the unit disk D 1 in 
$\mathbb{R}^d$ and hence can be contracted to a point inside 
$\mathbb{R}^d - S_1$. Similarly, S 2 must intersect P 1 outside S 1.
The intersection 
$S_2 \cap P_1$ is a sphere or a single point. Therefore, if S 2 intersects P 1 in more than two points, at least one in each component, then it intersects P 1 in at least a circle, which contains points in each connected component of 
$P_1 - S_1$ and hence intersects S 1.
Similarly, if 
$d \lt 2n$, then P 1 and P 2 must intersect in at least a 2-dimensional space, and so if S 2 intersects P 1 in one point in each component, it must intersect P 1 in at least a circle. Therefore, if S 2 intersects P 1 exactly once in each connected component, then 
$d= 2n+1$, and it is not difficult to see, by rotating P 2, that there is a homotopy bringing S 2 to the position of a standard generator for 
$\pi_n(\mathbb{R}^{2n+1} - S_1)$.
Now we are prepared to prove the main result. The proof splits naturally into two separate arguments, which we record as the following two lemmas.
Lemma 2. Suppose that there exists a unit n-sphere fibration of an open subset 
$E \subset \mathbb{R}^d$. Then 
$d \geq 2n+1$, and if 
$d=2n+1$, then the spheres are pairwise linked.
The final statement of Lemma 2, in the case n = 1, formalizes the intuitive idea that a 1-parameter family of unit circles cannot ‘pass through’ another unit circle in 
$\mathbb{R}^3$.
Lemma 3. Suppose that there exists a fibration of an open connected subset 
$E \subset \mathbb{R}^{2n+1}$ by unit, linked, n-spheres. Then 
$n \in \left\{1, 3, 7 \right\}$.
We prove the first lemma with a Borsuk–Ulam construction. For the second lemma we develop and apply a correspondence between unit sphere fibrations and skew fibrations.
Proof of Lemma 2
Let p be a point in the interior of E. We use the following notation: C is the unit n-sphere fibre through p, B is the open 
$(n+1)$-ball whose boundary is C, P is the 
$(n+1)$-plane containing C, and 
$D \subset E$ is a small 
$(d-n)$-dimensional disk centred at p, orthogonal to the tangent space to C at p, and small enough such that every fibre intersecting D does so transversely and exactly once; see Figure 2.
Fibration by unit circles: local depiction.
By compactness of C and openness of E, there exists ɛ > 0 such that the n-sphere 
$S \subset P$ concentric with C and of radius 
$1 + \varepsilon$ lies inside E. Suppose momentarily that for each 
$x \in S$, the fibre passing through x intersects the boundary 
$\partial D$ of the small transverse disk at p. This allows us to define 
$f \colon S \to \partial D \simeq S^{d-n-1}$ as the map which sends 
$x \in S$ to the intersection of the fibre through x with 
$\partial D$.
If 
$d \leq 2n$, i.e. if 
$d-n -1 \lt n$, then by the Borsuk–Ulam theorem, there exists x such that 
$f(x) = f(-x)$, and so the fibres through x and −x intersect. Since x and −x are distance 
$2 + 2\varepsilon$ apart, their containing fibres are different and hence do not intersect, a contradiction.
Now let 
$d = 2n+1$. Since linkedness is preserved by homotopy, there is a dichotomy: all fibres passing through 
$\partial D$ are linked with C, or all fibres are unlinked with C. We assume all fibres are unlinked and derive a contradiction. Since C is unlinked with the fibre through 
$x \in S$, and since 
$x \in S$ lies in the exterior of 
$C \subset P$, then by Lemma 1, the fibre through x does not intersect B. In particular, the image of f misses the point 
$\partial D \cap B$. Therefore, f may be considered as a map 
$S^n \to \mathbb{R}^n$, and the Borsuk–Ulam theorem applies in the same way.
Now we address the assumption that for each 
$x \in S$, the fibre passing through x intersects the boundary 
$\partial D$ of the small transverse disk at p. This need not occur, and so we modify the definition of S as follows. Since D is chosen to correspond to a local trivialization, the fibres passing through D generate a foliated neighbourhood 
$U \subset E$ of the fibre C, and the fibres passing through 
$\partial D$ form a fibered topological torus 
$T = S^{d-n-1} \times S^n$ which serves as the boundary of U. For 
$q \in C$, following fibres from 
$\partial D$ yields a homeomorphism between 
$\partial D$ and the boundary of a topological 
$(d-n)$-dimensional disk Dq which is transverse to C at q (explicitly, Dq is the intersection of U with the 
$(d-n)$-plane orthogonal to the tangent space TqC). Now let x be the point on the line segment 
$D_q \cap P$ which is farthest from the centre c of C. By compactness of 
$\partial D_q$, this line segment intersects 
$\partial D_q$ in finitely many points, and the line segment intersects 
$\partial D_q$ at least once on each side of q, since C lies in the interior connected component of 
$\mathbb{R}^d - T$. Therefore x is well-defined, and moreover it varies continuously with q. Said differently, if δq represents the maximum distance from 
$D_q \cap P$ to c, then 
$x = c + \delta_q q$. The map 
$q \mapsto c + \delta_q q = x$ is a homeomorphism from C onto its image Sʹ, and the proof goes through with this topological n-sphere Sʹ in place of the geometric sphere S.
Remark 3. It is tempting to apply Borsuk–Ulam to some natural map on the n-sphere 
$\partial D$, instead of on the topological n-sphere Sʹ obtained by fattening the fibre C, but we were unable to define a map which simultaneously captures both the unit assumption and the unlinkedness assumption. It would be interesting to know if the same obstruction could be found by such a map.
Proof of Lemma 3
We show for n > 0 that the fibration of E induces a fibration of some open set 
$U \subset \mathbb{R}^{2n+1}$ by skew affine copies of 
$\mathbb{R}^n$, which is known to only exist when 
$n \in \left\{1,3,7\right\}$ (see Remark 2).
For each spherical fibre F, consider the containing 
$(n+1)$-plane P, and then consider the orthogonal n-plane Q passing through the centre of F. Let Q 1 and Q 2 be two such n-planes obtained from fibres F 1 and F 2 with corresponding planes P 1 and P 2. We claim that if Q 1 and Q 2 are not skew, then the fibres F 1 and F 2 are unlinked.
Suppose first that Q 1 and Q 2 share a common line 
$\ell$, i.e. their linear representatives span fewer than 2n dimensions. Then 
$\ell$ is orthogonal to both P 1 and P 2, so P 1 and P 2 do not span a 
$(2n+1)$-dimensional space. If P 1 and P 2 do not intersect, then F 1 and F 2 are not linked. If P 1 and P 2 intersect, then F 1 and F 2 both lie in some 2n-dimensional space, hence are not linked by Lemma 1.
Suppose now that the linear representatives of Q 1 and Q 2 span a 2n-dimensional space, but that Q 1 and Q 2 intersect. Then the vector u connecting the centres of F 1 and F 2 lies in the linear span of Q 1 and Q 2. Moreover, P 1 and P 2 span a 
$(2n+1)$-dimensional space and intersect in some line 
$\ell$, and 
$\ell$ is orthogonal to u. To verify that F 1 and F 2 are unlinked, it suffices to understand the intersections of 
$F_1 \cap P_2$ and 
$F_2 \cap P_1$. If either intersection is empty, the spheres are unlinked. Each non-empty intersection occurs within 
$P_1 \cap P_2 = \ell$ and must consist of one or two points.
Let xi be the point on 
$\ell$ nearest to the centre of Fi. Note that xi can be obtained by taking the intersection of 
$\ell$ with the copy of the hyperplane orthogonal to 
$\ell$ passing through the centre of Fi. Since 
$\ell$ is orthogonal to u, the centres of F 1 and F 2 each lie on the same such hyperplane, and therefore 
$x_1 = x_2$.
Now the intersection of each Fi with 
$\ell$ is symmetric about x 1. Therefore, F 1 cannot intersect P 2 both outside and inside F 2, so the spheres are not linked by Lemma 1.
It remains to observe that continuity of the unit sphere fibration leads to continuity of the induced skew covering, and hence yields a bona fide fibration. Indeed, the skew plane corresponding to a point 
$p \in E$ is defined by taking the affine plane through the centre c(p) and with linear representative 
$\nu(p) \in \operatorname{Gr}_{n}(2n+1)$, each of which vary continuously with p. Hence the collection of affine planes varies continuously with p.
3. Example
In this section, we construct an example of a unit n-sphere fibration of 
$\mathbb{R}^{2n+1}$ for 
$n \in \left\{1, 3, 7 \right\}$. The construction was inspired by the relationship between unit sphere fibrations and skew fibrations exhibited in the proof of Lemma 3.
Example 1. We construct an example of a unit sphere fibration for n = 3. The same construction works for n = 1 and n = 7 with complex numbers or octonions in place of quaternions. Write 
$\mathbb{R}^7 = \mathbb{R}^4 \times \mathbb{R}^3$ and for each 
$(y,0) \in \mathbb{R}^4 \times \mathbb{R}^3$, let Py be the 3-plane spanned by 
$(iy,1,0,0)$, 
$(jy,0,1,0)$, 
$(ky, 0, 0, 1)$, where 
$\left\{i, j, k \right\}$ represents a basis of the imaginary unit quaternions. Let Qy be the 4-plane which passes through y and is orthogonal to Py, and let Sy be the unit 3-sphere in Qy centred at y. Now the centre map 
$y \mapsto (y,0) \in \mathbb{R}^7$ and the normal map 
$y \mapsto Q_y$ are continuous, and therefore the assignment 
$y \mapsto S_y$ defines a continuous unit sphere fibration provided that no two spheres intersect.
We will show that for distinct 
$y,z$ with 
$|y|,|z| \lt 1$, Sy and Sz do not intersect.
First note that for distinct y and z, the linear span of Py and Pz is 6-dimensional. Therefore, the affine 4-planes Qy and Qz intersect in precisely a line 
$\ell$, and it is enough to check that 
$S_y \cap \ell$ and 
$S_z \cap \ell$ do not intersect. Note that 
$\ell$ contains the origin, and it is easy to check, by taking inner products with the basis vectors of Py and Pz, that 
$\ell$ has direction:
\begin{equation*}
v = (y-z,- \langle iy, z \rangle, - \langle jy,z \rangle, - \langle ky, z \rangle).
\end{equation*}For future convenience, we compute,
\begin{align}
|v|^2 = |y-z|^2 + \langle iy, z \rangle^2 + \langle jy,z \rangle^2 + \langle ky, z \rangle^2 = |y-z|^2 + |z|^2|y|^2 - \langle y,z \rangle^2,
\end{align} where we have used the Parseval identity together with the fact that 
$\left\{y, iy, jy, ky \right\}$ forms an orthogonal basis of 
$\mathbb{R}^4$. The nearest point on 
$\ell$ to the point y is:
\begin{equation*}
c_y = \frac{\langle y, v \rangle }{|v|^2} v = \frac{\langle y, y-z \rangle}{|v|^2} v,
\end{equation*}and the squared distance from y to cy is:
\begin{equation*}
d_y^2 = \langle y - c_y, y - c_y \rangle = |y|^2 - 2\langle y, c_y \rangle + |c_y|^2 = |y|^2 - \frac{\langle y, y-z \rangle^2}{|v|^2}.
\end{equation*} Now, the intersection points of Sy with 
$\ell$ are the points 
$c_y \pm r_y \frac{v}{|v|}$, where 
$r_y = \sqrt{1 - d_y^2}$. Let d represent the distance from cy to cz. Then
\begin{equation*}
d = |c_y - c_z| = \frac{|y-z|^2}{|v|}.
\end{equation*} To check that Sy and Sz do not intersect, we must show that 
$c_y \pm r_y\frac{v}{|v|} \neq c_z \pm r_z\frac{v}{|v|}$. For this, it suffices to show that 
$d \lt r_y + r_z \lt \min\left\{2r_y + d, 2r_z + d \right\}$; see Figure 3.
Two non-intersecting spheres centred at y and z; the red points represent 
$c_y \pm r_y\frac{v}{|v|}$, and the blue points represent 
$c_z \pm r_z\frac{v}{|v|}$.
Lemma 4. In the notation above, we have
(a)
$d \lt r_y + r_z$,(b)
$r_z \lt r_y + d$,(c)
$r_y \lt r_z + d$.
Proof. The inequality in part (a) can be written explicitly:
\begin{equation*}
\frac{|y-z|^2}{|v|} \lt \sqrt{1 - |y|^2 + \frac{\langle y, y-z \rangle^2}{|v|^2}} + \sqrt{1 - |z|^2 + \frac{\langle z, y-z \rangle^2}{|v|^2}},
\end{equation*}or
\begin{align}
|y-z|^2 \lt \sqrt{|v|^2 - |v|^2|y^2| + \langle y, y - z \rangle^2} + \sqrt{|v|^2 - |v|^2|z^2| + \langle z, y - z \rangle^2}.
\end{align} Since 
$|y| \lt 1$, we have
\begin{equation*}
0 \lt |v|^2-|v|^2|y|^2,
\end{equation*}therefore,
\begin{equation*}
\langle y, y-z \rangle^2 \lt |v|^2 - |v|^2|y^2| + \langle y, y - z \rangle^2,
\end{equation*}and so,
\begin{align}
\langle y, y-z \rangle \lt \sqrt{|v|^2 - |v|^2|y^2| + \langle y, y - z \rangle^2}.
\end{align}Similarly,
\begin{align}
-\langle z, y-z \rangle \lt \sqrt{|v|^2 - |v|^2|z^2| + \langle z, y - z \rangle^2}.
\end{align}Now, summing (3.3) and (3.4) yields (3.2).
The inequality in part (b) can be written explicitly:
\begin{equation*}
\sqrt{1 - |z|^2 + \frac{\langle z, y-z \rangle^2}{|v|^2}} \lt \sqrt{1 - |y|^2 + \frac{\langle y, y-z \rangle^2}{|v|^2}} + \frac{|y-z|^2}{|v|}.
\end{equation*} Multiplying by 
$|v|$ and using (3.1) yields:
\begin{equation*}
\sqrt{|y-z|^2 - |z|^2(|y|^2|z|^2 - \langle y, z \rangle^2)} \lt \sqrt{|y-z|^2 - |y|^2(|y|^2|z|^2 - \langle y, z \rangle^2)} + |y-z|^2.
\end{equation*}Squaring both sides and rearranging, we obtain the inequality:
\begin{align}
(|y|^2-|z|^2)(|y|^2|z|^2 - \langle y, z \rangle^2) \lt 2\sqrt{|y-z|^2 - |y|^2(|y|^2|z|^2 - \langle y, z \rangle^2)}|y-z|^2 + |y-z|^4.
\end{align} For convenience let 
$a = |y|$, 
$b=|z|$ and 
$c = \cos \theta$, where θ is the angle between y and z. We write
\begin{align*}
0 & \geq -(a-bc)^2-(ac-b)^2 \\ & = -a^2-b^2-a^2c^2-b^2c^2+4abc \\
& = (a^2+b^2)(1-c^2) -2(a^2-2abc+b^2)\\ & \geq (a^3b^2+a^2b^3)(1-c^2)-2(a^2-2abc+b^2), \hspace{.2in} \mbox{since } a,b \lt 1.
\end{align*}Now, we may assume a > b, since otherwise (3.5) is immediate. Then, multiplying the previous inequality by a − b yields:
\begin{align*}
0 & \geq (a-b)(a^3b^2+a^2b^3)(1-c^2)-2(a-b)(a^2-2abc+b^2) \\
& = (a^4b^2-a^2b^4)(1-c^2) - 2(a-b)(a^2-2abc+b^2).
\end{align*}Therefore,
\begin{align*}
(a^4b^2-a^2b^4)(1-c^2)& \leq 2(a-b)(a^2-2abc+b^2) \\
& = 2(a-b)|y-z|^2 \\
& = 2\sqrt{a^2-2ab+b^2}|y-z|^2 \\
& \leq 2\sqrt{a^2-ab(1+2c-c^2)+b^2}|y-z|^2, \hspace{.5in} \mbox{since } -1 \leq c \leq 1,\\
& = 2\sqrt{a^2 - 2abc + b^2 -ab(1-c^2)}|y-z|^2\\
& \leq 2\sqrt{a^2 - 2abc + b^2 -a^4b^2(1-c^2)}|y-z|^2\\
& = 2\sqrt{|y-z|^2 -a^4b^2(1-c^2)}|y-z|^2\\
& \lt 2\sqrt{|y-z|^2 -a^4b^2(1-c^2)}|y-z|^2 + |y-z|^4,
\end{align*}and this is (3.5).
Part (c) is similar.
Remark 4. It is possible that the relationship between linked unit sphere fibrations and skew fibrations can be formalized using certain transformations of the ambient space, though such transformations would necessarily need to carry both geometric and topological data, or perhaps be dimension-specific. For example, applying such a transformation to a portion of a skew line fibration of 
$\mathbb{R}^5$ cannot yield a linked unit circle fibration since circles cannot be linked in 
$\mathbb{R}^5$, and applying such a transformation to an unlinked unit circle fibration in 
$\mathbb{R}^4$ cannot yield a skew line fibration of 
$\mathbb{R}^4$, since no such fibration exists. An anonymous referee suggested that Lie Sphere Geometry may contain the tools to formalize these relationships.
4. History
Questions regarding the possibility of fibering an open set in 
$\mathbb{R}^3$ by unit circles were raised by the first named author (DA) as a graduate student at UC Berkeley circa 1970. No one questioned had answers, but in the early 1990s DA discovered the fibering of the toroidal region depicted in Figure 1 (he called this region a ‘bialy’ in reference to the bagel-like cuisine). He gave a series of talks on his findings, which were published in a 1993 technical report ‘Hoops in 
$\mathbb{R}^3$’, though he later retracted this report after discovering an error.
In 2019, the other three authors (OA) were introduced to this type of problem at Carnegie Mellon, when colleague Anton Bernshteyn posed to them the existence question for Borel coverings of 
$\mathbb{R}^3$ by unit circles. OA became interested in the continuous case, and a literature search resulted in only a single finding: a 1995 blog post by Evelyn Sander, which she wrote to summarize a 1994 talk by DA at UIUC. Her post mentioned the announced results of DA, but OA were unable to extract any formal proofs, though through this post they did learn the importance of linkedness and the potential relationship to division algebras. They were unable to access the technical report nor find any contact information for DA.
OA later discovered that a mutual colleague had contact information for DA. OA shared a draft of this article and learned that DA also knew (though had not published) proofs of some of these results. This led to the present collaboration, which we hope provides some closure to a 50 year old question – though of course many interesting questions remain!
Acknowledgements
F. F. was supported by NSF grant DMS 1855591, NSF CAREER grant DMS 2042428 and a Sloan Research Fellowship. M. H. was supported by NSF grant DMS 1926686. W. P. was supported by NSF grant DMS 2054503.
