a field theory motivated approach to computer algebra

# The Kaluza-Klein example from section 2.5 of hep-th/0701238

This example shows how to use split_index in a somewhat more complicated setting. We first declare the indices that we will use.
{\mu,\nu,\rho,\sigma,\kappa,\lambda,\eta,\chi#}::Indices(full, position=independent); {m,n,p,q,r,s,t,u,v,w,x,y,z,m#}::Indices(subspace, position=independent, parent=full);
$$\displaystyle{}\text{Attached property Indices(position=independent) to }\left(\mu, \nu, \rho, \sigma, \kappa, \lambda, \eta, \chi\#\right).$$
$$\displaystyle{}\text{Attached property Indices(position=independent) to }\left(m, n, p, q, r, s, t, u, v, w, x, y, z, m\#\right).$$
Note the appearance of parent=full. This indicates that the indices in the second set span a subspace of the indices in the first set. Also note that we have declared the indices as position=independent, to prevent Cadabra from automatically raising or lowering them when canonicalising (as this does not really help us here). The remaining declarations are standard,
\partial{#}::PartialDerivative. g_{\mu\nu}::Metric. g^{\mu\nu}::InverseMetric. g_{\mu? \nu?}::Symmetric. g^{\mu? \nu?}::Symmetric. h_{m n}::Metric. h^{m n}::InverseMetric. \delta^{\mu?}_{\nu?}::KroneckerDelta. \delta_{\mu?}^{\nu?}::KroneckerDelta. F_{m n}::AntiSymmetric.
We will want to expand the Riemann tensor in terms of the metric. The following two substitution rules do the conversion from Riemann tensor to Christoffel symbol and from Christoffel symbol to metric (a library with common substitution rules like these is in preparation).
RtoG:= R^{\lambda?}_{\mu?\nu?\kappa?} -> - \partial_{\kappa?}{ \Gamma^{\lambda?}_{\mu?\nu?} } + \partial_{\nu?}{ \Gamma^{\lambda?}_{\mu?\kappa?} } - \Gamma^{\eta}_{\mu?\nu?} \Gamma^{\lambda?}_{\kappa?\eta} + \Gamma^{\eta}_{\mu?\kappa?} \Gamma^{\lambda?}_{\nu?\eta}; Gtog:= \Gamma^{\lambda?}_{\mu?\nu?} -> (1/2) * g^{\lambda?\kappa} ( \partial_{\nu?}{ g_{\kappa\mu?} } + \partial_{\mu?}{ g_{\kappa\nu?} } - \partial_{\kappa}{ g_{\mu?\nu?} } );
$$\displaystyle{}R^{\lambda?}\,_{\mu? \nu? \kappa?} \rightarrow -\partial_{\kappa?}{\Gamma^{\lambda?}\,_{\mu? \nu?}}% +\partial_{\nu?}{\Gamma^{\lambda?}\,_{\mu? \kappa?}}% -\Gamma^{\eta}\,_{\mu? \nu?} \Gamma^{\lambda?}\,_{\kappa? \eta}% +\Gamma^{\eta}\,_{\mu? \kappa?} \Gamma^{\lambda?}\,_{\nu? \eta}$$
$$\displaystyle{}\Gamma^{\lambda?}\,_{\mu? \nu?} \rightarrow \frac{1}{2}g^{\lambda? \kappa} \left(\partial_{\nu?}{g_{\kappa \mu?}}% +\partial_{\mu?}{g_{\kappa \nu?}}% -\partial_{\kappa}{g_{\mu? \nu?}}\right)$$
In this example we want to compute the Kaluza-Klein reduction of the $R_{m 4 n 4}$ component of the Riemann tensor. So we enter this component and do the substitution that takes the Riemann tensor to metrics. After each substitution, we distribute products over sums. We also apply the product rule to distribute derivatives over factors in a product.
todo:= g_{m1 m} R^{m1}_{4 n 4} + g_{4 m} R^{4}_{4 n 4};
$$\displaystyle{}g_{m1 m} R^{m1}\,_{4 n 4}% +g_{4 m} R^{4}\,_{4 n 4}$$
substitute(_, RtoG) substitute(_, Gtog) distribute(_) product_rule(_) distribute(_) sort_product(_);
$$\displaystyle{} - \frac{1}{2}\partial_{4}{g^{m1 \kappa}} \partial_{n}{g_{\kappa 4}} g_{m1 m}% - \frac{1}{2}\partial_{4 n}{g_{\kappa 4}} g_{m1 m} g^{m1 \kappa}% - \frac{1}{2}\partial_{4}{g_{\kappa n}} \partial_{4}{g^{m1 \kappa}} g_{m1 m}% - \frac{1}{2}\partial_{4 4}{g_{\kappa n}} g_{m1 m} g^{m1 \kappa}% +\frac{1}{2}\partial_{4}{g^{m1 \kappa}} \partial_{\kappa}{g_{4 n}} g_{m1 m}% +\frac{1}{2}\partial_{4 \kappa}{g_{4 n}} g_{m1 m} g^{m1 \kappa}% +\partial_{4}{g_{\kappa 4}} \partial_{n}{g^{m1 \kappa}} g_{m1 m}% +\partial_{n 4}{g_{\kappa 4}} g_{m1 m} g^{m1 \kappa}% - \frac{1}{2}\partial_{\kappa}{g_{4 4}} \partial_{n}{g^{m1 \kappa}} g_{m1 m}% - \frac{1}{2}\partial_{n \kappa}{g_{4 4}} g_{m1 m} g^{m1 \kappa}% - \frac{1}{4}\partial_{\eta}{g_{\mu 4}} \partial_{n}{g_{\kappa 4}} g^{\eta \kappa} g_{m1 m} g^{m1 \mu}% - \frac{1}{4}\partial_{4}{g_{\mu \eta}} \partial_{n}{g_{\kappa 4}} g^{\eta \kappa} g_{m1 m} g^{m1 \mu}% +\frac{1}{4}\partial_{\mu}{g_{4 \eta}} \partial_{n}{g_{\kappa 4}} g^{\eta \kappa} g_{m1 m} g^{m1 \mu}% - \frac{1}{4}\partial_{4}{g_{\kappa n}} \partial_{\eta}{g_{\mu 4}} g^{\eta \kappa} g_{m1 m} g^{m1 \mu}% - \frac{1}{4}\partial_{4}{g_{\mu \eta}} \partial_{4}{g_{\kappa n}} g^{\eta \kappa} g_{m1 m} g^{m1 \mu}% +\frac{1}{4}\partial_{4}{g_{\kappa n}} \partial_{\mu}{g_{4 \eta}} g^{\eta \kappa} g_{m1 m} g^{m1 \mu}% +\frac{1}{4}\partial_{\kappa}{g_{4 n}} \partial_{\eta}{g_{\mu 4}} g^{\eta \kappa} g_{m1 m} g^{m1 \mu}% +\frac{1}{4}\partial_{4}{g_{\mu \eta}} \partial_{\kappa}{g_{4 n}} g^{\eta \kappa} g_{m1 m} g^{m1 \mu}% - \frac{1}{4}\partial_{\mu}{g_{4 \eta}} \partial_{\kappa}{g_{4 n}} g^{\eta \kappa} g_{m1 m} g^{m1 \mu}% +\frac{1}{2}\partial_{4}{g_{\kappa 4}} \partial_{\eta}{g_{\mu n}} g^{\eta \kappa} g_{m1 m} g^{m1 \mu}% +\frac{1}{2}\partial_{4}{g_{\kappa 4}} \partial_{n}{g_{\mu \eta}} g^{\eta \kappa} g_{m1 m} g^{m1 \mu}% - \frac{1}{2}\partial_{4}{g_{\kappa 4}} \partial_{\mu}{g_{n \eta}} g^{\eta \kappa} g_{m1 m} g^{m1 \mu}% - \frac{1}{4}\partial_{\kappa}{g_{4 4}} \partial_{\eta}{g_{\mu n}} g^{\eta \kappa} g_{m1 m} g^{m1 \mu}% - \frac{1}{4}\partial_{\kappa}{g_{4 4}} \partial_{n}{g_{\mu \eta}} g^{\eta \kappa} g_{m1 m} g^{m1 \mu}% +\frac{1}{4}\partial_{\kappa}{g_{4 4}} \partial_{\mu}{g_{n \eta}} g^{\eta \kappa} g_{m1 m} g^{m1 \mu}% - \frac{1}{2}\partial_{4}{g^{4 \kappa}} \partial_{n}{g_{\kappa 4}} g_{4 m}% - \frac{1}{2}\partial_{4 n}{g_{\kappa 4}} g_{4 m} g^{4 \kappa}% - \frac{1}{2}\partial_{4}{g^{4 \kappa}} \partial_{4}{g_{\kappa n}} g_{4 m}% - \frac{1}{2}\partial_{4 4}{g_{\kappa n}} g_{4 m} g^{4 \kappa}% +\frac{1}{2}\partial_{4}{g^{4 \kappa}} \partial_{\kappa}{g_{4 n}} g_{4 m}% +\frac{1}{2}\partial_{4 \kappa}{g_{4 n}} g_{4 m} g^{4 \kappa}% +\partial_{4}{g_{\kappa 4}} \partial_{n}{g^{4 \kappa}} g_{4 m}% +\partial_{n 4}{g_{\kappa 4}} g_{4 m} g^{4 \kappa}% - \frac{1}{2}\partial_{\kappa}{g_{4 4}} \partial_{n}{g^{4 \kappa}} g_{4 m}% - \frac{1}{2}\partial_{n \kappa}{g_{4 4}} g_{4 m} g^{4 \kappa}% - \frac{1}{4}\partial_{\eta}{g_{\mu 4}} \partial_{n}{g_{\kappa 4}} g_{4 m} g^{4 \mu} g^{\eta \kappa}% - \frac{1}{4}\partial_{4}{g_{\mu \eta}} \partial_{n}{g_{\kappa 4}} g_{4 m} g^{4 \mu} g^{\eta \kappa}% +\frac{1}{4}\partial_{\mu}{g_{4 \eta}} \partial_{n}{g_{\kappa 4}} g_{4 m} g^{4 \mu} g^{\eta \kappa}% - \frac{1}{4}\partial_{4}{g_{\kappa n}} \partial_{\eta}{g_{\mu 4}} g_{4 m} g^{4 \mu} g^{\eta \kappa}% - \frac{1}{4}\partial_{4}{g_{\mu \eta}} \partial_{4}{g_{\kappa n}} g_{4 m} g^{4 \mu} g^{\eta \kappa}% +\frac{1}{4}\partial_{4}{g_{\kappa n}} \partial_{\mu}{g_{4 \eta}} g_{4 m} g^{4 \mu} g^{\eta \kappa}% +\frac{1}{4}\partial_{\kappa}{g_{4 n}} \partial_{\eta}{g_{\mu 4}} g_{4 m} g^{4 \mu} g^{\eta \kappa}% +\frac{1}{4}\partial_{4}{g_{\mu \eta}} \partial_{\kappa}{g_{4 n}} g_{4 m} g^{4 \mu} g^{\eta \kappa}% - \frac{1}{4}\partial_{\mu}{g_{4 \eta}} \partial_{\kappa}{g_{4 n}} g_{4 m} g^{4 \mu} g^{\eta \kappa}% +\frac{1}{2}\partial_{4}{g_{\kappa 4}} \partial_{\eta}{g_{\mu n}} g_{4 m} g^{4 \mu} g^{\eta \kappa}% +\frac{1}{2}\partial_{4}{g_{\kappa 4}} \partial_{n}{g_{\mu \eta}} g_{4 m} g^{4 \mu} g^{\eta \kappa}% - \frac{1}{2}\partial_{4}{g_{\kappa 4}} \partial_{\mu}{g_{n \eta}} g_{4 m} g^{4 \mu} g^{\eta \kappa}% - \frac{1}{4}\partial_{\kappa}{g_{4 4}} \partial_{\eta}{g_{\mu n}} g_{4 m} g^{4 \mu} g^{\eta \kappa}% - \frac{1}{4}\partial_{\kappa}{g_{4 4}} \partial_{n}{g_{\mu \eta}} g_{4 m} g^{4 \mu} g^{\eta \kappa}% +\frac{1}{4}\partial_{\kappa}{g_{4 4}} \partial_{\mu}{g_{n \eta}} g_{4 m} g^{4 \mu} g^{\eta \kappa}$$
None of the tensors will depend on the 4-th direction, so we set derivatives in that direction to zero.
substitute(_, $\partial_{4}{A??} -> 0$, repeat=True) substitute(_, $\partial_{4 m?}{A??} -> 0$, repeat=True) substitute(_, $\partial_{m? 4}{A??} -> 0$, repeat=True);
$$\displaystyle{} - \frac{1}{2}\partial_{\kappa}{g_{4 4}} \partial_{n}{g^{m1 \kappa}} g_{m1 m}% - \frac{1}{2}\partial_{n \kappa}{g_{4 4}} g_{m1 m} g^{m1 \kappa}% - \frac{1}{4}\partial_{\eta}{g_{\mu 4}} \partial_{n}{g_{\kappa 4}} g^{\eta \kappa} g_{m1 m} g^{m1 \mu}% +\frac{1}{4}\partial_{\mu}{g_{4 \eta}} \partial_{n}{g_{\kappa 4}} g^{\eta \kappa} g_{m1 m} g^{m1 \mu}% +\frac{1}{4}\partial_{\kappa}{g_{4 n}} \partial_{\eta}{g_{\mu 4}} g^{\eta \kappa} g_{m1 m} g^{m1 \mu}% - \frac{1}{4}\partial_{\mu}{g_{4 \eta}} \partial_{\kappa}{g_{4 n}} g^{\eta \kappa} g_{m1 m} g^{m1 \mu}% - \frac{1}{4}\partial_{\kappa}{g_{4 4}} \partial_{\eta}{g_{\mu n}} g^{\eta \kappa} g_{m1 m} g^{m1 \mu}% - \frac{1}{4}\partial_{\kappa}{g_{4 4}} \partial_{n}{g_{\mu \eta}} g^{\eta \kappa} g_{m1 m} g^{m1 \mu}% +\frac{1}{4}\partial_{\kappa}{g_{4 4}} \partial_{\mu}{g_{n \eta}} g^{\eta \kappa} g_{m1 m} g^{m1 \mu}% - \frac{1}{2}\partial_{\kappa}{g_{4 4}} \partial_{n}{g^{4 \kappa}} g_{4 m}% - \frac{1}{2}\partial_{n \kappa}{g_{4 4}} g_{4 m} g^{4 \kappa}% - \frac{1}{4}\partial_{\eta}{g_{\mu 4}} \partial_{n}{g_{\kappa 4}} g_{4 m} g^{4 \mu} g^{\eta \kappa}% +\frac{1}{4}\partial_{\mu}{g_{4 \eta}} \partial_{n}{g_{\kappa 4}} g_{4 m} g^{4 \mu} g^{\eta \kappa}% +\frac{1}{4}\partial_{\kappa}{g_{4 n}} \partial_{\eta}{g_{\mu 4}} g_{4 m} g^{4 \mu} g^{\eta \kappa}% - \frac{1}{4}\partial_{\mu}{g_{4 \eta}} \partial_{\kappa}{g_{4 n}} g_{4 m} g^{4 \mu} g^{\eta \kappa}% - \frac{1}{4}\partial_{\kappa}{g_{4 4}} \partial_{\eta}{g_{\mu n}} g_{4 m} g^{4 \mu} g^{\eta \kappa}% - \frac{1}{4}\partial_{\kappa}{g_{4 4}} \partial_{n}{g_{\mu \eta}} g_{4 m} g^{4 \mu} g^{\eta \kappa}% +\frac{1}{4}\partial_{\kappa}{g_{4 4}} \partial_{\mu}{g_{n \eta}} g_{4 m} g^{4 \mu} g^{\eta \kappa}$$
We then split the greek indices into roman indices and the 4-th direction. Notice the use of repeat=True which will keep doing this substitution until all greek indices have been split. After that we again have to set all derivatives in the 4-th direction to zero.
split_index(_, $\mu, m1, 4$, repeat=True) substitute(_, $\partial_{4}{A??} -> 0$, repeat=True) substitute(_, $\partial_{4 m?}{A??} -> 0$, repeat=True) substitute(_, $\partial_{m? 4}{A??} -> 0$, repeat=True) canonicalise(_);
$$\displaystyle{} - \frac{1}{2}\partial_{p}{g_{4 4}} \partial_{n}{g^{m1 p}} g_{m1 m}% - \frac{1}{2}\partial_{n p}{g_{4 4}} g_{m1 m} g^{m1 p}% - \frac{1}{4}\partial_{p}{g_{r 4}} \partial_{n}{g_{q 4}} g^{p q} g_{m1 m} g^{m1 r}% - \frac{1}{4}\partial_{p}{g_{4 4}} \partial_{n}{g_{q 4}} g^{p q} g_{m1 m} g^{m1 4}% - \frac{1}{4}\partial_{p}{g_{q 4}} \partial_{n}{g_{4 4}} g^{p 4} g_{m1 m} g^{m1 q}% - \frac{1}{4}\partial_{p}{g_{4 4}} \partial_{n}{g_{4 4}} g^{p 4} g_{m1 m} g^{m1 4}% +\frac{1}{4}\partial_{r}{g_{4 p}} \partial_{n}{g_{q 4}} g^{p q} g_{m1 m} g^{m1 r}% +\frac{1}{4}\partial_{q}{g_{4 p}} \partial_{n}{g_{4 4}} g^{p 4} g_{m1 m} g^{m1 q}% +\frac{1}{4}\partial_{q}{g_{4 4}} \partial_{n}{g_{p 4}} g^{4 p} g_{m1 m} g^{m1 q}% +\frac{1}{4}\partial_{p}{g_{4 4}} \partial_{n}{g_{4 4}} g^{4 4} g_{m1 m} g^{m1 p}% +\frac{1}{4}\partial_{q}{g_{4 n}} \partial_{p}{g_{r 4}} g^{p q} g_{m1 m} g^{m1 r}% +\frac{1}{4}\partial_{q}{g_{4 n}} \partial_{p}{g_{4 4}} g^{p q} g_{m1 m} g^{m1 4}% - \frac{1}{4}\partial_{r}{g_{4 p}} \partial_{q}{g_{4 n}} g^{p q} g_{m1 m} g^{m1 r}% - \frac{1}{4}\partial_{q}{g_{4 4}} \partial_{p}{g_{4 n}} g^{4 p} g_{m1 m} g^{m1 q}% - \frac{1}{4}\partial_{q}{g_{4 4}} \partial_{p}{g_{r n}} g^{p q} g_{m1 m} g^{m1 r}% - \frac{1}{4}\partial_{q}{g_{4 4}} \partial_{p}{g_{4 n}} g^{p q} g_{m1 m} g^{m1 4}% - \frac{1}{4}\partial_{q}{g_{4 4}} \partial_{n}{g_{r p}} g^{p q} g_{m1 m} g^{m1 r}% - \frac{1}{4}\partial_{q}{g_{4 4}} \partial_{n}{g_{4 p}} g^{p q} g_{m1 m} g^{m1 4}% - \frac{1}{4}\partial_{p}{g_{4 4}} \partial_{n}{g_{q 4}} g^{4 p} g_{m1 m} g^{m1 q}% - \frac{1}{4}\partial_{p}{g_{4 4}} \partial_{n}{g_{4 4}} g^{4 p} g_{m1 m} g^{m1 4}% +\frac{1}{4}\partial_{q}{g_{4 4}} \partial_{r}{g_{n p}} g^{p q} g_{m1 m} g^{m1 r}% +\frac{1}{4}\partial_{p}{g_{4 4}} \partial_{q}{g_{n 4}} g^{4 p} g_{m1 m} g^{m1 q}% - \frac{1}{2}\partial_{p}{g_{4 4}} \partial_{n}{g^{4 p}} g_{4 m}% - \frac{1}{2}\partial_{n p}{g_{4 4}} g_{4 m} g^{4 p}% - \frac{1}{4}\partial_{p}{g_{r 4}} \partial_{n}{g_{q 4}} g_{4 m} g^{4 r} g^{p q}% - \frac{1}{4}\partial_{p}{g_{4 4}} \partial_{n}{g_{q 4}} g_{4 m} g^{4 4} g^{p q}% - \frac{1}{4}\partial_{p}{g_{q 4}} \partial_{n}{g_{4 4}} g_{4 m} g^{4 q} g^{p 4}% - \frac{1}{4}\partial_{p}{g_{4 4}} \partial_{n}{g_{4 4}} g_{4 m} g^{4 4} g^{p 4}% +\frac{1}{4}\partial_{r}{g_{4 p}} \partial_{n}{g_{q 4}} g_{4 m} g^{4 r} g^{p q}% +\frac{1}{4}\partial_{q}{g_{4 p}} \partial_{n}{g_{4 4}} g_{4 m} g^{4 q} g^{p 4}% +\frac{1}{4}\partial_{q}{g_{4 4}} \partial_{n}{g_{p 4}} g_{4 m} g^{4 q} g^{4 p}% +\frac{1}{4}\partial_{p}{g_{4 4}} \partial_{n}{g_{4 4}} g_{4 m} g^{4 p} g^{4 4}% +\frac{1}{4}\partial_{q}{g_{4 n}} \partial_{p}{g_{r 4}} g_{4 m} g^{4 r} g^{p q}% +\frac{1}{4}\partial_{q}{g_{4 n}} \partial_{p}{g_{4 4}} g_{4 m} g^{4 4} g^{p q}% - \frac{1}{4}\partial_{r}{g_{4 p}} \partial_{q}{g_{4 n}} g_{4 m} g^{4 r} g^{p q}% - \frac{1}{4}\partial_{q}{g_{4 4}} \partial_{p}{g_{4 n}} g_{4 m} g^{4 q} g^{4 p}% - \frac{1}{4}\partial_{q}{g_{4 4}} \partial_{p}{g_{r n}} g_{4 m} g^{4 r} g^{p q}% - \frac{1}{4}\partial_{q}{g_{4 4}} \partial_{p}{g_{4 n}} g_{4 m} g^{4 4} g^{p q}% - \frac{1}{4}\partial_{q}{g_{4 4}} \partial_{n}{g_{r p}} g_{4 m} g^{4 r} g^{p q}% - \frac{1}{4}\partial_{q}{g_{4 4}} \partial_{n}{g_{4 p}} g_{4 m} g^{4 4} g^{p q}% - \frac{1}{4}\partial_{p}{g_{4 4}} \partial_{n}{g_{q 4}} g_{4 m} g^{4 q} g^{4 p}% - \frac{1}{4}\partial_{p}{g_{4 4}} \partial_{n}{g_{4 4}} g_{4 m} g^{4 4} g^{4 p}% +\frac{1}{4}\partial_{q}{g_{4 4}} \partial_{r}{g_{n p}} g_{4 m} g^{4 r} g^{p q}% +\frac{1}{4}\partial_{p}{g_{4 4}} \partial_{q}{g_{n 4}} g_{4 m} g^{4 q} g^{4 p}$$
Now comes the ansatz of the metric in terms of the Kaluza-Klein gauge field $A_\mu$ and the scalar $\phi$. We substitute this with the lines below. The output is not particularly enlightening...
substitute(_, $g_{4 4} -> \phi$ ) substitute(_, $g_{m 4} -> \phi A_{m}$ ) substitute(_, $g_{4 m} -> \phi A_{m}$ ) substitute(_, $g_{m n} -> \phi**{-1} h_{m n} + \phi A_{m} A_{n}$ ) substitute(_, $g^{4 4} -> \phi**{-1} + \phi A_{m} h^{m n} A_{n}$ ) substitute(_, $g^{m 4} -> - \phi h^{m n} A_{n}$ ) substitute(_, $g^{4 m} -> - \phi h^{m n} A_{n}$ ) substitute(_, $g^{m n} -> \phi h^{m n}$ );
$$\displaystyle{} - \frac{1}{2}\partial_{p}{\phi} \partial_{n}\left(\phi h^{m1 p}\right) \left(\phi^{-1} h_{m1 m}% +\phi A_{m1} A_{m}\right)% - \frac{1}{2}\partial_{n p}{\phi} \left(\phi^{-1} h_{m1 m}% +\phi A_{m1} A_{m}\right) \phi h^{m1 p}% - \frac{1}{4}\partial_{p}\left(\phi A_{r}\right) \partial_{n}\left(\phi A_{q}\right) \phi h^{p q} \left(\phi^{-1} h_{m1 m}% +\phi A_{m1} A_{m}\right) \phi h^{m1 r}% +\frac{1}{4}\partial_{p}{\phi} \partial_{n}\left(\phi A_{q}\right) \phi h^{p q} \left(\phi^{-1} h_{m1 m}% +\phi A_{m1} A_{m}\right) \phi h^{m1 r} A_{r}% +\frac{1}{4}\partial_{p}\left(\phi A_{q}\right) \partial_{n}{\phi} \phi h^{p r} A_{r} \left(\phi^{-1} h_{m1 m}% +\phi A_{m1} A_{m}\right) \phi h^{m1 q}% - \frac{1}{4}\partial_{p}{\phi} \partial_{n}{\phi} \phi h^{p q} A_{q} \left(\phi^{-1} h_{m1 m}% +\phi A_{m1} A_{m}\right) \phi h^{m1 r} A_{r}% +\frac{1}{4}\partial_{r}\left(\phi A_{p}\right) \partial_{n}\left(\phi A_{q}\right) \phi h^{p q} \left(\phi^{-1} h_{m1 m}% +\phi A_{m1} A_{m}\right) \phi h^{m1 r}% - \frac{1}{4}\partial_{q}\left(\phi A_{p}\right) \partial_{n}{\phi} \phi h^{p r} A_{r} \left(\phi^{-1} h_{m1 m}% +\phi A_{m1} A_{m}\right) \phi h^{m1 q}% - \frac{1}{4}\partial_{q}{\phi} \partial_{n}\left(\phi A_{p}\right) \phi h^{p r} A_{r} \left(\phi^{-1} h_{m1 m}% +\phi A_{m1} A_{m}\right) \phi h^{m1 q}% +\frac{1}{4}\partial_{p}{\phi} \partial_{n}{\phi} \left(\phi^{-1}% +\phi A_{q} h^{q r} A_{r}\right) \left(\phi^{-1} h_{m1 m}% +\phi A_{m1} A_{m}\right) \phi h^{m1 p}% +\frac{1}{4}\partial_{q}\left(\phi A_{n}\right) \partial_{p}\left(\phi A_{r}\right) \phi h^{p q} \left(\phi^{-1} h_{m1 m}% +\phi A_{m1} A_{m}\right) \phi h^{m1 r}% - \frac{1}{4}\partial_{q}\left(\phi A_{n}\right) \partial_{p}{\phi} \phi h^{p q} \left(\phi^{-1} h_{m1 m}% +\phi A_{m1} A_{m}\right) \phi h^{m1 r} A_{r}% - \frac{1}{4}\partial_{r}\left(\phi A_{p}\right) \partial_{q}\left(\phi A_{n}\right) \phi h^{p q} \left(\phi^{-1} h_{m1 m}% +\phi A_{m1} A_{m}\right) \phi h^{m1 r}% +\frac{1}{4}\partial_{q}{\phi} \partial_{p}\left(\phi A_{n}\right) \phi h^{p r} A_{r} \left(\phi^{-1} h_{m1 m}% +\phi A_{m1} A_{m}\right) \phi h^{m1 q}% - \frac{1}{4}\partial_{q}{\phi} \partial_{p}\left(\phi^{-1} h_{r n}% +\phi A_{r} A_{n}\right) \phi h^{p q} \left(\phi^{-1} h_{m1 m}% +\phi A_{m1} A_{m}\right) \phi h^{m1 r}% +\frac{1}{4}\partial_{q}{\phi} \partial_{p}\left(\phi A_{n}\right) \phi h^{p q} \left(\phi^{-1} h_{m1 m}% +\phi A_{m1} A_{m}\right) \phi h^{m1 r} A_{r}% - \frac{1}{4}\partial_{q}{\phi} \partial_{n}\left(\phi^{-1} h_{r p}% +\phi A_{r} A_{p}\right) \phi h^{p q} \left(\phi^{-1} h_{m1 m}% +\phi A_{m1} A_{m}\right) \phi h^{m1 r}% +\frac{1}{4}\partial_{q}{\phi} \partial_{n}\left(\phi A_{p}\right) \phi h^{p q} \left(\phi^{-1} h_{m1 m}% +\phi A_{m1} A_{m}\right) \phi h^{m1 r} A_{r}% +\frac{1}{4}\partial_{p}{\phi} \partial_{n}\left(\phi A_{q}\right) \phi h^{p r} A_{r} \left(\phi^{-1} h_{m1 m}% +\phi A_{m1} A_{m}\right) \phi h^{m1 q}% - \frac{1}{4}\partial_{p}{\phi} \partial_{n}{\phi} \phi h^{p r} A_{r} \left(\phi^{-1} h_{m1 m}% +\phi A_{m1} A_{m}\right) \phi h^{m1 q} A_{q}% +\frac{1}{4}\partial_{q}{\phi} \partial_{r}\left(\phi^{-1} h_{n p}% +\phi A_{n} A_{p}\right) \phi h^{p q} \left(\phi^{-1} h_{m1 m}% +\phi A_{m1} A_{m}\right) \phi h^{m1 r}% - \frac{1}{4}\partial_{p}{\phi} \partial_{q}\left(\phi A_{n}\right) \phi h^{p r} A_{r} \left(\phi^{-1} h_{m1 m}% +\phi A_{m1} A_{m}\right) \phi h^{m1 q}% +\frac{1}{2}\partial_{p}{\phi} \partial_{n}\left(\phi h^{p q} A_{q}\right) \phi A_{m}% +\frac{1}{2}\partial_{n p}{\phi} \phi A_{m} \phi h^{p q} A_{q}% +\frac{1}{4}\partial_{p}\left(\phi A_{r}\right) \partial_{n}\left(\phi A_{q}\right) \phi A_{m} \phi h^{r s} A_{s} \phi h^{p q}% - \frac{1}{4}\partial_{p}{\phi} \partial_{n}\left(\phi A_{q}\right) \phi A_{m} \left(\phi^{-1}% +\phi A_{r} h^{r s} A_{s}\right) \phi h^{p q}% - \frac{1}{4}\partial_{p}\left(\phi A_{q}\right) \partial_{n}{\phi} \phi A_{m} \phi h^{q s} A_{s} \phi h^{p r} A_{r}% +\frac{1}{2}\partial_{p}{\phi} \partial_{n}{\phi} \phi A_{m} \left(\phi^{-1}% +\phi A_{q} h^{q r} A_{r}\right) \phi h^{p s} A_{s}% - \frac{1}{4}\partial_{r}\left(\phi A_{p}\right) \partial_{n}\left(\phi A_{q}\right) \phi A_{m} \phi h^{r s} A_{s} \phi h^{p q}% +\frac{1}{4}\partial_{q}\left(\phi A_{p}\right) \partial_{n}{\phi} \phi A_{m} \phi h^{q s} A_{s} \phi h^{p r} A_{r}% +\frac{1}{4}\partial_{q}{\phi} \partial_{n}\left(\phi A_{p}\right) \phi A_{m} \phi h^{q r} A_{r} \phi h^{p s} A_{s}% - \frac{1}{4}\partial_{p}{\phi} \partial_{n}{\phi} \phi A_{m} \phi h^{p s} A_{s} \left(\phi^{-1}% +\phi A_{q} h^{q r} A_{r}\right)% - \frac{1}{4}\partial_{q}\left(\phi A_{n}\right) \partial_{p}\left(\phi A_{r}\right) \phi A_{m} \phi h^{r s} A_{s} \phi h^{p q}% +\frac{1}{4}\partial_{q}\left(\phi A_{n}\right) \partial_{p}{\phi} \phi A_{m} \left(\phi^{-1}% +\phi A_{r} h^{r s} A_{s}\right) \phi h^{p q}% +\frac{1}{4}\partial_{r}\left(\phi A_{p}\right) \partial_{q}\left(\phi A_{n}\right) \phi A_{m} \phi h^{r s} A_{s} \phi h^{p q}% - \frac{1}{4}\partial_{q}{\phi} \partial_{p}\left(\phi A_{n}\right) \phi A_{m} \phi h^{q r} A_{r} \phi h^{p s} A_{s}% +\frac{1}{4}\partial_{q}{\phi} \partial_{p}\left(\phi^{-1} h_{r n}% +\phi A_{r} A_{n}\right) \phi A_{m} \phi h^{r s} A_{s} \phi h^{p q}% - \frac{1}{4}\partial_{q}{\phi} \partial_{p}\left(\phi A_{n}\right) \phi A_{m} \left(\phi^{-1}% +\phi A_{r} h^{r s} A_{s}\right) \phi h^{p q}% +\frac{1}{4}\partial_{q}{\phi} \partial_{n}\left(\phi^{-1} h_{r p}% +\phi A_{r} A_{p}\right) \phi A_{m} \phi h^{r s} A_{s} \phi h^{p q}% - \frac{1}{4}\partial_{q}{\phi} \partial_{n}\left(\phi A_{p}\right) \phi A_{m} \left(\phi^{-1}% +\phi A_{r} h^{r s} A_{s}\right) \phi h^{p q}% - \frac{1}{4}\partial_{p}{\phi} \partial_{n}\left(\phi A_{q}\right) \phi A_{m} \phi h^{q r} A_{r} \phi h^{p s} A_{s}% - \frac{1}{4}\partial_{q}{\phi} \partial_{r}\left(\phi^{-1} h_{n p}% +\phi A_{n} A_{p}\right) \phi A_{m} \phi h^{r s} A_{s} \phi h^{p q}% +\frac{1}{4}\partial_{p}{\phi} \partial_{q}\left(\phi A_{n}\right) \phi A_{m} \phi h^{q r} A_{r} \phi h^{p s} A_{s}$$
We now need to distribute derivatives over sums and use the product rule (twice). The output of those operations is quite large, so we suppress it (by not writing semi-colons after the function calls).
distribute(_) product_rule(_) distribute(_) product_rule(_) canonicalise(_)
substitute(_, $\partial_{p}{h^{n m}} h_{q m} -> - \partial_{p}{h_{q m}} h^{n m}$ ) collect_factors(_) sort_product(_) converge(todo): substitute(_, $h_{m1 m2} h^{m3 m2} -> \delta_{m1}^{m3}$, repeat=True ) eliminate_kronecker(_) canonicalise(_) ;
$$\displaystyle{} - \frac{1}{4}\partial_{m}{\phi} \partial_{n}{\phi} \phi^{-1}% +\frac{1}{4}\partial_{m1}{\phi} \partial_{n}{h_{m p}} h^{m1 p}% - \frac{3}{4}A_{m} A_{m1} \partial_{n}{\phi} \partial_{p}{\phi} \phi h^{m1 p}% - \frac{1}{2}A_{m} A_{m1} \partial_{p}{\phi} \partial_{n}{h^{m1 p}} \phi^{2}% - \frac{1}{2}\partial_{m n}{\phi}% - \frac{1}{2}A_{m} A_{m1} \partial_{n p}{\phi} \phi^{2} h^{m1 p}% - \frac{1}{2}A_{m} \partial_{n}{A_{m1}} \partial_{q}{\phi} \phi^{2} h^{m1 q}% - \frac{1}{4}A_{m1} \partial_{q}{A_{m}} \partial_{n}{\phi} \phi^{2} h^{m1 q}% - \frac{1}{4}\partial_{n}{A_{m1}} \partial_{p}{A_{m}} \phi^{3} h^{m1 p}% - \frac{1}{4}A_{m} A_{m1} \partial_{n}{A_{p}} \partial_{q}{A_{r}} \phi^{5} h^{m1 r} h^{p q}% +\frac{3}{4}A_{m} A_{m1} \partial_{n}{\phi} \partial_{q}{\phi} \phi h^{m1 q}% +\frac{1}{2}A_{m} \partial_{n}{A_{p}} \partial_{q}{\phi} \phi^{2} h^{p q}% +\frac{1}{4}A_{m1} \partial_{r}{A_{m}} \partial_{n}{\phi} \phi^{2} h^{m1 r}% +\frac{1}{4}A_{m1} A_{q} \partial_{m}{\phi} \partial_{n}{\phi} \phi h^{m1 q}% +\frac{1}{4}A_{p} \partial_{n}{A_{q}} \partial_{m}{\phi} \phi^{2} h^{p q}% +\frac{1}{4}A_{p} \partial_{m}{A_{q}} \partial_{n}{\phi} \phi^{2} h^{p q}% +\frac{1}{4}\partial_{m}{A_{m1}} \partial_{n}{A_{q}} \phi^{3} h^{m1 q}% +\frac{1}{4}A_{m} A_{m1} \partial_{n}{A_{p}} \partial_{q}{A_{r}} \phi^{5} h^{m1 q} h^{p r}% - \frac{1}{4}A_{m1} A_{r} \partial_{m}{\phi} \partial_{n}{\phi} \phi h^{m1 r}% - \frac{1}{4}A_{p} \partial_{m}{A_{r}} \partial_{n}{\phi} \phi^{2} h^{p r}% - \frac{1}{4}A_{p} \partial_{n}{A_{r}} \partial_{m}{\phi} \phi^{2} h^{p r}% +\frac{1}{4}A_{m} \partial_{m1}{A_{n}} \partial_{p}{\phi} \phi^{2} h^{m1 p}% +\frac{1}{4}\partial_{m1}{A_{m}} \partial_{p}{A_{n}} \phi^{3} h^{m1 p}% +\frac{1}{4}A_{m} A_{m1} \partial_{p}{A_{n}} \partial_{q}{A_{r}} \phi^{5} h^{m1 r} h^{p q}% - \frac{1}{4}A_{p} \partial_{q}{A_{n}} \partial_{m}{\phi} \phi^{2} h^{p q}% - \frac{1}{4}A_{n} \partial_{m}{A_{p}} \partial_{q}{\phi} \phi^{2} h^{p q}% - \frac{1}{4}\partial_{m}{A_{m1}} \partial_{q}{A_{n}} \phi^{3} h^{m1 q}% - \frac{1}{4}A_{m} A_{m1} \partial_{p}{A_{n}} \partial_{q}{A_{r}} \phi^{5} h^{m1 q} h^{p r}% +\frac{1}{4}A_{p} \partial_{r}{A_{n}} \partial_{m}{\phi} \phi^{2} h^{p r}% +\frac{1}{4}\partial_{m1}{\phi} \partial_{p}{\phi} \phi^{-1} h_{m n} h^{m1 p}% - \frac{1}{4}\partial_{m1}{\phi} \partial_{p}{h_{m n}} h^{m1 p}% +\frac{1}{4}A_{m} A_{n} \partial_{p}{\phi} \partial_{q}{\phi} \phi h^{p q}% - \frac{1}{4}A_{m} A_{m1} \partial_{p}{\phi} \partial_{q}{h_{n r}} \phi^{2} h^{m1 r} h^{p q}% - \frac{1}{4}A_{m} \partial_{m1}{A_{n}} \partial_{q}{\phi} \phi^{2} h^{m1 q}% - \frac{1}{4}A_{m} A_{m1} \partial_{p}{\phi} \partial_{n}{h_{q r}} \phi^{2} h^{m1 q} h^{p r}% - \frac{1}{4}A_{m1} \partial_{n}{A_{m}} \partial_{q}{\phi} \phi^{2} h^{m1 q}% +\frac{1}{4}A_{m1} \partial_{n}{A_{m}} \partial_{r}{\phi} \phi^{2} h^{m1 r}% +\frac{1}{4}\partial_{m1}{\phi} \partial_{m}{h_{n q}} h^{m1 q}% - \frac{1}{4}A_{m} A_{p} \partial_{n}{\phi} \partial_{r}{\phi} \phi h^{p r}% +\frac{1}{4}A_{m} A_{m1} \partial_{p}{\phi} \partial_{q}{h_{n r}} \phi^{2} h^{m1 q} h^{p r}% +\frac{1}{4}A_{m1} \partial_{m}{A_{n}} \partial_{q}{\phi} \phi^{2} h^{m1 q}% +\frac{1}{4}A_{n} \partial_{m}{A_{m1}} \partial_{q}{\phi} \phi^{2} h^{m1 q}% - \frac{1}{4}A_{m1} \partial_{m}{A_{n}} \partial_{r}{\phi} \phi^{2} h^{m1 r}% +\frac{1}{2}A_{m} A_{p} \partial_{q}{\phi} \partial_{n}{h^{p q}} \phi^{2}% +\frac{1}{2}A_{m} A_{p} \partial_{n q}{\phi} \phi^{2} h^{p q}% +\frac{1}{4}A_{m} A_{p} \partial_{n}{A_{q}} \partial_{r}{A_{s}} \phi^{5} h^{p s} h^{q r}% +\frac{1}{4}A_{m} A_{p} \partial_{n}{\phi} \partial_{s}{\phi} \phi h^{p s}% - \frac{1}{4}A_{m} A_{p} \partial_{n}{A_{q}} \partial_{r}{A_{s}} \phi^{5} h^{p r} h^{q s}% - \frac{1}{4}A_{m} A_{p} \partial_{q}{A_{n}} \partial_{r}{A_{s}} \phi^{5} h^{p s} h^{q r}% +\frac{1}{4}A_{m} A_{p} \partial_{q}{A_{n}} \partial_{r}{A_{s}} \phi^{5} h^{p r} h^{q s}% - \frac{1}{4}A_{m} A_{n} \partial_{q}{\phi} \partial_{r}{\phi} \phi h^{q r}% +\frac{1}{4}A_{m} A_{p} \partial_{q}{\phi} \partial_{r}{h_{n s}} \phi^{2} h^{p s} h^{q r}% - \frac{1}{4}A_{m} A_{p} \partial_{n}{\phi} \partial_{q}{\phi} \phi h^{p q}% +\frac{1}{4}A_{m} A_{p} \partial_{q}{\phi} \partial_{n}{h_{r s}} \phi^{2} h^{p r} h^{q s}% +\frac{1}{4}A_{m} A_{q} \partial_{n}{\phi} \partial_{s}{\phi} \phi h^{q s}% - \frac{1}{4}A_{m} A_{p} \partial_{q}{\phi} \partial_{r}{h_{n s}} \phi^{2} h^{p r} h^{q s}$$
We can now finally rewrite the derivative of the gauge potential in terms of the field strength to obtain the result.
substitute(_, $\partial_{n}{A_{m}} -> 1/2*\partial_{n}{A_{m}} + 1/2*F_{n m} + 1/2*\partial_{m}{A_{n}}$ ) distribute(_) sort_product(_) canonicalise(_) rename_dummies(_);
$$\displaystyle{} - \frac{1}{4}\partial_{m}{\phi} \partial_{n}{\phi} \phi^{-1}% +\frac{1}{4}\partial_{p}{\phi} \partial_{n}{h_{m q}} h^{p q}% - \frac{1}{2}\partial_{m n}{\phi}% +\frac{1}{4}F_{m p} F_{n q} \phi^{3} h^{p q}% +\frac{1}{4}\partial_{p}{\phi} \partial_{q}{\phi} \phi^{-1} h_{m n} h^{p q}% - \frac{1}{4}\partial_{p}{\phi} \partial_{q}{h_{m n}} h^{p q}% +\frac{1}{4}\partial_{p}{\phi} \partial_{m}{h_{n q}} h^{p q}$$