Harmonic $Sp(2)$-invariant $G_2$-structures on the 7-sphere

We describe the $10$-dimensional space of $Sp(2)$-invariant $G_2$-structures on the homogeneous $7$-sphere $S^7=Sp(2)/Sp(1)$ as $\mathbb{R}^+\times Gl^+(3,\mathbb{R})$. In those terms, we formulate a general Ansatz for $G_2$-structures, which realises representatives in each of the $7$ possible isom