The hamiltonian is very long, but the idea is that I have, products of bosons operator, for example: a*a + a*b + a*c + b*b + bb + b*c and so on...
Where a, b an c are boson operators.
As there are "non-diagonal" terms, the idea is use Bogoliubov Transformation for make the hamiltonian a...