[step:Show that the Jacobi field is index-orthogonal to endpoint-vanishing fields]Let $g$ denote the Riemannian metric on $M$, let $R$ denote the Riemann curvature tensor with convention
\begin{align*}
R(X,Y)Z=\nabla_X\nabla_YZ-\nabla_Y\nabla_XZ-\nabla_{[X,Y]}Z,
\end{align*}
and let $\mathcal{L}^1$ denote one-dimensional [Lebesgue measure](/page/Lebesgue%20Measure) on $\mathbb{R}$. For continuous piecewise smooth vector fields $X,Y:[a,b]\to TM$ along $\gamma$, meaning $X(t),Y(t)\in T_{\gamma(t)}M$, define the index form by
\begin{align*}
I_\gamma[X,Y]
:=
\int_a^b
\left(
g_{\gamma(t)}(\nabla_t X(t),\nabla_t Y(t))
-
g_{\gamma(t)}(R(X(t),\dot{\gamma}(t))\dot{\gamma}(t),Y(t))
\right)
\,d\mathcal{L}^1(t),
\end{align*}
where $\nabla_t X$ denotes covariant differentiation of $X$ along $\gamma$ and $\dot{\gamma}$ denotes the velocity field of $\gamma$. We prove that
\begin{align*}
I_\gamma[J,W]=0.
\end{align*}
Let $a=t_0<t_1<\cdots<t_m=b$ be a partition such that $W$ is smooth on each closed subinterval $[t_{k-1},t_k]$. On each such subinterval, metric compatibility of the Levi-Civita connection gives
\begin{align*}
\frac{d}{dt}g_{\gamma(t)}(\nabla_t J(t),W(t))
=
g_{\gamma(t)}(\nabla_t\nabla_t J(t),W(t))
+
g_{\gamma(t)}(\nabla_t J(t),\nabla_t W(t)).
\end{align*}
Therefore,
\begin{align*}
\int_{t_{k-1}}^{t_k}
g_{\gamma(t)}(\nabla_t J(t),\nabla_t W(t))
\,d\mathcal{L}^1(t)
&=
g_{\gamma(t_k)}(\nabla_t J(t_k),W(t_k))
-
g_{\gamma(t_{k-1})}(\nabla_t J(t_{k-1}),W(t_{k-1}))\\
&\quad -
\int_{t_{k-1}}^{t_k}
g_{\gamma(t)}(\nabla_t\nabla_t J(t),W(t))
\,d\mathcal{L}^1(t).
\end{align*}
Summing over $k=1,\dots,m$, the interior boundary terms telescope because $J$ is smooth and $W$ is continuous at the break points. Since $W$ is continuous at each interior break point and $W(a)=W(b)=0$, the endpoint terms vanish. Hence
\begin{align*}
\int_a^b
g_{\gamma(t)}(\nabla_t J(t),\nabla_t W(t))
\,d\mathcal{L}^1(t)
=
-
\int_a^b
g_{\gamma(t)}(\nabla_t\nabla_t J(t),W(t))
\,d\mathcal{L}^1(t).
\end{align*}
Because $J$ is a Jacobi field along $\gamma$, it satisfies
\begin{align*}
\nabla_t\nabla_t J(t)+R(J(t),\dot{\gamma}(t))\dot{\gamma}(t)=0.
\end{align*}
Substituting this identity into the definition of the index form gives
\begin{align*}
I_\gamma[J,W]
&=
\int_a^b
\left(
g_{\gamma(t)}(\nabla_t J(t),\nabla_t W(t))
-
g_{\gamma(t)}(R(J(t),\dot{\gamma}(t))\dot{\gamma}(t),W(t))
\right)
\,d\mathcal{L}^1(t)\\
&=
-\int_a^b
g_{\gamma(t)}(\nabla_t\nabla_t J(t),W(t))
\,d\mathcal{L}^1(t)
-
\int_a^b
g_{\gamma(t)}(R(J(t),\dot{\gamma}(t))\dot{\gamma}(t),W(t))
\,d\mathcal{L}^1(t)\\
&=
-\int_a^b
g_{\gamma(t)}(\nabla_t\nabla_t J(t)+R(J(t),\dot{\gamma}(t))\dot{\gamma}(t),W(t))
\,d\mathcal{L}^1(t)\\
&=0.
\end{align*}[/step]