The first and second derivatives of the distance function (either the full $d:M\times M\to \mathbb{R}$ function or the $d(p,\cdot):M\to \mathbb{R}$ function) as well as the derivative of the exponential map (again both of the full $\exp:TM\to M$ map and of the map $\exp\_p:T\_pM\to M$) may be computed with the aid of Jacobi fields, i.e, solutions to Jacobi's equation.


I have a scenario where I need second derivatives of the `full' exponential map $\exp:TM\to M$. That is, denoting the pushforward of a differentiable map by a '$ \_\*$', I would like to know a thing or two about $\nabla\_X\exp\_\*\mathcal{V}$ (where $\mathcal{V}$ is a section of $TTM$ and $X$ is an appropriate vector field). In particular, I think I will require some comparison techniques analagous to those for Jacobi fields (e.g., Rauch's comparison theorems).


Can anyone point me in the right direction?



Higher derivatives than Jacobi fields
The higher derivatives of the exponential map satisfy the corresponding higher derivative of the Jacobi equation (because the first derivative satisfies the Jacobi equation itself), which is just an inhomogeneous Jacobi equation, where the homogeneous part is the original Jacobi equation, and the inhomogeneous term involves lower order covariant derivatives of the Jacobi field and covariant derivatives of the curvature tensor. So you would proceed recursively, bootstrapping pointwise bounds on lower derivatives, as well as pointwise bounds on the curvature tensor and its covariant derivatives, into a pointwise bound of the derivative of the Jacobi field. You'll need to figure out how get pointwise bounds for a solution to an inhomogeneous self-adjoint linear second order ODE. I'm sure this has been done before, probably for exactly the same purpose as here, but I don't know or remember where.


