Let $(M,g)$ be some smooth, Riemannian manifold. Let $d$ be the exterior derivative and $\delta$ the codifferential on forms. For a smooth vector field $X$, let $L\_X$ be the Lie derivative associated to $X$. We know from Cartan formula that $L\_X = d \iota\_X + \iota\_X d$ where $\iota\_X$ is the interior derivative associated to the vector field $X$. So it is well-known that $L\_X$ and $d$ commute: for any arbitrary form $\omega$, we have that $L\_Xd\omega = dL\_X\omega$. 


This is, of course, not true for codifferentials. In general $[L\_X,\delta]\neq 0$. For certain cases the answer is well known: if $X$ is a Killing vector ($L\_Xg = 0$) then since it leaves the metric structure in variant, it commutes with the Hodge star operator, and so $L\_X$ commutes with $\delta$. Another useful case is when $X$ is conformally Killing with constant conformal factor ($L\_X g = k g$ with $dk = 0$). In this case conformality implies that the commutator $[L\_X,\*] = k^\alpha \*$ where $\alpha$ is some numerical power depending on the rank of the form it is acting on (I think... correct me if I am wrong), so we have that $[L\_X,\delta] \propto \delta$. 


So my question is: "Is there a general nice formula for the commutator $[L\_X,\delta]$?" If it is written down somewhere, a reference will be helpful. (In the Riemannian case, by working with suitable symmetrisations of metric connection one can get a fairly ugly answer by doing something like $\delta \omega \propto \mathop{tr}\_{g^{-1}} \nabla\omega$ and use that the commutators $[L\_X, g^{-1}]$ and $[L\_X,\nabla]$ are fairly well known [the latter giving a second-order deformation tensor measuring affine-Killingness]. But this formula is the same for the divergence of arbitrary covariant tensors. I am wondering if there is a better formula for forms in particular.) 


A simple explanation of why what I am asking is idiotic would also be welcome. 



Commutator of Lie derivative and codifferential?
Doing along what Deane and Jose suggest one can make the following observations. First is that the Hodge star operator can be expressed in terms of the inverse metric $g^{-1}$ and volume form $\epsilon$. That is to say, formally for a form $\omega$ we have
$$ \*\omega = \epsilon \cdot (g^{-1}) \cdot \omega $$


Take $L\_X$, when it falls on $\epsilon$ we have by definition of divergence
$$ L\_X \epsilon = (\mathrm{div} X) \epsilon.$$ 
Where $\mathrm{div} X$ is also half the metric trace of the deformation tensor $L\_X g$. We also have the formula
$$ L\_X g^{-1} = g^{-1}(L\_X g) g^{-1} .$$
Plugging this in we have that schematically 
$$ L\_X \*\omega = (\mathrm{div} X) \*\omega + \* [(L\_X g)\cdot \omega] + \* L\_X \omega$$
from which the computation can be completed. 


