If we have a map from a Riemannian manifold $(M,g)$ let's call it $\phi$ to reals, with $x^i$ the local coordinates on $M$. The curved Laplace equation (Laplace-Beltrami) looks like $$ g^{ij}\left(\phi_{ij}-\Gamma^k{}_{ij}\phi_{k}\right)=0 $$ Where $\Gamma^k_{ij}$ are Christoffel symbols of $g$ and I denote derivatives by subscripts. Looking at the structure we can notice that adding to $\Gamma^k_{ij}$ a tensor let's say $S^k_{ij}$ which has the property that the vector $g^{ij}S^k{}_{ij}=0$, this change doesn't affect the equation. This would correspond to adding a metric $h$ to $g$ such that $S$ are Christoffel symbols of $h$.
This would for example mean there are more Riemannian spaces $(M,g)$ for which the Laplace equation looks the same, also there are Riemannian spaces $(M,g)$ for which the Laplace equation looks exactly like a flat one.
First of all, I would like to ask if there is some mistake in my reasoning? Secondly, I would like to ask for some references where this is explored, I tried googling but I didn't quite know what to write and that what I tried didn't quite find any papers or something like that.
Maybe this could be used to find the simplest metric to a particular solution.
Thank you very much.