Clavelli-Shapiro Trick

Computing the expected energy for a harmonic oscillator in a heat bath isn't bad. Other observables aren't too bad either, but here is a trick that will generalize.

First introduce new b operators obeying the same commutation relations.

To compute the thermal expectation value with parameter y ( fugacity goes with chemical potential, I mean the one that corresponds to temperature), you take the regular expectation value of your operator, but replace all the a's and their adjoints by

respectively.

3 Proofs in A.6 at least in a useful special case.

This seems insanely useful depending on how far it is pushed.