Why Hubbard model cannot be solvable exactly?

Question

The Hubbard model is a model to describe electrons in a lattice. In general, the Hubbard model Hamiltonian $H$ contains two terms:

1. The kinetic term:$$T=-t\sum_{\langle ij\rangle\sigma} [c_{i\sigma}^\dagger c_{j\sigma} + h.c.] $$
2. The onsite Coulomb interaction term:

$$U=u\sum_{i=1}^N n_{i\uparrow}n_{i\downarrow}$$

So my question is: why can we not diagonalize the Hamiltonian: $H=T+U$? Some books attribute the reason that $T$ doesn't commute with $U$, therefore we need to formulate a perturbation theory. In particular, I want to know whether the space of solution of $T$ has the same dimension compared to the space of solution of $U$ due to $[T,U] \neq 0$?

Edit:

For $T$ operator we have the following eigenequation: $$T|n\rangle=T_n|n\rangle$$ For $U$ operator we have another one eigenequation: $$U|\alpha\rangle=U_\alpha|\alpha\rangle$$ Due to $[T,U] \neq 0$, I am wandering whether $|n\rangle$ has the same dimension compared to $|\alpha\rangle$? And why $$[T+U]|N\rangle \overset{?}{=} H_N|N\rangle.$$

Answer 1

As for the question of dimensionality: Yes, the two sets of solutions for $T$ and $H$ have the same dimensions (Because they act on the same Hilbert space).

The pedestrian reason why you cannot solve the Hubbard model analytically is that the $U$ term contains quartic interactions ($c^\dagger c^\dagger c c$), and it is only in special cases that these can be diagonalized exactly.