Where my questions stem from:
When we study the rotations in a plane or of some specific higher dimensions, there exists a neat approach to represent all the rotations as a spheres $\mathbb S^i$, for appropriate $i$.
For instance, rotations of the complex plane is $\mathbb S^1$, with a group structure over that, where $\mathbb S^3$ is $\{a+bi|a^2+b^2=1\}$. Similarly, all rotations of the 3-dimensional pure imaginary subspace $\mathbb Ri+ \mathbb Rj+ \mathbb Rk$ in the 4-dimensional space $\mathbb H$ is $\mathbb S^3$, using the correspondence between rotations of the subspace and quaternions of absolute value one. Namely, $\mathbb S^3$ is $\{a+bi+cj+dk|a^2+b^2+c^2+d^2=1\}$, equipped with a group structure to encodes the rotations.
In order to construct the group structures which we put on $\mathbb S^1$ and $\mathbb S^3$ as the groups of rotations, we use the Brahmagupta–Fibonacci identity and Euler's four-square identity, respectively, which give us the identities we need to insure the multiplications of two of the elements from each of those spheres would be another element of absolute value one, and consequently they form a group.
Having this pattern in mind, here are my main questions:
I have seen an argument which tries to give a sense of the main obstacle we could not tackle to use the same approach for any higher dimension, which seems to be the absence of similar identities for an arbitrary $n$ sum of squares. For instance, there is a product of two numbers, each of which is a sum of seven squares, while itself is not a sum of seven squares. So, the multiplication fails since we don't have the same fact about the size of two elements of $\mathbb S^6$.
However, there are analogous identities, in particular, Degen's eight-square idnetity and Pfister's sixteen-square identity, which could potentially provide us with the chance to imitate the same method, at least, for some specific higher dimensions. However, presumably we deserve less chance for case $n=16$, since, unlike the cases $n=1,2,4,8$ (proved by Hurwitz), for which we have bilinear identities, for $n=16$ the Pfister's identity is not bilinear.
I was wondering if there is any resolution for this issue and we have any weaker version of the argument I alluded to above. More importantly, is that the only reason we could not step forward and extend the idea? Or there are some other issues, which we cannot resolve for even case $n=16$ and possibly $n=8$!