This is probably a very naive question, but so far I could not find an answer:
Let $(M,g)$ be a Riemannian manifold. Can we always find "orthogonal coordinates" locally?
More precisely, I am asking if for every $p \in M$ there exists a neighbourhood $U$ and a diffeomorphism $\phi:\mathbb{R}^n \to U$, such that $g_{ij}=g(d\phi(e_i),d\phi(e_j))=0$ for $i \neq j$.
Clarification: Note that I want $g_{ij}=0$ on all $U$, not just at $p$. Also, I allow $g_{ii} \neq g_{jj}$ for $i \neq j$ (the special case where $g_{ii}$ is independent of $i$ is called isothermal coordinates-and corresponds to conformal flatness of $U$).
Of course, this is weaker than requiring $M$ to be conformally flat, since a (linear) map which maps an orthogonal basis to an orthogonal basis does not need to be conformal.