Partial derviative indices issues

Question

As a new user of Cadabra I am a bit confused what is going on with the following code. I am trying to repeat the steps done in Schutz's GR book for equations 6.63 - 6.67 as shown below.

Note I have edited my original post and added print outputs for when I substitute Gammpd into the equation for R. Cut and paste of my outputs produced garbage.

{a,b,c,d,e,f,g,h,i,j,k,l,m,n,o,p,q,r,s,t,u#}::Indices(name=fixed,position=fixed).
{\alpha,\beta,\nu,\mu,\sigma}::Indices(name=fixed,position=fixed);
{\rho,\delta,\kappa,\theta}::Indices(name=fixed,position=fixed);

\partial{#}::PartialDerivative.

\Gamma^{a}_{b c}::TableauSymmetry(shape={2}, indices={1,2});

g_{a b}::Metric.



;::Symbol;
            # see https://cadabra.science/qa/473/is-this-legal-syntax
            # this code works with and without this trick

R := R^{\alpha}_{\beta \mu \nu} -> \partial_{\mu}{\Gamma^{\alpha}_{\beta \nu}}  
                                -\partial_{\nu}{\Gamma^{\alpha}_{\beta \mu}}  
                                +\Gamma^{\alpha}_{\sigma \mu} \Gamma^{\sigma}_{\beta \nu}  
                                -\Gamma^{\alpha}_{\sigma \nu} \Gamma^{\sigma}_{\beta \mu};

 exp := R^{\alpha}_{\beta \mu \nu};

 substitute(exp,R);

 Gammapd := \partial_{\sigma}{\Gamma^{\alpha}_{\mu \nu}} -> (1/2) g^{\alpha \beta}  
            (\partial_{\nu \sigma}{g_{\beta \mu}} + \partial_{\mu \sigma}{g_{\beta \nu}}  
            -\partial_{\beta \sigma}{g_{\mu \nu}});

 substitute(exp,Gammapd);

Doing the above results in :

1/2 g^{α β}(\partial{ν μ}(g{β β})) - 1/2 g^{α β}(\partial{μ ν}(g{β β})) + Γ^{α}{σ μ} Γ^{σ}{β ν}-Γ^{α}{σ ν} Γ^{σ}{β μ}

Using the following for Gammapd I get what I think is correct.

Gammapd := \partial_{\sigma}{\Gamma^{\alpha}_{\mu \nu}} -> (1/2) g^{\alpha a}  
            (\partial_{\nu \sigma}{g_{a \mu}} + \partial_{\mu \sigma}{g_{a \nu}}  
            -\partial_{a \sigma}{g_{\mu \nu}});

1/2 g^{α a}(\partial{ν μ}(g{a β}) + \partial{β μ}(g{a ν})-\partial{a μ}(g{β ν})) - 1/2 g^{α a}(\partial{μ ν}(g{a β}) + \partial{β ν}(g{a μ})-\partial{a ν}(g{β μ})) + Γ^{α}{σ μ} Γ^{σ}{β ν}-Γ^{α}{σ ν} Γ^{σ}{β μ}

Note here that that \partial{ν μ}(g{a β}) should cancel out with \partial{μ ν}(g{a β}. In Shutz we are considering the inertial frame so all the Gammas = 0.

From Schutz: (the alpha on R is contravariant as are the metric's alpha sigma indices)

Rα βμν = 12 gασ (gσβ,νμ + gσ ν,βμ − gβν,σ μ − gσβ,μν − gσ μ,βν + gβμ,σ ν ). (6.65)

I would paste pictures of the Schutz equation and of my output but I can't seem to do that here and I'm not sure why?

I have also used:

Gammapd := \partial_{\sigma}{\Gamma^{\alpha}_{\mu \nu}} -> (1/2) g^{a b}  
            (\partial_{c d}{g_{b e}} + \partial_{e d}{g_{b c}}  
            -\partial_{b d}{g_{e c}});

for my substitution and for that I get: Γ^{α}{σ μ} Γ^{σ}{β ν}-Γ^{α}{σ ν} Γ^{σ}{β μ}

only the Gamma terms?

Do I somehow have to use Depends or something to get this to work. Why are the terms with the double partial derivatives not canceling when the indices are symmetric?