Taylor expansion of an integral (for Thermodynamic Perturbations)

[mSSM]

In thermodynamic perturbation theory (chapter 32 in Landau's Statistical Physics) for the Gibbs (= canonical) distribution, we have E = E_0 + V, where V is the perturbation of our energy.

When we want to calculate the free energy, we have:
e^{-F/T} = \int e^{-(E+V)/T} \mathrm{d}\Gamma

We can expand the expression inside the integral for small values of V to second order we obtain:
e^{-F/T} = \int e^{E/T} (1 - V/T + (V/T)^2) \mathrm{d}\Gamma

Now, the next step is to logarithmize this expression and then expand it in terms of V again, to second order:
F = -T \log \bigl( \int e^{E/T} (1 - V/T + (V/T)^2) \mathrm{d}\Gamma \bigr)

But now I am stuck. I don't understand how to expand this in a series. The zeroth order still works great, and I obtain -T \log \int e^{E/T} = F, which is just what I want. But from the first order onwards I am stuck:
What is the rule for Taylor expanding an integral? Especially, e.g., if I have f'(a) (x-a) for first order, where exactly do I 'put' the x? Inside the integral?

 

[mSSM]

So, I thought about this a little with a friend. I believe that I could see the integral as a sum, with a lot of different V_i, albeit them lying infinitesimally close to each other. Therefore, one could view the expansion with respect to V as expansions with respect to a lot of different V_i in this 'sum'. After this expansion one could imagine making the sum into an integral once more, such that the orders of V are inside the integrals.