Appendix 2 — Derivation of the Derivative of the Christoffel Symbols in General Form
It is shown how the Christoffel symbol depends solely on the components of the metric tensor and its derivatives. This is particularly useful for implementation in spreadsheets or computer programs.
Christoffel symbol
\begin{equation}
\begin{aligned}
\Gamma^{\rho}{}_{\mu\nu}
=
\frac{1}{2} g^{\rho\alpha}
\left(
\frac{\partial g_{\nu\alpha}}{\partial x^{\mu}}
+
\frac{\partial g_{\mu\alpha}}{\partial x^{\nu}}
-
\frac{\partial g_{\mu\nu}}{\partial x^{\alpha}}
\right)
\end{aligned}
\end{equation}
Derivative of the Christoffel symbol
\begin{equation}
\begin{aligned}
\frac{\partial \Gamma^{\rho}{}_{\mu\nu}}{\partial x^{\gamma}}
&=
\frac{1}{2}
\frac{\partial g^{\rho\alpha}}{\partial x^{\gamma}}
\left(
\frac{\partial g_{\nu\alpha}}{\partial x^{\mu}}
+
\frac{\partial g_{\mu\alpha}}{\partial x^{\nu}}
-
\frac{\partial g_{\mu\nu}}{\partial x^{\alpha}}
\right)
\\
&\quad
+
\frac{1}{2} g^{\rho\alpha}
\left(
\frac{\partial^{2} g_{\nu\alpha}}{\partial x^{\mu}\partial x^{\gamma}}
+
\frac{\partial^{2} g_{\mu\alpha}}{\partial x^{\nu}\partial x^{\gamma}}
-
\frac{\partial^{2} g_{\mu\nu}}{\partial x^{\alpha}\partial x^{\gamma}}
\right)
\end{aligned}
\end{equation}
Derivative of the inverse metric
\begin{equation}
\begin{aligned}
\frac{\partial g^{\rho\alpha}}{\partial x^{\gamma}}
=
\frac{\partial \frac{1}{g_{\rho\alpha}}}{\partial x^{\gamma}}
=
\frac{-1}{g^2_{\rho\alpha}}\cdot\frac{\partial g_{\rho\alpha}}{\partial x^{\lambda}}
=
- (g^{\rho\alpha})^{2}
\frac{\partial g_{\rho\alpha}}{\partial x^{\gamma}}
\end{aligned}
\end{equation}
Full derivative of the Christoffel symbol
\begin{equation}
\begin{aligned}
\frac{\partial \Gamma^{\rho}_{\mu\nu}}{\partial x^{\gamma}}
&=
-\frac{1}{2}
(g^{\rho\alpha})^{2}
\frac{\partial g_{\rho\alpha}}{\partial x^{\gamma}}
\left(
\frac{\partial g_{\nu\alpha}}{\partial x^{\mu}}
+
\frac{\partial g_{\mu\alpha}}{\partial x^{\nu}}
-
\frac{\partial g_{\mu\nu}}{\partial x^{\alpha}}
\right)
\\
&\quad
+
\frac{1}{2} g^{\rho\alpha}
\left(
\frac{\partial^{2} g_{\nu\alpha}}{\partial x^{\mu}\partial x^{\gamma}}
+
\frac{\partial^{2} g_{\mu\alpha}}{\partial x^{\nu}\partial x^{\gamma}}
-
\frac{\partial^{2} g_{\mu\nu}}{\partial x^{\alpha}\partial x^{\gamma}}
\right)
\end{aligned}
\end{equation}
Compact form
\begin{equation}
\begin{aligned}
\frac{\partial \Gamma^{\rho}_{\mu\nu}}{\partial x^{\gamma}}
=
- g^{\rho\alpha}
\frac{\partial g_{\rho\alpha}}{\partial x^{\gamma}}
\Gamma^{\rho}_{\mu\nu}
+
\frac{1}{2} g^{\rho\alpha}
\left(
\frac{\partial^{2} g_{\nu\alpha}}{\partial x^{\mu}\partial x^{\gamma}}
+
\frac{\partial^{2} g_{\mu\alpha}}{\partial x^{\nu}\partial x^{\gamma}}
-
\frac{\partial^{2} g_{\mu\nu}}{\partial x^{\alpha}\partial x^{\gamma}}
\right)
\end{aligned}
\end{equation}