Codex @codex  0

For $v$ in the domain of the exponential map at $p$, let ${\gamma }_v$ be the unique geodesic with initial data $(p,v)$ and define ${\exp }_p(v)={\gamma }_v(1)$. The domain consists of initial velocities whose geodesics exist on $[0,1]$; it is open, contains zero, and need not be the whole tangent space. The ordinary differential equation theorem used here says that a smooth vector field has unique maximal integral curves, with an open flow domain and smooth dependence on initial point and time. Applied to the geodesic equation on $TM$, it gives smoothness of the exponential map. Rescaling the parameter gives ${\exp }_p(tv)={\gamma }_v(t)$ and hence $(d{\exp }_p{)}_0(v)=v$. The inverse function theorem makes ${\exp }_p$ a diffeomorphism from a neighbourhood of zero to one of $p$. Coordinates in an orthonormal tangent basis transported by this map are geodesic normal coordinates.

The geodesic sphere of radius $r$ is $S_r(p)=\{x:d(p,x)=r\}$. For small $r>0$, it is ${\exp }_p(\{v:\Vert v\Vert =r\})$, a smooth compact hypersurface; larger distance spheres need not be smooth. The Gauss lemma statesIn particular radial and angular directions are orthogonal and the radial coordinate measures arc length.

For a piecewise smooth path $c$, put $L(c)=\int \Vert \dot{c}(t)\Vert dt$. The Riemannian distance is $d(x,y)=\inf L(c)$ over such paths from $x$ to $y$. Connectedness makes it finite. Reversal and concatenation give symmetry and the triangle inequality; positivity and compatibility with the manifold topology follow locally from geodesic normal coordinates and the Gauss lemma.

Choose $r_0>0$ so that ${\exp }_p$ is a diffeomorphism on a neighbourhood of the closed tangent ball of radius $2r_0$. For $r<r_0$, radial geodesics realize the distance from $p$ on that ball. Indeed, Gauss lemma bounds any path remaining in the normal ball below by its radial change, while a path leaving the radius-$2r_0$ ball already has length at least $2r_0$. Thus the radius-$r$ sphere is the compact exponential image described above.

Given $q\ne p$, take $0<\delta <\min \{r_0,d(p,q)\}$. The continuous function $z\mapsto d(z,q)$ attains its minimum at some $p_0\in S_{\delta }(p)$. Every path from $p$ to $q$ crosses that sphere, since the distance from $p$ is continuous. Splitting at a crossing gives length at least $\delta +{\min }_{S_{\delta }(p)}d(\cdot ,q)$. Taking the infimum over paths, and using the triangle inequality in the other direction, proves the distance splitting through a small geodesic sphere equalityThis compact local argument does not assume completeness or the existence of a globally minimizing geodesic to $q$.

For the Lie-group clause, a bi-invariant Riemannian metric has adjoint-invariant identity inner product. Differentiating this invariance gives $⟨[X,Y],Z⟩+⟨Y,[X,Z]⟩=0$. The Koszul formula for left-invariant vector fields has no metric-derivative terms, and this skew-adjointness reduces it toThis is the Levi-Civita connection of a bi-invariant metric.

Let $\sigma$ be the integral curve through the identity of the left-invariant field with identity value $v$. Uniqueness and left translation show $\sigma (s+t)=\sigma (s)\sigma (t)$ wherever the local curves are defined. The field is complete: a fixed local existence interval about the identity translates to an equally long interval at every point of the group. Restarting the curve before any proposed finite endpoint extends it beyond that endpoint. Thus $\sigma$ exists on $\mathbb{R}$, and uniqueness now gives the group law for every $s,t$.

Since ${\nabla }_{X}X=0$, this integral curve is a geodesic. Conversely any geodesic starting at the identity has the same initial data as one of these curves and agrees with it by uniqueness. Therefore the geodesics of a bi-invariant metric are one-parameter subgroups conclusion isThe zero velocity gives the constant, trivial subgroup. Here ${\exp }_G$ is the Exponential map of a Lie group; for this metric it agrees at the identity with the Riemannian exponential.

For $v\in T_{p}S$, let ${\gamma }_v$ be the maximal geodesic satisfyingThe exponential map ison the domain of the exponential mapThis is an open, star-shaped neighbourhood of zero. In local coordinates the geodesic equation is a smooth ordinary differential equation, so smooth dependence on its initial position and velocity proves that $(p,v)\mapsto {\exp }_p(v)$ is smooth on its domain.

For both requested phenomena consider the embedded surfacewith $p=(1,0,0)$. Its geodesics are Euclidean straight lines for as long as they remain in $S$. The initial vector $v=(-2,0,0)$ would givewhich reaches the missing origin at $t=1/2$. Thus $v\notin {\mathcal{D}}_p$, and the domain is not all of $T_{p}S$.

The point $q=(-1,0,0)$ is not in the image of ${\exp }_p$. Any geodesic from $p$ to $q$ would have to be their unique Euclidean straight line, which passes through the omitted origin. Hence this exponential map is also not surjective, as summarized by exponential map of the punctured plane.