Determinant of a propagator Say I have a path integral $\int D \phi \exp(i S_0)$. $S_0$ is the usual free action $$S_0=\frac{1}{2}\int\phi (-\Box-m^2) \phi=\frac{1}{2}\int \phi G^{-1} \phi,$$ 
and at the moment I'm not interested in interactions.
When I do this integration, I just get a term like $(\det G)^{-1/2}$, which is the zero point energy. Normally I don't care about it and it is absorbed in the normalization of the path integral.
But say that this propagator depends on a parameter, and I would like to actually keep this determinant. How do I actually compute it?
I know I can do the trick
$$(\det G)^{-1/2}=\exp(-\frac{1}{2}\text{Tr} \log G) $$
but then I would have to compute something like (after Wick rotating and going to momentum space)
$$\int \frac{d^d k}{(2 \pi)^d} \log(\frac{1}{k^2+m^2})$$
which of course diverges (because I have infinitely many harmonic oscillators) and requires some kind of renormalization. Hence my question: how do you compute such a determinant?
 A: You can simply compute the integral using your preferred regularization method (cut-off, dimensional, Pauli-Villars...), and if all goes well (which is not guaranteed), the divergences will not depend on your parameter and they will eventually disappear when you compute physical stuff. If this does not happen, maybe your theory is simply ill-defined.
And as for adding an extra factor of $k^2$ as mentioned in the comments, we know that adding the trace of $\ln k^2$ would add a (for physical applications) irrelevant constant (which in dimensional regularization is furthermore zero), such that you can add it whenever you feel like it.
A: As far as I know, the easiest way to do the integration by hand is for integer dimensions (i.e. Mathematica gives you the integral in terms of hypergeometric functions, but that's not really helpful). You can derive the integral (with a hard cut-off $\Lambda$) with respect to $m^2$, perform the integral, expand in $m/\Lambda$, and then integrate back. It might be helpful to subtract a constant $\int_k \ln k^2$. 
In the case $d=3$, you should get $\frac{\Lambda}{2\pi^2}m^2-\frac{m^3}{6\pi}$. Note that the first term is not (necessarily) problematic, and can have a very physical interpretation. For example, in the context of $N$ relativistic bosons in the limit $N\to\infty$ (a well known model in condensed matter), this term correspond to the renormalization of the parameter driving the system through the phase transition between the ordered and disordered phase. 
