Is there any interacting fermionic model (i.e. with more than quadratic terms) which can be analytically diagonalised? Even the "simple looking" Hubbard model seems to lack an analytical solution. If no, is it even possible to diagonalise an interacting (fermionic) theory?
I'm asking this from a condensed matter perspective but other points of view are more than welcomed (HEP...etc).