Let  $X$  be  a  non vanishing real analytic vector  field  on an open  manifold $M$. What  kind  of  obstructions would  appear when we search for  a Riemannian metric  on $M$ such that the space  of  harmonic  functions would  be  invariant under the derivation operator $f \mapsto X.f$?