Hi everyone,
I use Wick's theorem to decompose expectation values of a string of bosonic creation and annihilation operators evaluated at the vacuum state. This can only be done when the time evolution is driven by a Hamiltonian of the form:
H=\sum_{i,j}{\epsilon_{i,j} c^{\dagger}_{i}c_{j}}...