Given any collection of complementary, oriented (2D) planes in n-dimensional space, and an angle associated with each one, there is a unique rotation of the whole space which restricts to rotations in each plane by their given angles. Conversely, it's not impossible to prove all rotations of space arise in this way. Which leads to the question: can a single rotation arise in this way but from two different sets of complementary planes? If the attached angles are all distinct, I think the answer is no, which means the collection of stable planes is an invariant of "almost all" rotations (in the sense that the set of exceptions are of positive codimension in the group of all rotations).
If any two of the attached angles are the same, though, then there are an infinite number of planes which are stable under the rotation. Indeed, it is easy to verify that such a rotation moves every ray from the origin by the same angle. (Compare with rotation around an axis in 3D: if a ray is perpendicular to the axis then it is rotated by the full angle, but otherwise if the ray is closer to the axis, the before and after rays will be less than the full angle apart from each other.)
For simplicity, lets just take n=4 so our isoclinic rotations are rotations in two orthogonal planes by the same angle. Depending on if combining the orientations of the planes agrees with the orientation of the whole space or not, we can call them left isoclinic or right isoclinic rotations, because they correspond to multiplying quaternions (which form a 4-dimensional space) either on the left or on the right by unit quaternions. As such, the subgroups of left and right isoclinic rotations are both copies of S^3 within SO(4).
Now my curiosity is in which planes are stable under an isoclinic rotation, or without loss of generality a one-parameter subgroup of them. Any 2D plane has a unique orthogonal complement, and to the two resulting planes we get a one-parameter subgroup, which can be identified with a line in the corresponding lie algbera so(3) of dimension 3. Therefore we have some map from the Grassmanian manifold of 2D planes to the projective plane. Characterizing stable planes of a given one-parameter group of isoclinic rotations is then equivalent to finding the fibers of this map, which should all be identical by symmetery. So we have some fiber bundle
$$F\to \mathrm{Gr}(2,4)\to\mathbb{P}^2 $$
Is there a good explicit description of the fibers? Or just its isomorphism type?
My initial thoughts: Given the fiber $F\subseteq\mathrm{Gr}(2,4)$ of a point in $\mathbb{P}^2$ corresponding to some one-parameter subgroup $H$ of isoclinic rotations, if we have any point $x\in\mathbb{S}^3\subset\mathbb{R}^4$ the orbit $Hx$ is a circle in $\mathbb{S}^3$, which determines a plane $\pi$ in $\mathbb{R}^4$, and indeed $\pi\in F$. So we get a map $\mathbb{S}^3\to F$, and the fibers of this map should be the circles in $\mathbb{S}^3$ which intersect planes $\pi\in F$. Thus we have a fiber bundle $\mathbb{S}^1\to\mathbb{S}^3\to F$, which naively suggests to me $F\cong\mathbb{S}^2$ based on the Hopf fibration.
Some folks in chat though suggested $F$ should be $\mathbb{P}^2$ or a couple disjoint copies of it. Also mentioned was the fact that $\widetilde{\mathrm{Gr}}(2,4)$, the oriented Grassmanian which is a double cover of the usual $\mathrm{Gr}(2,4)$, is $\cong \mathbb{S}^2\times\mathbb{S}^2$, although I don't know why or for sure how to use this.
Ideas?
$\newcommand{\Reals}{\mathbf{R}}\newcommand{\Cpx}{\mathbf{C}}$Consider a consistently-oriented pair of planes, regarded as complex lines. In fact, you may as well fix an identification of $\Reals^{4}$ with $\Cpx^{2}$ and take the complex coordinate axes as $2$-planes. A left isoclinic rotation in these planes is precisely multiplication by a complex scalar, so it stabilizes every complex line through the origin, i.e., an $S^{2}$.
Alternatively, the standard complex structure $J$ is the left isoclinic rotation by $\pi/2$, and commutes with every left isoclinic rotation about the complex coordinate axes.
Particularly, this picture makes apparent that a left isoclinic rotation by some angle $\theta$, namely scalar multiplication by $e^{i\theta}$, is the same for every pair of orthogonal complex lines.
(Incidentally, rotation by a half-turn (a.k.a. $-I$) obviously stabilizes every plane through the origin, and this is the only non-trivial rotation to stabilize arbitrary $2$-planes. I'm guessing this was too trivial to mention.:)