Quantum systems with real structure

A lot has been said as to why quantum mechanics needs complex numbers.

However, all measurements produce real values. Expectation values are real, the observables form a real Lie-algebra (use $-\mathrm i[·,·]$ as Lie bracket instead of the raw commutator) and any state $\mathbf z$ is uniquely determined by the real-valued function $P_{\mathbf z}(·) = |\langle \mathbf z,·\rangle|^2$ on the space of states.

The space of states itself, however, is a complex projective space, and even though phase factors are irrelevant when taking expectation values, they do matter when combining states - linear superposition $α\mathbf z + β\mathbf z'$ depends on the ration $|α|/|β|$ as well as the phase difference $\arg(α) - \arg(β)$.

But there is at least one important quantum system which can be modeled as a real system: the qubit, ie the space of spin states of a single electron.

In addition to the standard representation as $\mathbf z=(z_1,z_2)\in\mathbb C^2,\langle \mathbf z,\mathbf z\rangle =1$, states can also be represented as $\vec z\in\mathbb R^3,\vec z^2=1$ via $$ \vec z = \begin{pmatrix} \langle\mathbf z,\pmb σ_1\mathbf z\rangle\\ \langle\mathbf z,\pmb σ_2\mathbf z\rangle\\ \langle\mathbf z,\pmb σ_3\mathbf z\rangle \end{pmatrix} = \begin{pmatrix} 2\Re(\bar z_1z_2)\\ 2\Im(\bar z_1z_2)\\ |z_1|^2 - |z_2|^2 \end{pmatrix} $$ where $\pmb σ_1, \pmb σ_2, \pmb σ_3$ are the Pauli-matrices.

This is the well-known Hopf fibration.

What makes this useful is the fact that $P_z$ takes a particular simple form: $$ P_{\mathbf z}(\mathbf z') = \frac 12 (1 + \vec z·\vec z') = \frac 12 (1 + \cos\measuredangle(\vec z,\vec z')) $$

This also means that states $\mathbf z,\mathbf z'$ are orthogonal in $\mathbb C^2$ exactly when $\measuredangle(\vec z,\vec z')=180°$ in $\mathbb R^3$, ie $\vec z'=-\vec z$

The expection value of any observable $\mathbf A$ can be expanded in terms of $P_{\mathbf a_i}$ after a choice of orthonormal eigenvectors $\{\mathbf a_i\}$, ie $$ \langle\mathbf z,\mathbf{Az}\rangle = \sum_i α_iP_{\mathbf a_i}(\mathbf z) $$ where $\{α_i\}$ denote the eigenvalues.

Substituting our definition of $P_{\mathbf a_i}$ and using the fact that $\mathbf a_1=-\mathbf a_2$, this is simply $$ \langle\mathbf z,\mathbf{Az}\rangle = \frac 12(α_1 + α_2) + \frac 12(α_1-α_2)\,\vec a_1·\vec z $$

Restricting ourselves to $\mathfrak{su}(2)$, ie the traceless matrices with $α_2=-α_1$, yields $$ \langle\mathbf z,\mathbf{Az}\rangle = \vec A·\vec z $$ where $\vec A=α_1\,\vec a_1$ is the representation of our observable.

We arrived at a quantum system which can be described with nothing more than school math.

After this lengthy exposition, finally my question: Are there any other quantum system with a similar real structure? If not, is there a particular reason why?

Update

To get a better understanding of the problem, I read up on complex manifolds, and it seems that the qubit is indeed special.

In particular, the projective spaces $P_n\mathbb C$ are complex manifolds, but spheres $S^k$ do not even permit almost-complex structures for $k\not= 2,6$ (Borel, Serre 1951) and are out as state spaces in the general case.

The finite-dimensional complex projective spaces can be realized as different homogeneous spaces. In particular, we have the obvious $$P_n\mathbb C \cong \mathbb C^{n+1}\setminus\{0\}\:\big/\:\mathbb C^*$$ the not-so-obvious $$P_n\mathbb C \cong U(n+1)\:\big/\:U(n)\times U(1)$$ and the one which is probably the most clear-cut $$P_n\mathbb C \cong S^{2n+1}\:\big/\:S^1$$ but does not convey that we are dealing with complex spaces. This quotient also only ever results in a sphere for $n=1$, ie our qubit case.

The Wikipedia page on Hopf fibration links to this paper, where the fibration $$S^7\xrightarrow{S^3}S^4$$ is used to model the state space $P_3\mathbb C$ of two qubits. Even though there might be some insight into the structure of the state space in this fibration, it is less natural (and less useful, for that matter) than the single qubit Hopf fibration $$S^3\xrightarrow{S^1}S^2$$ where actually $S^2 \cong P_1\mathbb C$, whereas obviously $S^4\not\cong P_3\mathbb C$ as they do not even have the same dimensions when considered as real manifolds.