I am currently studying functional integration in the context of classical and quantum equilibrium thermodynamics. However one thing puzzles me:
In the book "Phase Transitions and Renormalization Group" by Zinn-Justin he explains several times, for example in 8.9 that corrections to a Gaussian model described by the partition function $$ \mathcal{Z}=\int \Pi_id\sigma_i\;\exp\left(-H_0(\sigma)-b_4 \sum_i\sigma_i^4+\sum_iH_i\sigma_i \right)\tag{p.203} $$ can be either derived by expanding the quartic term like this: $$ \mathcal{Z}=\int \Pi_id\sigma_i\;(\sum \frac{(-1)^l}{l!}(b_4 \sum_i\sigma_i^4)^l)\;\exp\left(-H_0(\sigma)+\sum_iH_i\sigma_i \right)\tag{p.203} $$ which allows you to compute the vertex function at vanishing magnetic field $\Gamma^{(2)}(M=0)$, in other words the heat capacity perturbatively (see Eq.8.75). You can do this for example until $l=1$ which is first order perturbation theory. Alternatively, he shows that the same result can be obtained by the steepest descent method which gives you directly the generating functional $$ \Gamma(M)=\Gamma(M)_{\text{Gaussian}}+\frac{1}{2}\text{tr}\log(\frac{\partial^2H(M)}{\partial M\partial M}).\tag{6.31} $$ The second derivative with respect to $M$ then gives the same result for the heat capacity. I take the message that steepest descent gives the same result as first order perturbation theory.
My question is: Is this always the case? In general, I thought the saddle point approximation is accurate when the pre-factor of the action in the exponent is large, and perturbation theory is the way to go when the parameter in front of the quartic term is small. In this case here, I understand that when the parameter $b_4$ is small, the distribution is close to Gaussian so I see why the saddle point approximation gives the same result as simple perturbation theory. So, here this is the same limit. But this is not a general result right? When are these two approaches a different limit?