dbeta[la, uj]*dbeta[lj, ub] - em4phi*dalpha[la]*dalpha[lj]*gtu[ub, uj] + 
 (alpha*em4phi*dalpha[lj]*dphi[la]*gtu[ub, uj])/3 - 
 (alpha*em4phi*dalpha[la]*dphi[lj]*gtu[ub, uj])/3 + 
 (alpha^2*em4phi*dphi[la]*dphi[lj]*gtu[ub, uj])/18 - 
 alpha*em4phi*dalpha[la]*dgt[lj, li, lk]*gtu[ub, uj]*gtu[ui, uk] + 
 (alpha^2*em4phi*dgt[lj, li, lk]*dphi[la]*gtu[ub, uj]*gtu[ui, uk])/6 + 
 alpha*em4phi*dalpha[lj]*dgt[la, li, lk]*gtu[ub, ui]*gtu[uj, uk] + 
 2*alpha*em4phi*dalpha[la]*dgt[lj, li, lk]*gtu[ub, ui]*gtu[uj, uk] - 
 (alpha^2*em4phi*dgt[lj, li, lk]*dphi[la]*gtu[ub, ui]*gtu[uj, uk])/3 + 
 (alpha^2*em4phi*dgt[la, li, lk]*dphi[lj]*gtu[ub, ui]*gtu[uj, uk])/6 - 
 alpha^2*em4phi*dgt[la, lj, li]*dgt[lk, ll, lm]*gtu[ub, ul]*gtu[ui, um]*
  gtu[uj, uk] + (alpha^2*em4phi*dgt[la, lj, li]*dgt[lk, ll, lm]*gtu[ub, uk]*
   gtu[ui, um]*gtu[uj, ul])/2 - alpha^2*em4phi*dgt[la, lj, li]*
  dgt[lk, ll, lm]*gtu[ub, uj]*gtu[ui, ul]*gtu[uk, um] + 
 (alpha^2*em4phi*dgt[la, lj, li]*dgt[lk, ll, lm]*gtu[ub, uj]*gtu[ui, uk]*
   gtu[ul, um])/2 - alpha*em4phi*gtu[ub, uj]*PD[dalpha[la], lj] + 
 beta[uj]*PD[dbeta[la, ub], lj] - (alpha^2*em4phi*gtu[ub, uj]*gtu[ui, uk]*
   PD[dgt[la, li, lk], lj])/2 + alpha^2*em4phi*gtu[ub, uj]*gtu[ui, uk]*
  PD[dgt[la, lj, li], lk] - (alpha^2*em4phi*gtu[ub, uj]*PD[dphi[la], lj])/6
