Let $M$ be an $n$-dimensional complete Riemannian manifold and $r$ is the distance function to a fixed point.
The Hessian comparison theorem says that if the sectional curvature of $M$ is bounded (precisely $k\le \operatorname{sec}\le K$), then the Hessian of $r$ is bounded by the Hessians of the distance function for the space form with constant sectional curvature $k$ and $K$ (precisely $\operatorname{Hess}_Kr\le\operatorname{Hess}r\le \operatorname{Hess}_kr$).
My question: What is the difference between the Hessian of $r$ and the Hessian of the distance function for the space form with constant sectional curvature 0? It's better to have a Taylor expansion.
Following the notation in Petersen's Riemannian geometry, the Hessian of the distance function for the space form with constant sectional curvature $k$ is $\operatorname{Hess}_kr=\frac{\operatorname{sn}_k'(r)}{\operatorname{sn}_k(r)}g_r$ where $g_r=g-(\operatorname{d}r)^2$ and $\frac{\operatorname{sn}_k'(r)}{\operatorname{sn}_k(r)}=\frac{1}{r}-\frac{k}{3}r+O(r^2)$. Note that $\operatorname{Hess}_0r=\frac{1}{r}g_r$. Based on the Hessian comparison theorem, I guess that $\operatorname{Hess}r-\operatorname{Hess}_0r=-\frac{\operatorname{sec}}{3}rg_r+O(r^2)$.
I originally guessed that the difference $\operatorname{Hess}r-\operatorname{Hess}_0r$ is a multiplication of the Ricci tensor, big thanks to MySheperd who excluded this possibility and suggested that the difference is $O(r)$ in the answer below. Then I started to ask the above question.
Thank you!