Se pueden derivar las propiedades de los cuaternarios a través del álgebra de Clifford y el producto geométrico de los vectores.
El producto geométrico funciona así:
$$e_a e_b = \begin {cases} 1, & a = b \\ -e_b e_a, & a \neq b \end {cases}$$
donde $a, b$ puede ser $x$ , $y$ o $z$ como siempre. Esto captura tanto el trabajo del producto cruzado como el producto puntual en un producto de vectores base.
Entonces puede identificar
$$ \begin {align*} i &= -e_y e_z \\ j &= -e_z e_x \\ k &= -e_x e_y \end {align*}$$
Y entonces las propiedades de los cuaterniones siguen naturalmente.
$$i^2 = (-e_y e_z)(-e_y e_z) = (e_y e_z)(e_y e_z) = -e_y (e_z e_z) e_y = -e_y e_y = -1$$
Y de manera similar para $j^2$ y $k^2$ así como el $ijk$ producto:
$$ijk = (-e_y e_z)(-e_z e_x)(-e_x e_y) = - e_y e_z e_z e_x e_x e_y = -e_y (e_z e_z)(e_x e_x) e_y = -e_y e_y = -1$$
Esto te permite interpretar cuaterniones en un muy geométrico manera: el $i,j,k$ hacer no representan vectores, sino más bien planos orientados. Es sólo que en 3d cada plano tiene un único vector normal, por lo que a menudo abusamos de esta dualidad.