Is there a "canonical" quantum system whose complex Hilbert space is $n$-dimensional? (I'm not a physicist.)
It seems from the little I've read that the "canonical" example (or maybe just easiest or first example everyone goes to) of a quantum system whose Hilbert space has 2 dimensions comes from thinking about the spin of a particle; the 2 dimensions comes from the spin-up and spin-down states being a basis.
My question is very simple: Is there a "canonical" (or easy or first) example of a quantum system whose complex Hilbert space is $n$ dimensional for arbitrary $n$?
 A: The set of $n$-dimensional harmonic oscillator states of the type $\{a_k^\dagger\vert 0\rangle, k=1,\ldots,n\}$ span the $n$-dimensional defining representation of $u(n)$.  The elements of $u(n)$ are of the form
$$
C_{ij}=a_i^\dagger a_j\, ,\quad i,j=1,\ldots, n\, .
$$
For $n=2$, one recovers angular momentum by the identification
$$
C_{12}\to L_+\, ,\qquad C_{21}\to L_-\, ,\qquad
\frac{1}{2}\left(C_{11}-C_{22}\right)\to L_z
$$
with $C_{11}+C_{22}$ commuting with all the others when acting on state with fixed total number of excitations.
This is about as simple as it gets.  
Note that, when the total number of excitations is greater than one, v.g. the set of states of the type $a_k^\dagger a_m^\dagger \vert 0\rangle$, with $k,m=1,\ldots,n$, the Hilbert is larger and has the dimensionality of the $(2,0,\ldots,0)$ irrep of $su(n)$.  This is the simple generalization of higher spin states living in a larger Hilbert space.
A: Not particularly.
The reason that spin is a good model is because any nontrivial hamiltonian (and indeed any hermitian operator) can essentially be mapped onto the canonical spin matrix
$$
\sigma_z=\begin{pmatrix}1&0\\0&-1\end{pmatrix}
$$
by adding and multiplying with constants. Moreover, you have a clean basis $\{\sigma_x,\sigma_y,\sigma_y\}$ for representing all nontrivial hermitian operators on that basis, with clear geometrical understanding underpinning all of them ─ but, more importantly, with all of them included in the angular-momentum algebra of the quantum spin. This is why e.g. the Jaynes-Cummings model pretends that its two-level system is a spin, when it's most(?) often an electronic transition between states with a different spatial wavefunction.
This is no longer the case in higher dimensions, because the spin-$j$ algebra is too constrained to match the full generality of hermitian operators on dimension $2j+1$, and indeed it struggles begin with diagonal operators and move on from there.
And as for non-spin systems, there's simply no likely candidates.
