Suppose $M$ is a Riemannian manifold. Consider flow $\frac{\partial}{\partial t}g_{ij}=-2(R_{ij}+\nabla_i \nabla_j f)$, where $f$ is a time-dependent function. I would like to prove that flows of this form are equivalent, up to diffeomorphism, to the Ricci flow $\frac{\partial}{\partial t}g_{ij}=-2R_{ij}$, that is:
By defining a 1-parameter family of diffeomorphism $\Psi(t):M\to M$ by
$$\frac{d}{dt}\Psi(t)=\nabla_{g(t)}f(t),$$ $$\Psi(0)=id_M$$
I want to show that $\bar{g}(t):=\Psi(t)^*g(t)$ satisfy $$\frac{\partial}{\partial t}\bar{g}_{ij}=-2 \bar{R}_{ij}.$$
My problem is that I don't know how to calculate $\frac{\partial}{\partial t}\Psi(t)^*g(t)$. I know that $\frac{\partial}{\partial t}\Psi(t)^* \alpha=\mathcal{L}_{\nabla f} \alpha$, where $\mathcal{L}$ is Lie derivative and $\alpha$ is a time-independent object, but I faced problem when $\alpha$ is a time-dependent object.
Can someone point me in the right direction? Thanks in advance for your time.