4D N=4 CV Specific GL2R Central Charges.nb
Calculations for the specific GL(2,R) choice given in 1812.04236 
1) Demonstrates closure of unhatted transformation laws (for N=2 Vector multiplet)
2) Demonstrates closure of hatted transformation laws 
3) Solves for N=4 vector multiplet central charges.

Setting Rep and Importing Adinkra.wl

$RecursionLimit=Infinity

Rep = N4CV;

<<Adinkra.wl

BuildDate[Adinkra]

FunctionList[Adinkra]

FieldLists

FieldList:= list of all fields

FieldListDynamicalBosons = Join[{A[t,x,y,z],B[t,x,y,z]},Table[Am[mu][t,x,y,z],{mu,0,3}],{Ah[t,x,y,z],Bh[t,x,y,z],Ab[t,x,y,z],Bb[t,x,y,z]}]

FieldListAuxiliaryBosons = Join[{F[t,x,y,z],G[t,x,y,z],d[t,x,y,z]},{Fh[t,x,y,z],Gh[t,x,y,z],Fb[t,x,y,z],Gb[t,x,y,z]}]

FieldListFermions = Join[(Table[Ψ[ii][aa][t,x,y,z],{ii,1,2},{aa,1,4}]//Flatten),(Table[Ψh[ii][aa][t,x,y,z],{ii,1,2},{aa,1,4}]//Flatten)]

FieldList = Join[FieldListDynamicalBosons,FieldListAuxiliaryBosons,FieldListFermions]

FieldListWithDerivatives:= list of all fields and all derivatives of fields, up to 2nd order

FieldListWithDerivatives = Join[FieldList,Table[D[FieldList,Ϛ[mu]],{mu,0,3}]//Flatten,Table[D[FieldList,Ϛ[mu],Ϛ[nu]],{mu,0,3},{nu,0,3}]//Flatten];

FieldProductList := List of all combinations of boson*fermion w/ derivatives on the bosons.

FieldProductsListDynamicalBosons = Table[D[FieldListDynamicalBosons[[ii]],Ϛ[mu],Ϛ[nu]]FieldListFermions[[jhat]],{mu,0,3},{nu,mu,3},{ii,1,Length[FieldListDynamicalBosons]},{jhat,1,Length[FieldListFermions]}]//Flatten;

FieldProductsListAuxiliaryBosons = Table[D[FieldListAuxiliaryBosons[[ii]],Ϛ[mu]]FieldListFermions[[jhat]],{mu,0,3},{ii,1,Length[FieldListAuxiliaryBosons]},{jhat,1,Length[FieldListFermions]}]//Flatten;

FieldProductsList = Join[FieldProductsListDynamicalBosons,FieldProductsListAuxiliaryBosons];

Symmetries, Constraints, and raised indices such as Amup, etc.

Ψup[i][a] := (Ψ^ia)

Do[Ψup[ii][aa][t,x,y,z] = RaiseFermionIndex[Ψ[ii][aa][t,x,y,z]],{ii,1,2},{aa,1,4}]

InverseCmetric//MatrixForm

Table[Ψup[ii][aa][t,x,y,z] ,{ii,1,2},{aa,1,4}]

Ψhup[i][a] := ((OverHat[Ψ])^ia)

Do[Ψhup[ii][aa][t,x,y,z] = RaiseFermionIndex[Ψh[ii][aa][t,x,y,z]],{ii,1,2},{aa,1,4}]

InverseCmetric//MatrixForm

Table[Ψhup[ii][aa][t,x,y,z] ,{ii,1,2},{aa,1,4}]

The Field strength, covariant and contravariant:  Fmn[μ,ν]:=(Subscript[F, μν]) = (∂_μ)(Subscript[A, ν]) - (∂_ν)(Subscript[A, μ])  ,  Fmnup[μ,ν][t,x,y,z] (F^(μν)) = (∂^μ)(A^ν)-(∂^ν)(A^μ)

Do[Fmn[μ,ν][t,x,y,z] = D[Am[ν][t,x,y,z],Ϛ[μ]] - D[Am[μ][t,x,y,z],Ϛ[ν]], {μ,0,3},{ν,0,3}]

Do[Fmnup[μ,ν][t,x,y,z] =  RaiseSTIndex[FmnHolder[μ,ν][t,x,y,z],UP,UP]/.FmnHolder->Fmn,{μ,0,3},{ν,0,3}]

Table[Fmn[μ,ν][t,x,y,z],{μ,0,3},{ν,0,3}]//MatrixForm

Table[Fmnup[μ,ν][t,x,y,z],{μ,0,3},{ν,0,3}]//MatrixForm

Check UD[Field_,Ϛ[μ]] := (∂^μ)Field   ,   Lap[Field_] := (∂^μ)(∂_μ)Field defined in SpaceTime2.nb. Doesn’t matter in which order these are taken when strung together.

D[UD[Am[0][t,x,y,z],Ϛ[0]],Ϛ[0]]

Lap[D[Am[1][t,x,y,z],Ϛ[0]]]

UD[D[Lap[Am[2][t,x,y,z]],Ϛ[0]],Ϛ[0]]

Lap[Lap[A[t,x,y,z]]]

Lap[Ψ[1][1][t,x,y,z]]

Lap[Ψh[2][1][t,x,y,z]]

Lap[Ψhup[2][1][t,x,y,z]]

4D, N=2 Vector transformation laws
Untilded, unhatted transformation laws Dia[1][ii][aa]

Do[Dia[1][ii][aa][A[t,x,y,z]]=  Ψ[ii][aa][t,x,y,z] ,{ii,1,2},{aa,1,4}]

Dia[1][1][1][A[t,x,y,z]]

Dia[1][2][3][A[t,x,y,z]]

Do[Dia[1][ii][aa][B[t,x,y,z]]= I Sum[γ5[[aa,bb]]Ψ[ii][bb][t,x,y,z],{bb,1,4}] ,{ii,1,2},{aa,1,4}]

Dia[1][2][1][B[t,x,y,z]]

Do[Dia[1][ii][aa][F[t,x,y,z]]= Sum[γ[mu][[aa,bb]]D[Ψ[ii][bb][t,x,y,z],Ϛ[mu]],{mu,0,3},{bb,1,4}] ,{ii,1,2},{aa,1,4}]

Dia[1][1][2][F[t,x,y,z]]

Do[Dia[1][ii][aa][G[t,x,y,z]]= I Sum[sigma[3][[ii,jj]]γ5γ[mu][[aa,bb]]D[Ψ[jj][bb][t,x,y,z],Ϛ[mu]],{mu,0,3},{bb,1,4},{jj,1,2}]//Simplify ,{ii,1,2},{aa,1,4}]

Dia[1][1][2][G[t,x,y,z]]

Do[Dia[1][ii][aa][Am[mu][t,x,y,z]] = I Sum[sigma[2][[ii,jj]]γstdown[mu][[aa,bb]]Ψ[jj][bb][t,x,y,z],{jj,1,2},{bb,1,4}],{ii,1,2},{aa,1,4},{mu,0,3}]

Dia[1][2][3][Am[1][t,x,y,z]]

Do[Dia[1][ii][aa][d[t,x,y,z]]= I Sum[sigma[1][[ii,jj]]γ5γ[mu][[aa,bb]]D[Ψ[jj][bb][t,x,y,z],Ϛ[mu]],{mu,0,3},{bb,1,4},{jj,1,2}]//Simplify ,{ii,1,2},{aa,1,4}]

Dia[1][2][2][d[t,x,y,z]]

Do[Dia[1][ii][aa][Ψ[jj][bb][t,x,y,z]] =KroneckerDelta[ii,jj]( Sum[I γdown[mu][[aa,bb]]D[A[t,x,y,z],Ϛ[mu]] -γ5γdown[mu][[aa,bb]]D[B[t,x,y,z],Ϛ[mu]],{mu,0,3}] - I Cmetric[[aa,bb]]F[t,x,y,z]) + sigma[3][[ii,jj]]γ5down[[aa,bb]]G[t,x,y,z] + sigma[1][[ii,jj]]γ5down[[aa,bb]]d[t,x,y,z] +Sum[ (1/4)sigma[2][[ii,jj]]CommuteGammadown[mu,nu][[aa,bb]]Fmn[mu, nu][t,x,y,z],{mu,0,3},{nu,0,3}],{ii,1,2},{aa,1,4},{jj,1,2},{bb,1,4}]

4D, N=2 WF Hypermultiplet (Chiral-Chiral) transformation laws
Untilded, unhatted transformation laws Dia[1][ii][aa]

Do[Dia[1][ii][aa][Ah[t,x,y,z]]=  Sum[sigma[3][[ii,jj]]Ψh[jj][aa][t,x,y,z],{jj,1,2}] ,{ii,1,2},{aa,1,4}]

Dia[1][1][1][Ah[t,x,y,z]]

Dia[1][1][3][Ah[t,x,y,z]]

Do[Dia[1][ii][aa][Bh[t,x,y,z]]= I Sum[γ5[[aa,bb]]Ψh[ii][bb][t,x,y,z],{bb,1,4}] ,{ii,1,2},{aa,1,4}]

Dia[1][2][1][Bh[t,x,y,z]]

Do[Dia[1][ii][aa][Fh[t,x,y,z]]= Sum[γ[mu][[aa,bb]]sigma[3][[ii,jj]]D[Ψh[jj][bb][t,x,y,z],Ϛ[mu]],{mu,0,3},{bb,1,4},{jj,1,2}] ,{ii,1,2},{aa,1,4}]

Dia[1][1][2][Fh[t,x,y,z]]

Do[Dia[1][ii][aa][Gh[t,x,y,z]]= I Sum[γ5γ[mu][[aa,bb]]D[Ψh[ii][bb][t,x,y,z],Ϛ[mu]],{mu,0,3},{bb,1,4}]//Simplify ,{ii,1,2},{aa,1,4}]

Dia[1][1][2][G[t,x,y,z]]

Do[Dia[1][ii][aa][Ab[t,x,y,z]]=  Sum[sigma[1][[ii,jj]]Ψh[jj][aa][t,x,y,z],{jj,1,2}] ,{ii,1,2},{aa,1,4}]

Dia[1][1][1][Ab[t,x,y,z]]

Dia[1][1][3][Ab[t,x,y,z]]

Do[Dia[1][ii][aa][Bb[t,x,y,z]]= -Sum[sigma[2][[ii,jj]]γ5[[aa,bb]]Ψh[jj][bb][t,x,y,z],{bb,1,4},{jj,1,2}] ,{ii,1,2},{aa,1,4}]

Dia[1][2][1][Bb[t,x,y,z]]

Do[Dia[1][ii][aa][Fb[t,x,y,z]]= Sum[γ[mu][[aa,bb]]sigma[1][[ii,jj]]D[Ψh[jj][bb][t,x,y,z],Ϛ[mu]],{mu,0,3},{bb,1,4},{jj,1,2}] ,{ii,1,2},{aa,1,4}]

Dia[1][1][2][Fb[t,x,y,z]]

Do[Dia[1][ii][aa][Gb[t,x,y,z]]= -Sum[sigma[2][[ii,jj]]γ5γ[mu][[aa,bb]]D[Ψh[jj][bb][t,x,y,z],Ϛ[mu]],{mu,0,3},{bb,1,4},{jj,1,2}]//Simplify ,{ii,1,2},{aa,1,4}]

Dia[1][1][2][G[t,x,y,z]]

Do[Dia[1][ii][aa][Ψh[jj][bb][t,x,y,z]] =I sigma[3][[ii,jj]](Sum[γdown[mu][[aa,bb]]D[Ah[t,x,y,z],Ϛ[mu]],{mu,0,3}] - Cmetric[[aa,bb]]Fh[t,x,y,z]) + KroneckerDelta[ii,jj](Sum[-γ5γdown[mu][[aa,bb]]D[Bh[t,x,y,z],Ϛ[mu]],{mu,0,3}] + γ5down[[aa,bb]]Gh[t,x,y,z]) + I sigma[1][[ii,jj]](Sum[γdown[mu][[aa,bb]]D[Ab[t,x,y,z],Ϛ[mu]],{mu,0,3}] - Cmetric[[aa,bb]]Fb[t,x,y,z]) + I sigma[2][[ii,jj]](Sum[-γ5γdown[mu][[aa,bb]]D[Bb[t,x,y,z],Ϛ[mu]],{mu,0,3}] + γ5down[[aa,bb]]Gb[t,x,y,z]),{ii,1,2},{aa,1,4},{jj,1,2},{bb,1,4}]

Dia[1][1][3][Ψh[1][1][t,x,y,z]]

Vector-Chiral hatted transformation laws Dia[0][ii][aa]

Do[Dia[0][ii][aa][A[t,x,y,z]]=  Sum[W[1][ii,jj]Ψh[jj][aa][t,x,y,z] ,{jj,1,2}],{ii,1,2},{aa,1,4}]

Dia[0][1][1][A[t,x,y,z]]

Dia[0][2][3][A[t,x,y,z]]

Do[Dia[0][ii][aa][B[t,x,y,z]]= I Sum[W[2][ii,jj]γ5[[aa,bb]]Ψh[jj][bb][t,x,y,z],{jj,1,2},{bb,1,4}] ,{ii,1,2},{aa,1,4}]

Dia[0][2][1][B[t,x,y,z]]

Do[Dia[0][ii][aa][F[t,x,y,z]]= Sum[W[3][ii,jj]γ[mu][[aa,bb]]D[Ψh[jj][bb][t,x,y,z],Ϛ[mu]],{jj,1,2},{mu,0,3},{bb,1,4}] ,{ii,1,2},{aa,1,4}]

Dia[0][1][2][F[t,x,y,z]]

Do[Dia[0][ii][aa][G[t,x,y,z]]= I Sum[W[4][ii,jj]γ5γ[mu][[aa,bb]]D[Ψh[jj][bb][t,x,y,z],Ϛ[mu]],{mu,0,3},{bb,1,4},{jj,1,2}]//Simplify ,{ii,1,2},{aa,1,4}]

Dia[0][1][2][G[t,x,y,z]]

Do[Dia[0][ii][aa][Am[mu][t,x,y,z]] = Sum[W[5][ii,jj]γstdown[mu][[aa,bb]]Ψh[jj][bb][t,x,y,z],{jj,1,2},{bb,1,4}],{ii,1,2},{aa,1,4},{mu,0,3}]

Dia[0][2][3][Am[1][t,x,y,z]]

Do[Dia[0][ii][aa][d[t,x,y,z]]= I Sum[W[6][ii,jj]γ5γ[mu][[aa,bb]]D[Ψh[jj][bb][t,x,y,z],Ϛ[mu]],{mu,0,3},{bb,1,4},{jj,1,2}]//Simplify ,{ii,1,2},{aa,1,4}]

Dia[0][2][2][d[t,x,y,z]]

Do[Dia[0][ii][aa][Ψ[jj][bb][t,x,y,z]] = I W[7][ii,jj]Sum[γdown[mu][[aa,bb]]D[Ah[t,x,y,z],Ϛ[mu]],{mu,0,3}] - I W[8][ii,jj] Cmetric[[aa,bb]]Fh[t,x,y,z] - W[9][ii,jj]Sum[γ5γdown[mu][[aa,bb]]D[Bh[t,x,y,z],Ϛ[mu]],{mu,0,3}] +W[10][ii,jj] γ5down[[aa,bb]]Gh[t,x,y,z]  +  I W[11][ii,jj]Sum[γdown[mu][[aa,bb]]D[Ab[t,x,y,z],Ϛ[mu]],{mu,0,3}] - I W[12][ii,jj] Cmetric[[aa,bb]]Fb[t,x,y,z] - W[13][ii,jj]Sum[γ5γdown[mu][[aa,bb]]D[Bb[t,x,y,z],Ϛ[mu]],{mu,0,3}] +W[14][ii,jj] γ5down[[aa,bb]]Gb[t,x,y,z] ,{ii,1,2},{aa,1,4},{jj,1,2},{bb,1,4}]

Dia[0][1][3][Ψ[1][1][t,x,y,z]]

Do[Dia[0][ii][aa][Ah[t,x,y,z]]=  Sum[X[1][ii,jj]Ψ[jj][aa][t,x,y,z],{jj,1,2}] ,{ii,1,2},{aa,1,4}]

Dia[0][1][1][Ah[t,x,y,z]]

Dia[0][1][3][Ah[t,x,y,z]]

Do[Dia[0][ii][aa][Bh[t,x,y,z]]= I Sum[X[2][ii,jj]γ5[[aa,bb]]Ψ[jj][bb][t,x,y,z],{jj,1,2},{bb,1,4}] ,{ii,1,2},{aa,1,4}]

Dia[0][2][1][Bh[t,x,y,z]]

Do[Dia[0][ii][aa][Fh[t,x,y,z]]= Sum[X[3][ii,jj]γ[mu][[aa,bb]]D[Ψ[jj][bb][t,x,y,z],Ϛ[mu]],{jj,1,2},{mu,0,3},{bb,1,4}] ,{ii,1,2},{aa,1,4}]

Dia[0][1][2][Fh[t,x,y,z]]

Do[Dia[0][ii][aa][Gh[t,x,y,z]]= I Sum[X[4][ii,jj]γ5γ[mu][[aa,bb]]D[Ψ[jj][bb][t,x,y,z],Ϛ[mu]],{jj,1,2},{mu,0,3},{bb,1,4}]//Simplify ,{ii,1,2},{aa,1,4}]

Dia[0][1][2][Gh[t,x,y,z]]

Do[Dia[0][ii][aa][Ab[t,x,y,z]]=  Sum[X[5][ii,jj]Ψ[jj][aa][t,x,y,z],{jj,1,2}] ,{ii,1,2},{aa,1,4}]

Dia[0][1][1][Ab[t,x,y,z]]

Dia[0][1][3][Ab[t,x,y,z]]

Do[Dia[0][ii][aa][Bb[t,x,y,z]]= I Sum[X[6][ii,jj]γ5[[aa,bb]]Ψ[jj][bb][t,x,y,z],{jj,1,2},{bb,1,4}] ,{ii,1,2},{aa,1,4}]

Dia[0][2][1][Bb[t,x,y,z]]

Do[Dia[0][ii][aa][Fb[t,x,y,z]]= Sum[X[7][ii,jj]γ[mu][[aa,bb]]D[Ψ[jj][bb][t,x,y,z],Ϛ[mu]],{jj,1,2},{mu,0,3},{bb,1,4}] ,{ii,1,2},{aa,1,4}]

Dia[0][1][2][Fb[t,x,y,z]]

Do[Dia[0][ii][aa][Gb[t,x,y,z]]= I Sum[X[8][ii,jj]γ5γ[mu][[aa,bb]]D[Ψ[jj][bb][t,x,y,z],Ϛ[mu]],{jj,1,2},{mu,0,3},{bb,1,4}]//Simplify ,{ii,1,2},{aa,1,4}]

Dia[0][1][2][Gb[t,x,y,z]]

Do[Dia[0][ii][aa][Ψh[jj][bb][t,x,y,z]] =Sum[I X[9][ii,jj] γdown[mu][[aa,bb]]D[A[t,x,y,z],Ϛ[mu]] -X[10][ii,jj]γ5γdown[mu][[aa,bb]]D[B[t,x,y,z],Ϛ[mu]],{mu,0,3}] - I X[11][ii,jj]Cmetric[[aa,bb]]F[t,x,y,z] + X[12][ii,jj]γ5down[[aa,bb]]G[t,x,y,z] + X[13][ii,jj]γ5down[[aa,bb]]d[t,x,y,z] -I (1/4)Sum[ X[14][ii,jj]CommuteGammadown[mu,nu][[aa,bb]]Fmn[mu, nu][t,x,y,z],{mu,0,3},{nu,0,3}],{ii,1,2},{aa,1,4},{jj,1,2},{bb,1,4}]

Dia[0][1][3][Ψh[1][1][t,x,y,z]]

Lagrangian Invariance:  ℒ = (Subscript[ℒ, (2VS) ])+(Subscript[ℒ, (2TS)])

Lagrangian Invariance: DL[tt] w/ tt= tilde index (0 = hat, 1 = no tilde or hat) (just the CC and VS portions)
(Subscript[D, a]^i)ℒ = 0  = ((Subscript[D, a]^i)A) □A + ((Subscript[D, a]^i)B) □B +F (Subscript[D, a]^i)F +G (Subscript[D, a]^i)G+d (Subscript[D, a]^i)d+ ((Subscript[D, a]^i)(Subscript[A, ν]))(∂_μ)(F^(μν)) - i (Subscript[Ψ, b]^j)((((γ^μ)))^bc)(∂_μ)((Subscript[D, a]^i)(Subscript[Ψ, c]^j)) + ((Subscript[D, a]^i)(OverHat[A])) □(OverHat[A]) + ((Subscript[D, a]^i)(OverHat[B])) □(OverHat[B]) +(OverHat[F]) (Subscript[D, a]^i)(OverHat[F])  +(OverHat[G])  (Subscript[D, a]^i)(OverHat[G]) + ((Subscript[D, a]^i)(Overscript[A, ⌣])) □(Overscript[A, ⌣]) + ((Subscript[D, a]^i)(Overscript[B, ⌣])) □(Overscript[B, ⌣]) +(Overscript[F, ⌣]) (Subscript[D, a]^i)(Overscript[F, ⌣])  +(Overscript[G, ⌣])  (Subscript[D, a]^i)(Overscript[G, ⌣])  - i (Subscript[(OverHat[Ψ]), b]^j)((((γ^μ)))^bc)(∂_μ)((Subscript[D, a]^i)(Subscript[(OverHat[Ψ]), c]^j)) 
(Subscript[(OverHat[D]), a]^i)ℒ = 0  = ((Subscript[(OverHat[D]), a]^i)A) □A + ((Subscript[(OverHat[D]), a]^i)B) □B +F (Subscript[(OverHat[D]), a]^i)F +G (Subscript[(OverHat[D]), a]^i)G+d (Subscript[D, a]^i)d+ ((Subscript[(OverHat[D]), a]^i)(Subscript[A, ν]))(∂_μ)(F^(μν)) - i (Subscript[Ψ, b]^j)((((γ^μ)))^bc)(∂_μ)((Subscript[(OverHat[D]), a]^i)(Subscript[Ψ, c]^j)) + ((Subscript[(OverHat[D]), a]^i)(OverHat[A])) □(OverHat[A]) + ((Subscript[(OverHat[D]), a]^i)(OverHat[B])) □(OverHat[B]) +(OverHat[F]) (Subscript[(OverHat[D]), a]^i)(OverHat[F])  +(OverHat[G])  (Subscript[D, a]^i)(OverHat[G]) + ((Subscript[(OverHat[D]), a]^i)(Overscript[A, ⌣])) □(Overscript[A, ⌣]) + ((Subscript[(OverHat[D]), a]^i)(Overscript[B, ⌣])) □(Overscript[B, ⌣]) +(Overscript[F, ⌣]) (Subscript[(OverHat[D]), a]^i)(Overscript[F, ⌣])  +(Overscript[G, ⌣])  (Subscript[(OverHat[D]), a]^i)(Overscript[G, ⌣])  - i (Subscript[(OverHat[Ψ]), b]^j)((((γ^μ)))^bc)(∂_μ)((Subscript[(OverHat[D]), a]^i)(Subscript[(OverHat[Ψ]), c]^j))

Integration by parts commands: Put all derivatives on the bosons

ByParts={ Derivative[0,0,0,1][Ψ[ii_][xp__]][t,x,y,z]efield__[t,x,y,z]-> -  Ψ[ii][xp][t,x,y,z]D[efield[t,x,y,z],z],
Derivative[0,0,1,0][Ψ[ii_][xp__]][t,x,y,z]efield__[t,x,y,z]-> - Ψ[ii][xp][t,x,y,z]D[ efield[t,x,y,z],y],Derivative[0,1,0,0][Ψ[ii_][xp__]][t,x,y,z]efield__[t,x,y,z]-> -  Ψ[ii][xp][t,x,y,z]D[ efield[t,x,y,z],x],Derivative[1,0,0,0][Ψ[ii_][xp__]][t,x,y,z]efield__[t,x,y,z]-> -  Ψ[ii][xp][t,x,y,z]D[ efield[t,x,y,z],t],Derivative[0,0,0,1][Ψh[ii_][xp__]][t,x,y,z]efield__[t,x,y,z]-> -  Ψh[ii][xp][t,x,y,z]D[efield[t,x,y,z],z],
Derivative[0,0,1,0][Ψh[ii_][xp__]][t,x,y,z]efield__[t,x,y,z]-> - Ψh[ii][xp][t,x,y,z]D[ efield[t,x,y,z],y],Derivative[0,1,0,0][Ψh[ii_][xp__]][t,x,y,z]efield__[t,x,y,z]-> -  Ψh[ii][xp][t,x,y,z]D[ efield[t,x,y,z],x],Derivative[1,0,0,0][Ψh[ii_][xp__]][t,x,y,z]efield__[t,x,y,z]-> -  Ψh[ii][xp][t,x,y,z]D[ efield[t,x,y,z],t]}

DLA[tt] = {((Subscript[D, a]^i)A) □A,((Subscript[(OverTilde[D]), a]^i)A) □A}

Do[DLA[tt] = Table[Dia[tt][ii][aa][A[t,x,y,z]]Lap[A[t,x,y,z]],{ii,1,2},{aa,1,4}],{tt,0,1}]

DLB[tt] = {((Subscript[D, a]^i)B) □B,((Subscript[(OverTilde[D]), a]^i)B) □B}

Do[DLB[tt] = Table[Dia[tt][ii][aa][B[t,x,y,z]]Lap[B[t,x,y,z]],{ii,1,2},{aa,1,4}],{tt,0,1}]

DLF[tt] = {((Subscript[D, a]^i)F) F,((Subscript[(OverTilde[D]), a]^i)F) F}

Do[DLF[tt] = Expand[Table[Dia[tt][ii][aa][F[t,x,y,z]]F[t,x,y,z],{ii,1,2},{aa,1,4}]]/.ByParts,{tt,0,1}]

DLG[tt] = {((Subscript[D, a]^i)G) G,((Subscript[(OverTilde[D]), a]^i)G) G}

Do[DLG[tt] =( Table[Dia[tt][ii][aa][G[t,x,y,z]]G[t,x,y,z],{ii,1,2},{aa,1,4}]//Expand)/.ByParts,{tt,0,1}]

DLd[tt] = {((Subscript[D, a]^i)d) d,((Subscript[(OverTilde[D]), a]^i)d) d}

Do[DLd[tt] = Expand[Table[Dia[tt][ii][aa][d[t,x,y,z]]d[t,x,y,z],{ii,1,2},{aa,1,4}]]/.ByParts,{tt,0,1}]

DLAm[tt] = {((Subscript[D, a]^i)(Subscript[A, ν]))(∂_μ)(F^(μν)),((Subscript[(OverTilde[D]), a]^i)(Subscript[A, ν]))(∂_μ)(F^(μν))}

Do[DLAm[tt] = Table[Sum[Dia[tt][ii][aa][Am[nu][t,x,y,z]]D[Fmnup[mu,nu][t,x,y,z],Ϛ[mu]],{mu,0,3},{nu,0,3}],{ii,1,2},{aa,1,4}],{tt,0,1}]

DLPsi[tt] =  {- i (Subscript[(OverTilde[Ψ]), b]^j)((((γ^μ)))^bc)(∂_μ)((Subscript[D, a])(Subscript[(OverTilde[Ψ]), c]^j)),- i (Subscript[(OverTilde[Ψ]), b]^j)((((γ^μ)))^bc)(∂_μ)((Subscript[(OverTilde[D]), a]^i)(Subscript[(OverTilde[Ψ]), c]^j))}

Do[DLPsi[tt] = Table[ - I Sum[Ψ[jj][bb][t,x,y,z]γup[mu][[bb,cc]]D[Dia[tt][ii][aa][Ψ[jj][cc][t,x,y,z]],Ϛ[mu]],{jj,1,2},{bb,1,4},{cc,1,4},{mu,0,3}],{ii,1,2},{aa,1,4}]/.ByParts,{tt,0,1}]

DLPsi[1][[2,2]]

DLAh[tt] = {((Subscript[D, a]^i)Ah) □Ah,((Subscript[(OverTilde[D]), a]^i)Ah) □Ah}

Do[DLAh[tt] = Table[Dia[tt][ii][aa][Ah[t,x,y,z]]Lap[Ah[t,x,y,z]],{ii,1,2},{aa,1,4}],{tt,0,1}]

DLBh[tt] = {((Subscript[D, a]^i)Bh) □Bh,((Subscript[(OverTilde[D]), a]^i)Bh) □Bh}

Do[DLBh[tt] = Table[Dia[tt][ii][aa][Bh[t,x,y,z]]Lap[Bh[t,x,y,z]],{ii,1,2},{aa,1,4}],{tt,0,1}]

DLFh[tt] = {((Subscript[D, a]^i)Fh) Fh,((Subscript[(OverTilde[D]), a]^i)Fh) Fh}

Do[DLFh[tt] = Expand[Table[Dia[tt][ii][aa][Fh[t,x,y,z]]Fh[t,x,y,z],{ii,1,2},{aa,1,4}]]/.ByParts,{tt,0,1}]

DLGh[tt] = {((Subscript[D, a]^i)Gh) Gh,((Subscript[(OverTilde[D]), a]^i)Gh) Gh}

Do[DLGh[tt] = (Table[Dia[tt][ii][aa][Gh[t,x,y,z]]Gh[t,x,y,z],{ii,1,2},{aa,1,4}]//Expand)/.ByParts,{tt,0,1}]

DLAb[tt] = {((Subscript[D, a]^i)Ab) □Ab,((Subscript[(OverTilde[D]), a]^i)Ab) □Ab}

Do[DLAb[tt] = Table[Dia[tt][ii][aa][Ab[t,x,y,z]]Lap[Ab[t,x,y,z]],{ii,1,2},{aa,1,4}],{tt,0,1}]

DLBb[tt] = {((Subscript[D, a]^i)Bb) □Bb,((Subscript[(OverTilde[D]), a]^i)Bb) □Bb}

Do[DLBb[tt] = Table[Dia[tt][ii][aa][Bb[t,x,y,z]]Lap[Bb[t,x,y,z]],{ii,1,2},{aa,1,4}],{tt,0,1}]

DLFb[tt] = {((Subscript[D, a]^i)Fb) Fb,((Subscript[(OverTilde[D]), a]^i)Fb) Fb}

Do[DLFb[tt] = Expand[Table[Dia[tt][ii][aa][Fb[t,x,y,z]]Fb[t,x,y,z],{ii,1,2},{aa,1,4}]]/.ByParts,{tt,0,1}]

DLGb[tt] = {((Subscript[D, a]^i)Gb) Gb,((Subscript[(OverTilde[D]), a]^i)Gb) Gb}

Do[DLGb[tt] = (Table[Dia[tt][ii][aa][Gb[t,x,y,z]]Gb[t,x,y,z],{ii,1,2},{aa,1,4}]//Expand)/.ByParts,{tt,0,1}]

DLPsih[tt] =  {- i (Subscript[(OverHat[Ψ]), b]^j)((((γ^μ)))^bc)(∂_μ)((Subscript[D, a])(Subscript[(OverHat[Ψ]), c]^j)),- i (Subscript[(OverHat[Ψ]), b]^j)((((γ^μ)))^bc)(∂_μ)((Subscript[(OverTilde[D]), a]^i)(Subscript[(OverHat[Ψ]), c]^j))}

Do[DLPsih[tt] = Table[ - I Sum[Ψh[jj][bb][t,x,y,z]γup[mu][[bb,cc]]D[Dia[tt][ii][aa][Ψh[jj][cc][t,x,y,z]],Ϛ[mu]],{jj,1,2},{bb,1,4},{cc,1,4},{mu,0,3}],{ii,1,2},{aa,1,4}]/.ByParts,{tt,0,1}]

DLPsih[0]

Lagrangian invariance

Do[DLTotal[tt]=(DLA[tt] +DLB[tt] +DLF[tt] +DLG[tt] +DLd[tt] +DLAm[tt] + DLPsi[tt] +DLAh[tt] +DLBh[tt] +DLFh[tt] +DLGh[tt]  +DLAb[tt] +DLBb[tt] +DLFb[tt] +DLGb[tt]+ DLPsih[tt] ),{tt,0,1}];

DLTotal[1]

Do[DLFinal[tt] = Collect[DLTotal[tt],FieldListWithDerivatives,Simplify],{tt,0,1}]

Lagrangian invariance: untilded DLFinal[1]

DLFinal[1]

Lagrangian Parameters

XList = Table[X[Index][ii,jj],{Index,1,14},{ii,1,2},{jj,1,2}]//Flatten

CoefficientEquations[0] = Table[Coefficient[DLFinal[0][[ii,aa]],FieldProductsList[[FieldProductsIndex]]]==0,{ii,1,2},{aa,1,4},{FieldProductsIndex,1,Length[FieldProductsList]}]//Flatten;

LagrangianParameters = Solve[CoefficientEquations[0],XList]//Flatten

Lagrangian Invariance: tilded

DLFinal[0]/.LagrangianParameters

We have 
X[1][i,j] =W[7][i,j] 
X[2][i,j] =W[9][i,j] 
X[3][i,j] =W[8][i,j] 
X[4][i,j]=W[10][i,j] 
X[5][i,j] =W[11][i,j] 
X[6][i,j] =W[13][i,j] 
X[7][i,j] =W[12][i,j]  
X[8][i,j] =W[14][i,j]  
X[9][i,j] =W[1][i,j]  
X[10][i,j] =W[2][i,j]  
X[11][i,j] =W[3][i,j] 
X[12][i,j] =W[4][i,j] 
X[13][i,j] =W[6][i,j] 
X[14][i,j] =W[5][i,j]

XWPairings = {{1,7},{2,9},{3,8},{4,10},{5,11},{6,13},{7,12},{8,14},{9,1},{10,2},{11,3},{12,4},{13,6},{14,5}};

Table[(X[XWPairings[[uvpair,1]]][ii,jj]/.LagrangianParameters)==W[XWPairings[[uvpair,2]]][ii,jj],{uvpair,1,14},{ii,1,2},{jj,1,2}] == Table[True,{uvpair,1,14},{ii,1,2},{jj,1,2}]

Lagrangian Summary

DLFinal[1]==(DLFinal[0]/.LagrangianParameters)=={{0,0,0,0},{0,0,0,0}}

Algebra {(Subscript[D, a]^i),(Subscript[D, b]^j)} and {(Subscript[(OverHat[D]), a]^i),(Subscript[(OverHat[D]), b]^j)}

DD

SUSYTransformFields[Operator_]:=Table[FieldList[[FieldIndex]]->Operator[FieldList[[FieldIndex]]],{FieldIndex,1,Length[FieldList]}]

SUSYTransformJustB[Operator_]:=Table[FieldList[[FieldIndex]]->Operator[FieldList[[FieldIndex]]],{FieldIndex,1,4}]

Do[SUSYTransform[Dia[tt][ii][aa]] = Join[SUSYTransformFields[Dia[tt][ii][aa]],D[SUSYTransformFields[Dia[tt][ii][aa]],t],D[SUSYTransformFields[Dia[tt][ii][aa]],x],D[SUSYTransformFields[Dia[tt][ii][aa]],y],D[SUSYTransformFields[Dia[tt][ii][aa]],z]],{tt,0,1},{ii,1,2},{aa,1,4}]

DD[LeftOperator_,RightOperator_][Field_]:=RightOperator[Field]/.SUSYTransform[LeftOperator]

AlgebraDDTerms[LeftOperator_,RightOperator_][Field_]:= DD[LeftOperator,RightOperator][Field] + DD[RightOperator,LeftOperator][Field]

AlgebraDD[ttl_,ttr_][Field_] := Table[AlgebraDDTerms[Dia[ttl][ii][aa],Dia[ttr][jj][bb]][Field],{ii,1,2},{aa,1,4},{jj,1,2},{bb,1,4}]  (* Convention is left D upper index, left D lower index, right D upper index, right D lower index *)

HoloraumyDDTerms[LeftOperator_,RightOperator_][Field_]:= DD[LeftOperator,RightOperator][Field] - DD[RightOperator,LeftOperator][Field]

HoloraumyDD[ttl_,ttr_][Field_]:=Table[HoloraumyDDTerms[Dia[ttl][ii][aa],Dia[ttr][jj][bb]][Field],{ii,1,2},{aa,1,4},{jj,1,2},{bb,1,4}]

ClosureTerms[ii_,jj_][aa_,bb_][Field_]:=Sum[2 I KroneckerDelta[ii,jj]γdown[mu][[aa,bb]]D[Field,Ϛ[mu]],{mu,0,3}]

Closure[Field_]:=Table[ClosureTerms[ii,jj][aa,bb][Field],{ii,1,2},{aa,1,4},{jj,1,2},{bb,1,4}]

Zeros = Table[0,{ii,1,2},{aa,1,4},{jj,1,2},{bb,1,4}];

Test it out

AlgebraDD[tt1,tt2][Ψh[1][1][t,x,y,z]]

Closure[A[t,x,y,z]]

AlgebraDD[1,1][A[t,x,y,z]]

Table[ClosureTerms[ii,jj][aa,bb][Ψ[ll][kk][t,x,y,z]],{ii,1,2},{aa,1,4},{jj,1,2},{bb,1,4},{ll,1,2},{kk,1,4}]

AllgebraDDResults[tt1,tt2][FieldIndex] = True if the algebra closes fo Field FieldIndex, False if not. tt1,tt2 are tilded indices: , tt1,tt2=0 for hatted, tt1,tt2=1 for untilded and unhatted. set tt1 = tt2 = 1 for the next few unhatted subsections.

tt1 = 1;
tt2 = 1;

4D, N=2 Vector Algebra, both unhatted {(Subscript[D, a]^i),(Subscript[D, b]^j)}

{(Subscript[D, a]^i),(Subscript[D, b]^j)}A= 2 i (Subscript[(δ^ij)((γ^μ)), ab])(∂_μ)A , AlgebraDDResult[1,1][FieldNumber], FieldNumber = 1

AlgebraDDResult[tt1,tt2][1]=(AlgebraDD[tt1,tt2][A[t,x,y,z]] ==Closure[A[t,x,y,z]])//Simplify

{(Subscript[D, a]^i),(Subscript[D, b]^j)}B= 2 i (Subscript[(δ^ij)((γ^μ)), ab])(∂_μ)B , FieldNumber = 2

AlgebraDDResult[tt1,tt2][2]=(AlgebraDD[tt1,tt2][B[t,x,y,z]] ==Closure[B[t,x,y,z]])//Simplify

{(Subscript[D, a]^i),(Subscript[D, b]^j)}F= 2 i (Subscript[(δ^ij)((γ^μ)), ab])(∂_μ)F, FieldNumber = 3

AlgebraDDResult[tt1,tt2][3]=(AlgebraDD[tt1,tt2][F[t,x,y,z]] ==Closure[F[t,x,y,z]])//Simplify

{(Subscript[D, a]^i),(Subscript[D, b]^j)}G= 2 i (Subscript[(δ^ij)((γ^μ)), ab])(∂_μ)G, FieldNumber = 4

AlgebraDDResult[tt1,tt2][4]=(AlgebraDD[tt1,tt2][G[t,x,y,z]] ==Closure[G[t,x,y,z]])//Simplify

{(Subscript[D, a]^i),(Subscript[D, b]^j)}d= 2 i (Subscript[(δ^ij)((γ^μ)), ab])(∂_μ)d, FieldNumber = 5

AlgebraDDResult[tt1,tt2][5]=(AlgebraDD[tt1,tt2][d[t,x,y,z]] ==Closure[d[t,x,y,z]])//Simplify

{(Subscript[D, a]^i),(Subscript[D, b]^j)}(Subscript[Ψ, c]^k)= 2 i (Subscript[(δ^ij)((γ^μ)), ab])(∂_μ)(Subscript[Ψ, c]^k), FieldNumber = 6

AlgebraDDResult[tt1,tt2][6]=Table[AlgebraDD[tt1,tt2][Ψ[kk][cc][t,x,y,z]],{kk,1,2},{cc,1,4}] ==Table[Closure[Ψ[kk][cc][t,x,y,z]],{kk,1,2},{cc,1,4}]//Simplify

{(Subscript[D, a]^i),(Subscript[D, b]^j)}(Subscript[A, μ])= 2 i (Subscript[(δ^ij)((γ^ν)), ab])(∂_ν)(Subscript[A, μ]) -(∂_μ)(Subscript[λ, ab]^ij) , FieldNumber = 7
(Subscript[λ, ab]^ij)= 2 i (Subscript[(δ^ij)((γ^ν)), ab])(Subscript[A, ν]) - 2 i ((((σ^2)))^ij)( i (Subscript[C, ab])A - (Subscript[((γ^5)), ab])B)

λ[tt1,tt2][ii_,jj_][aa_,bb_][t,x,y,z] := 2 I KroneckerDelta[ii,jj]Sum[γdown[nu][[aa,bb]]Am[nu][t,x,y,z],{nu,0,3}] - 2I sigma[2][[ii,jj]](I Cmetric[[aa,bb]] A[t,x,y,z] - γ5down[[aa,bb]]B[t,x,y,z])

λTable[tt1,tt2][t,x,y,z] = Table[λ[tt1,tt2][ii,jj][aa,bb][t,x,y,z],{ii,1,2},{aa,1,4},{jj,1,2},{bb,1,4}];

AlgebraDDResult[tt1,tt2][7]=(Table[AlgebraDD[tt1,tt2][Am[mu][t,x,y,z]],{mu,0,3}] ==Table[Closure[Am[mu][t,x,y,z]]- D[λTable[tt1,tt2][t,x,y,z],Ϛ[mu]],{mu,0,3}] )//Simplify

4D, N=2 Chiral-Chiral, both unhatted {(Subscript[D, a]^i),(Subscript[D, b]^j)}

, FieldNumber = 8

AlgebraDDResult[tt1,tt2][8]=(AlgebraDD[tt1,tt2][Ah[t,x,y,z]] ==Closure[Ah[t,x,y,z]] + Table[I sigma[2][[ii,jj]] 2 I Cmetric[[aa,bb]]Fb[t,x,y,z],{ii,1,2},{aa,1,4},{jj,1,2},{bb,1,4}])//Simplify

, FieldNumber = 9

AlgebraDDResult[tt1,tt2][9]=(AlgebraDD[tt1,tt2][Ab[t,x,y,z]] ==Closure[Ab[t,x,y,z]] - Table[I sigma[2][[ii,jj]] 2 I Cmetric[[aa,bb]]Fh[t,x,y,z],{ii,1,2},{aa,1,4},{jj,1,2},{bb,1,4}])//Simplify

, FieldNumber = 10

AlgebraDDResult[tt1,tt2][10]=(AlgebraDD[tt1,tt2][Bh[t,x,y,z]] ==Closure[Bh[t,x,y,z]] +  Table[I sigma[2][[ii,jj]] 2 I Cmetric[[aa,bb]]Gb[t,x,y,z],{ii,1,2},{aa,1,4},{jj,1,2},{bb,1,4}])//Simplify

, FieldNumber = 11

AlgebraDDResult[tt1,tt2][11]=(AlgebraDD[tt1,tt2][Bb[t,x,y,z]] ==Closure[Bb[t,x,y,z]] -  Table[I sigma[2][[ii,jj]] 2 I Cmetric[[aa,bb]]Gh[t,x,y,z],{ii,1,2},{aa,1,4},{jj,1,2},{bb,1,4}])//Simplify

, FieldNumber = 12

AlgebraDDResult[tt1,tt2][12]=(AlgebraDD[tt1,tt2][Fh[t,x,y,z]] ==Closure[Fh[t,x,y,z]] + Table[I sigma[2][[ii,jj]]2 I Cmetric[[aa,bb]]Lap[Ab[t,x,y,z]],{ii,1,2},{aa,1,4},{jj,1,2},{bb,1,4}])//Simplify

, FieldNumber = 13

AlgebraDDResult[tt1,tt2][13]=(AlgebraDD[tt1,tt2][Fb[t,x,y,z]] ==Closure[Fb[t,x,y,z]] - Table[I sigma[2][[ii,jj]]2 I Cmetric[[aa,bb]]Lap[Ah[t,x,y,z]],{ii,1,2},{aa,1,4},{jj,1,2},{bb,1,4}])//Simplify

, FieldNumber = 14

AlgebraDDResult[tt1,tt2][14]=(AlgebraDD[tt1,tt2][Gh[t,x,y,z]] ==Closure[Gh[t,x,y,z]] + Table[I sigma[2][[ii,jj]]2 I Cmetric[[aa,bb]]Lap[Bb[t,x,y,z]],{ii,1,2},{aa,1,4},{jj,1,2},{bb,1,4}])//Simplify

, FieldNumber = 15

AlgebraDDResult[tt1,tt2][15]=(AlgebraDD[tt1,tt2][Gb[t,x,y,z]] ==Closure[Gb[t,x,y,z]] - Table[I sigma[2][[ii,jj]]2 I Cmetric[[aa,bb]]Lap[Bh[t,x,y,z]],{ii,1,2},{aa,1,4},{jj,1,2},{bb,1,4}])//Simplify

, FieldNumber = 16

Closure[Ψh[1][2][t,x,y,z]][[2,4,2,4]]

Table[Closure[Ψh[kk][cc][t,x,y,z]],{kk,1,2},{cc,1,4}][[1,2,2,4,2,4]]

AlgebraDDResult[tt1,tt2][16]=(Table[AlgebraDD[tt1,tt2][Ψh[kk][cc][t,x,y,z]],{kk,1,2},{cc,1,4}] ==Table[Closure[Ψh[kk][cc][t,x,y,z]],{kk,1,2},{cc,1,4}] - Table[sigma[2][[ii,jj]]Sum[sigma[2][[kk,rr]]2I Cmetric[[aa,bb]]γ[mu][[cc,dd]]D[Ψh[rr][dd][t,x,y,z],Ϛ[mu]],{mu,0,3},{dd,1,4},{rr,1,2}],{kk,1,2},{cc,1,4},{ii,1,2},{aa,1,4},{jj,1,2},{bb,1,4}])//Simplify

4D, N=2 Summary, unhatted (tt1=tt2=1)

Table[AlgebraDDResult[tt1,tt2][Index],{Index,1,16}]==Table[True,{Index,1,16}]
tt1
tt2

tt1 = 0, tt2 = 0 for both hatted D’s

Clear[tt1,tt2];
tt1 = 0;
tt2 = 0;

4D, N=2 Vector Algebra, both hatted {(Subscript[(OverHat[D]), a]^i),(Subscript[(OverHat[D]), b]^j)}

Choose

which implies for Lagrangian invariance

SetW = Table[{W[1][ii,jj]-> KroneckerDelta[ii,jj], W[2][ii,jj]->KroneckerDelta[ii,jj],  W[3][ii,jj]->KroneckerDelta[ii,jj], W[4][ii,jj]->sigma[3][[ii,jj]], W[5][ii,jj]->I sigma[2][[ii,jj]], W[6][ii,jj]->sigma[1][[ii,jj]], W[7][ii,jj]->sigma[3][[ii,jj]], W[8][ii,jj]->sigma[3][[ii,jj]], W[9][ii,jj]->KroneckerDelta[ii,jj], W[10][ii,jj]->KroneckerDelta[ii,jj], W[11][ii,jj]->sigma[1][[ii,jj]], W[12][ii,jj]->sigma[1][[ii,jj]], W[13][ii,jj]->I sigma[2][[ii,jj]], W[14][ii,jj]->I sigma[2][[ii,jj]]},{ii,1,2},{jj,1,2}]//Flatten;

Table[PrintSigmaProduct[Table[W[nn][ii,jj],{ii,1,2},{jj,1,2}]/.LagrangianParameters/.SetW],{nn,1,14}]

Table[PrintSigmaProduct[Table[X[nn][ii,jj],{ii,1,2},{jj,1,2}]/.LagrangianParameters/.SetW],{nn,1,14}]

{(Subscript[(OverHat[D]), a]^i),(Subscript[(OverHat[D]), b]^j)}A= 2 i (Subscript[(δ^ij)((γ^μ)), ab])(∂_μ)A , FieldNumber = 1

LocalField = A[t,x,y,z];
LocalFieldNumber = 1;

AlgebraDDResult[tt1,tt2][LocalFieldNumber]=(Collect[(AlgebraDD[tt1,tt2][LocalField] ==Closure[LocalField])/.LagrangianParameters/.SetW,FieldListWithDerivatives,Simplify])

Clear[LocalField,LocalFieldNumber]

{(Subscript[(OverHat[D]), a]^i),(Subscript[(OverHat[D]), b]^j)}B= 2 i (Subscript[(δ^ij)((γ^μ)), ab])(∂_μ)B , FieldNumber = 2

LocalField = B[t,x,y,z];
LocalFieldNumber = 2;

AlgebraDDResult[tt1,tt2][LocalFieldNumber]=(Collect[(AlgebraDD[tt1,tt2][LocalField] ==Closure[LocalField])/.LagrangianParameters/.SetW,FieldListWithDerivatives,Simplify])

Clear[LocalField,LocalFieldNumber]

{(Subscript[(OverHat[D]), a]^i),(Subscript[(OverHat[D]), b]^j)}F= 2 i (Subscript[(δ^ij)((γ^μ)), ab])(∂_μ)F , FieldNumber = 3

LocalField = F[t,x,y,z];
LocalFieldNumber = 3;

AlgebraDDResult[tt1,tt2][LocalFieldNumber]=(Collect[(AlgebraDD[tt1,tt2][LocalField] ==Closure[LocalField])/.LagrangianParameters/.SetW,FieldListWithDerivatives,Simplify])

Clear[LocalField,LocalFieldNumber]

{(Subscript[(OverHat[D]), a]^i),(Subscript[(OverHat[D]), b]^j)}G= 2 i (Subscript[(δ^ij)((γ^μ)), ab])(∂_μ)G , FieldNumber = 4

LocalField = G[t,x,y,z];
LocalFieldNumber = 4;

AlgebraDDResult[tt1,tt2][LocalFieldNumber]=(Collect[(AlgebraDD[tt1,tt2][LocalField] ==Closure[LocalField])/.LagrangianParameters/.SetW,FieldListWithDerivatives,Simplify])

Clear[LocalField,LocalFieldNumber]

{(Subscript[(OverHat[D]), a]^i),(Subscript[(OverHat[D]), b]^j)}d= 2 i (Subscript[(δ^ij)((γ^μ)), ab])(∂_μ)d , FieldNumber = 5

LocalField = d[t,x,y,z];
LocalFieldNumber = 5;

AlgebraDDResult[tt1,tt2][LocalFieldNumber]=(Collect[(AlgebraDD[tt1,tt2][LocalField] ==Closure[LocalField])/.LagrangianParameters/.SetW,FieldListWithDerivatives,Simplify])

Clear[LocalField,LocalFieldNumber]

{(Subscript[(OverHat[D]), a]^i),(Subscript[(OverHat[D]), b]^j)}(Subscript[A, μ])= 2 i (Subscript[(δ^ij)((γ^ν)), ab])(∂_ν)(Subscript[A, μ]) -(∂_μ)(Subscript[(OverHat[(OverHat[λ])]), ab]^ij)  , 
FieldNumber = 7
(Subscript[(OverHat[(OverHat[λ])]), ab]^ij)= 2 i (δ^ij) (Subscript[((γ^ν)), ab])(Subscript[A, ν]) - 2 i ((((σ^2)))^ij)( i (Subscript[C, ab])A - (Subscript[((γ^5)), ab])B)

LocalFieldNumber = 7;

λ[tt1,tt2][ii_,jj_][aa_,bb_][t,x,y,z] := 2 I  KroneckerDelta[ii,jj] Sum[γdown[nu][[aa,bb]]Am[nu][t,x,y,z],{nu,0,3}] - 2I sigma[2][[ii,jj]](I Cmetric[[aa,bb]] A[t,x,y,z] -γ5down[[aa,bb]]B[t,x,y,z])

λTable[tt1,tt2][t,x,y,z] = Table[λ[tt1,tt2][ii,jj][aa,bb][t,x,y,z],{ii,1,2},{aa,1,4},{jj,1,2},{bb,1,4}];

AlgebraDDResult[tt1,tt2][LocalFieldNumber]=Collect[(Table[AlgebraDD[tt1,tt2][Am[mu][t,x,y,z]],{mu,0,3}] ==Table[Closure[Am[mu][t,x,y,z]]- D[λTable[tt1,tt2][t,x,y,z],Ϛ[mu]],{mu,0,3}] )/.LagrangianParameters/.SetW,FieldListWithDerivatives,Simplify]

Clear[LocalFieldNumber]

{(Subscript[(OverHat[D]), a]^i),(Subscript[(OverHat[D]), b]^j)}(Subscript[(OverHat[Ψ]), c]^k)= 2 i (Subscript[(δ^ij)((γ^μ)), ab])(∂_μ)(Subscript[(OverHat[Ψ]), c]^k)  , FieldNumber = 16

LocalFieldNumber = 16;

AlgebraDDResult[tt1,tt2][LocalFieldNumber]=Collect[(Table[AlgebraDD[tt1,tt2][Ψh[kk][cc][t,x,y,z]],{kk,1,2},{cc,1,4}] ==Table[Closure[Ψh[kk][cc][t,x,y,z]],{kk,1,2},{cc,1,4}])/.LagrangianParameters/.SetW,FieldListWithDerivatives,Simplify]

tt2
tt2

Clear[LocalFieldNumber]

Closure summary of N=2 Vector hat-hat algebra

Table[AlgebraDDResult[tt1,tt2][FieldIndex],{FieldIndex,{1,2,3,4,5,7,16}}]==Table[True,{FieldIndex,{1,2,3,4,5,7,16}}]

tt1
tt2

4D, N=2 Chiral-Chiral Algebra, both hatted {(Subscript[(OverHat[D]), a]^i),(Subscript[(OverHat[D]), b]^j)},

, FieldNumber = 8

tt1
tt2

AlgebraDDResult[tt1,tt2][8]=(AlgebraDD[tt1,tt2][Ah[t,x,y,z]] ==Closure[Ah[t,x,y,z]] + Table[I sigma[2][[ii,jj]] 2 I Cmetric[[aa,bb]]Fb[t,x,y,z],{ii,1,2},{aa,1,4},{jj,1,2},{bb,1,4}])/.LagrangianParameters/.SetW//Simplify

, FieldNumber = 9

AlgebraDDResult[tt1,tt2][9]=(AlgebraDD[tt1,tt2][Ab[t,x,y,z]] ==Closure[Ab[t,x,y,z]] - Table[I sigma[2][[ii,jj]] 2 I Cmetric[[aa,bb]]Fh[t,x,y,z],{ii,1,2},{aa,1,4},{jj,1,2},{bb,1,4}])/.LagrangianParameters/.SetW//Simplify

, FieldNumber = 10

AlgebraDDResult[tt1,tt2][10]=(AlgebraDD[tt1,tt2][Bh[t,x,y,z]] ==Closure[Bh[t,x,y,z]] +  Table[I sigma[2][[ii,jj]] 2 I Cmetric[[aa,bb]]Gb[t,x,y,z],{ii,1,2},{aa,1,4},{jj,1,2},{bb,1,4}])/.LagrangianParameters/.SetW//Simplify

, FieldNumber = 11

AlgebraDDResult[tt1,tt2][11]=(AlgebraDD[tt1,tt2][Bb[t,x,y,z]] ==Closure[Bb[t,x,y,z]] -  Table[I sigma[2][[ii,jj]] 2 I Cmetric[[aa,bb]]Gh[t,x,y,z],{ii,1,2},{aa,1,4},{jj,1,2},{bb,1,4}])/.LagrangianParameters/.SetW//Simplify

, FieldNumber = 12

AlgebraDDResult[tt1,tt2][12]=(AlgebraDD[tt1,tt2][Fh[t,x,y,z]] ==Closure[Fh[t,x,y,z]] + Table[I sigma[2][[ii,jj]]2 I Cmetric[[aa,bb]]Lap[Ab[t,x,y,z]],{ii,1,2},{aa,1,4},{jj,1,2},{bb,1,4}])/.LagrangianParameters/.SetW//Simplify

, FieldNumber = 13

AlgebraDDResult[tt1,tt2][13]=(AlgebraDD[tt1,tt2][Fb[t,x,y,z]] ==Closure[Fb[t,x,y,z]] - Table[I sigma[2][[ii,jj]]2 I Cmetric[[aa,bb]]Lap[Ah[t,x,y,z]],{ii,1,2},{aa,1,4},{jj,1,2},{bb,1,4}])/.LagrangianParameters/.SetW//Simplify

, FieldNumber = 14

AlgebraDDResult[tt1,tt2][14]=(AlgebraDD[tt1,tt2][Gh[t,x,y,z]] ==Closure[Gh[t,x,y,z]] + Table[I sigma[2][[ii,jj]]2 I Cmetric[[aa,bb]]Lap[Bb[t,x,y,z]],{ii,1,2},{aa,1,4},{jj,1,2},{bb,1,4}])/.LagrangianParameters/.SetW//Simplify

, FieldNumber = 15

AlgebraDDResult[tt1,tt2][15]=(AlgebraDD[tt1,tt2][Gb[t,x,y,z]] ==Closure[Gb[t,x,y,z]] - Table[I sigma[2][[ii,jj]]2 I Cmetric[[aa,bb]]Lap[Bh[t,x,y,z]],{ii,1,2},{aa,1,4},{jj,1,2},{bb,1,4}])/.LagrangianParameters/.SetW//Simplify

, FieldNumber = 6

AlgebraDDResult[tt1,tt2][6]=(Table[AlgebraDD[tt1,tt2][Ψ[kk][cc][t,x,y,z]],{kk,1,2},{cc,1,4}] ==Table[Closure[Ψ[kk][cc][t,x,y,z]],{kk,1,2},{cc,1,4}] - Table[sigma[2][[ii,jj]]Sum[sigma[2][[kk,rr]]2I Cmetric[[aa,bb]]γ[mu][[cc,dd]]D[Ψ[rr][dd][t,x,y,z],Ϛ[mu]],{mu,0,3},{dd,1,4},{rr,1,2}],{kk,1,2},{cc,1,4},{ii,1,2},{aa,1,4},{jj,1,2},{bb,1,4}])/.LagrangianParameters/.SetW//Simplify

Closure summary of N=2 chiral-chiral hat-hat algebra

Table[AlgebraDDResult[tt1,tt2][FieldIndex],{FieldIndex,{8,9,10,11,12,13,14,15,6}}]==Table[True,{FieldIndex,{8,9,10,11,12,13,14,15,6}}]

tt1
tt2

Closure summary of N=2 chiral-vector hat-hat algebra (tt1=tt2=0)

Table[AlgebraDDResult[tt1,tt2][FieldIndex],{FieldIndex,1,16}]==Table[True,{FieldIndex,1,16}]

tt1
tt2

Algebra {(Subscript[D, a]^i),(Subscript[(OverHat[D]), b]^j)} MUST have central charges as proved in “4D N=4 CV 181120 No Closure.nb”

tt1 = 1, tt2 = 0, Zgen[lfn,fn][i,j] := ((Z[lfn,fn])^ij)  for central charge coding ease

Clear[tt1,tt2,ss]
tt1 = 1
tt2 = 0

ZGen[lfn,fn][i_,j_]:= z[lfn,fn,0] (δ^ij) + z[lfn,fn,1] ((((σ^1)))^ij)+ i z[lfn,fn,2] ((((σ^2)))^ij)+z[lfn,fn,3] ((((σ^3)))^ij)

ZGen[lfn_,fn_][ii_,jj_] := z[lfn,fn,0] KroneckerDelta[ii,jj] + z[lfn,fn,1] sigma[1][[ii,jj]] + I z[lfn,fn,2] sigma[2][[ii,jj]]+z[lfn,fn,3] sigma[3][[ii,jj]]

{(Subscript[D, a]^i),(Subscript[(OverHat[D]), b]^j)}A= 2i ((((σ^3)))^ij)(Subscript[((γ^μ)), ab])(∂_μ)(OverHat[A])  +2i ((((σ^1)))^ij)(Subscript[((γ^μ)), ab])(∂_μ)(Overscript[A, ⌣]) -2 i ((((σ^2)))^ij)(Subscript[((γ^5)(γ^μ)), ab])(∂_μ)(Overscript[B, ⌣]) +2 i ((((σ^2)))^ij)(Subscript[((γ^5)), ab])(Overscript[G, ⌣])  
FieldNumber = 1

Set Local Field

LocalField = A[t,x,y,z];
LocalFieldNumber = 1;
LocalNumZCharges =8;
CentralChargeList[LocalFieldNumber] = Table[z[LocalFieldNumber,FieldNumber,sigmaNumber],{FieldNumber,1,LocalNumZCharges},{sigmaNumber,0,3}]//Flatten

Solving for the Central Charges

ZAnsatzTerms[LocalFieldNumber][ii_,jj_][aa_,bb_]:=2 I ZGen[LocalFieldNumber,1][ii,jj]Sum[γdown[mu][[aa,bb]]D[Ah[t,x,y,z],Ϛ[mu]],{mu,0,3}]+2 ZGen[LocalFieldNumber,2][ii,jj]Sum[γ5γdown[mu][[aa,bb]]D[Bh[t,x,y,z],Ϛ[mu]],{mu,0,3}] + 2 I ZGen[LocalFieldNumber,3][ii,jj]Cmetric[[aa,bb]]Fh[t,x,y,z]+ 2 ZGen[LocalFieldNumber,4][ii,jj]γ5down[[aa,bb]]Gh[t,x,y,z]+2 I ZGen[LocalFieldNumber,5][ii,jj]Sum[γdown[mu][[aa,bb]]D[Ab[t,x,y,z],Ϛ[mu]],{mu,0,3}]+2 ZGen[LocalFieldNumber,6][ii,jj]Sum[γ5γdown[mu][[aa,bb]]D[Bb[t,x,y,z],Ϛ[mu]],{mu,0,3}] + 2 I ZGen[LocalFieldNumber,7][ii,jj]Cmetric[[aa,bb]]Fb[t,x,y,z]+ 2 ZGen[LocalFieldNumber,8][ii,jj]γ5down[[aa,bb]]Gb[t,x,y,z]

ZAnsatz[LocalFieldNumber] = Table[ZAnsatzTerms[LocalFieldNumber][ii,jj][aa,bb],{ii,1,2},{aa,1,4},{jj,1,2},{bb,1,4}];

AlgebraDDMinusClosure[tt1,tt2][LocalFieldNumber]=(Collect[(AlgebraDD[tt1,tt2][LocalField] - ZAnsatz[LocalFieldNumber])/.LagrangianParameters/.SetW,FieldListWithDerivatives,Simplify]);

AlgebraDDEqs[tt1,tt2][LocalFieldNumber] = DeleteCases[Table[(Coefficient[AlgebraDDMinusClosure[tt1,tt2][LocalFieldNumber][[ii,aa,jj,bb]],FieldListWithDerivatives[[FieldListIndex]]])==0,{ii,1,2},{aa,1,4},{jj,1,2},{bb,1,4},{FieldListIndex,1,Length[FieldListWithDerivatives]}]//Simplify//Flatten,True]//DeleteDuplicates;

Length[AlgebraDDEqs[tt1,tt2][LocalFieldNumber]]

ZSolutions[LocalFieldNumber] = Solve[AlgebraDDEqs[tt1,tt2][LocalFieldNumber],CentralChargeList[LocalFieldNumber]]//Flatten;

Final Results

DeleteCases[Table[Z[LocalFieldNumber,nn]==PrintSigmaProduct[Table[ZGen[LocalFieldNumber,nn][ii,jj],{ii,1,2},{jj,1,2}]/.ZSolutions[LocalFieldNumber]],{nn,1,LocalNumZCharges}],Z[_,_]==0]

ZTerms[LocalFieldNumber][ii_,jj_][aa_,bb_]:=2 I sigma[3][[ii,jj]]Sum[γdown[mu][[aa,bb]]D[Ah[t,x,y,z],Ϛ[mu]],{mu,0,3}]+2 I sigma[1][[ii,jj]]Sum[γdown[mu][[aa,bb]]D[Ab[t,x,y,z],Ϛ[mu]],{mu,0,3}]-2 I sigma[2][[ii,jj]]Sum[γ5γdown[mu][[aa,bb]]D[Bb[t,x,y,z],Ϛ[mu]],{mu,0,3}] + 2 I sigma[2][[ii,jj]]γ5down[[aa,bb]]Gb[t,x,y,z]

Z[LocalFieldNumber] = Table[ZTerms[LocalFieldNumber][ii,jj][aa,bb],{ii,1,2},{aa,1,4},{jj,1,2},{bb,1,4}];

AlgebraDDResult[tt1,tt2][LocalFieldNumber] =(Collect[(AlgebraDD[tt1,tt2][LocalField] == Z[LocalFieldNumber])/.LagrangianParameters/.SetW,FieldListWithDerivatives,Simplify])

Clear[LocalField,LocalFieldNumber,LocalNumZCharges]

{(Subscript[D, a]^i),(Subscript[(OverHat[D]), b]^j)}B= 2i (Subscript[(δ^ij)((γ^μ)), ab])(∂_μ)(OverHat[B])  -2 ((((σ^2)))^ij)(Subscript[C, ab])(Overscript[G, ⌣])  
FieldNumber = 2

Set Local Field

LocalField = B[t,x,y,z];
LocalFieldNumber = 2;
LocalNumZCharges =8;
CentralChargeList[LocalFieldNumber] = Table[z[LocalFieldNumber,FieldNumber,sigmaNumber],{FieldNumber,1,LocalNumZCharges},{sigmaNumber,0,3}]//Flatten

Solving for the Central Charges

ZAnsatzTerms[LocalFieldNumber][ii_,jj_][aa_,bb_]:=2 ZGen[LocalFieldNumber,1][ii,jj]Sum[γ5γdown[mu][[aa,bb]]D[Ah[t,x,y,z],Ϛ[mu]],{mu,0,3}]+2 I ZGen[LocalFieldNumber,2][ii,jj]Sum[γdown[mu][[aa,bb]]D[Bh[t,x,y,z],Ϛ[mu]],{mu,0,3}] + 2 ZGen[LocalFieldNumber,3][ii,jj]γ5down[[aa,bb]]Fh[t,x,y,z]+ 2I ZGen[LocalFieldNumber,4][ii,jj]Cmetric[[aa,bb]]Gh[t,x,y,z]+2 ZGen[LocalFieldNumber,5][ii,jj]Sum[γ5γdown[mu][[aa,bb]]D[Ab[t,x,y,z],Ϛ[mu]],{mu,0,3}]+2 I ZGen[LocalFieldNumber,6][ii,jj]Sum[γdown[mu][[aa,bb]]D[Bb[t,x,y,z],Ϛ[mu]],{mu,0,3}] + 2 ZGen[LocalFieldNumber,7][ii,jj]γ5down[[aa,bb]]Fb[t,x,y,z]+ 2I ZGen[LocalFieldNumber,8][ii,jj]Cmetric[[aa,bb]]Gb[t,x,y,z]

ZAnsatz[LocalFieldNumber] = Table[ZAnsatzTerms[LocalFieldNumber][ii,jj][aa,bb],{ii,1,2},{aa,1,4},{jj,1,2},{bb,1,4}];

AlgebraDDMinusClosure[tt1,tt2][LocalFieldNumber]=(Collect[(AlgebraDD[tt1,tt2][LocalField] - ZAnsatz[LocalFieldNumber])/.LagrangianParameters/.SetW,FieldListWithDerivatives,Simplify]);

AlgebraDDEqs[tt1,tt2][LocalFieldNumber] = DeleteCases[Table[(Coefficient[AlgebraDDMinusClosure[tt1,tt2][LocalFieldNumber][[ii,aa,jj,bb]],FieldListWithDerivatives[[FieldListIndex]]])==0,{ii,1,2},{aa,1,4},{jj,1,2},{bb,1,4},{FieldListIndex,1,Length[FieldListWithDerivatives]}]//Simplify//Flatten,True]//DeleteDuplicates;

Length[AlgebraDDEqs[tt1,tt2][LocalFieldNumber]]

ZSolutions[LocalFieldNumber] = Solve[AlgebraDDEqs[tt1,tt2][LocalFieldNumber],CentralChargeList[LocalFieldNumber]]//Flatten;

Final Results

DeleteCases[Table[Z[LocalFieldNumber,nn]==PrintSigmaProduct[Table[ZGen[LocalFieldNumber,nn][ii,jj],{ii,1,2},{jj,1,2}]/.ZSolutions[LocalFieldNumber]],{nn,1,LocalNumZCharges}],Z[_,_]==0]

ZTerms[LocalFieldNumber][ii_,jj_][aa_,bb_]:=2 I KroneckerDelta[ii,jj]Sum[γdown[mu][[aa,bb]]D[Bh[t,x,y,z],Ϛ[mu]],{mu,0,3}] - 2sigma[2][[ii,jj]]Cmetric[[aa,bb]]Gb[t,x,y,z]

Z[LocalFieldNumber] = Table[ZTerms[LocalFieldNumber][ii,jj][aa,bb],{ii,1,2},{aa,1,4},{jj,1,2},{bb,1,4}];

AlgebraDDResult[tt1,tt2][LocalFieldNumber] =(Collect[(AlgebraDD[tt1,tt2][LocalField] == Z[LocalFieldNumber])/.LagrangianParameters/.SetW,FieldListWithDerivatives,Simplify])

Clear[LocalField,LocalFieldNumber,LocalNumZCharges]

{(Subscript[D, a]^i),(Subscript[(OverHat[D]), b]^j)}F= 2i ((((σ^3)))^ij)(Subscript[((γ^μ)), ab])(∂_μ)(OverHat[F])  +2i ((((σ^2)))^ij)(Subscript[((γ^5)), ab])□(Overscript[B, ⌣]) +2 i ((((σ^1)))^ij)(Subscript[((γ^μ)), ab])(∂_μ)(Overscript[F, ⌣]) -2 i ((((σ^2)))^ij)(Subscript[((γ^5)(γ^μ)), ab])(∂_μ)(Overscript[G, ⌣])  
FieldNumber = 3

Set Local Field

LocalField = F[t,x,y,z];
LocalFieldNumber = 3;
LocalNumZCharges =8;
CentralChargeList[LocalFieldNumber] = Table[z[LocalFieldNumber,FieldNumber,sigmaNumber],{FieldNumber,1,LocalNumZCharges},{sigmaNumber,0,3}]//Flatten

Solving for the Central Charges

ZAnsatzTerms[LocalFieldNumber][ii_,jj_][aa_,bb_]:=2 I ZGen[LocalFieldNumber,1][ii,jj]Cmetric[[aa,bb]]Lap[Ah[t,x,y,z]]+2 ZGen[LocalFieldNumber,2][ii,jj]γ5down[[aa,bb]]Lap[Bh[t,x,y,z]]+ 2 I ZGen[LocalFieldNumber,3][ii,jj]Sum[γdown[mu][[aa,bb]]D[Fh[t,x,y,z],Ϛ[mu]],{mu,0,3}]+2 ZGen[LocalFieldNumber,4][ii,jj]Sum[γ5γdown[mu][[aa,bb]]D[Gh[t,x,y,z],Ϛ[mu]],{mu,0,3}]+2 I ZGen[LocalFieldNumber,5][ii,jj]Cmetric[[aa,bb]]Lap[Ab[t,x,y,z]]+2 ZGen[LocalFieldNumber,6][ii,jj]γ5down[[aa,bb]]Lap[Bb[t,x,y,z]]+ 2 I ZGen[LocalFieldNumber,7][ii,jj]Sum[γdown[mu][[aa,bb]]D[Fb[t,x,y,z],Ϛ[mu]],{mu,0,3}]+2 ZGen[LocalFieldNumber,8][ii,jj]Sum[γ5γdown[mu][[aa,bb]]D[Gb[t,x,y,z],Ϛ[mu]],{mu,0,3}]

ZAnsatz[LocalFieldNumber] = Table[ZAnsatzTerms[LocalFieldNumber][ii,jj][aa,bb],{ii,1,2},{aa,1,4},{jj,1,2},{bb,1,4}];

AlgebraDDMinusClosure[tt1,tt2][LocalFieldNumber]=(Collect[(AlgebraDD[tt1,tt2][LocalField] - ZAnsatz[LocalFieldNumber])/.LagrangianParameters/.SetW,FieldListWithDerivatives,Simplify]);

AlgebraDDEqs[tt1,tt2][LocalFieldNumber] = DeleteCases[Table[(Coefficient[AlgebraDDMinusClosure[tt1,tt2][LocalFieldNumber][[ii,aa,jj,bb]],FieldListWithDerivatives[[FieldListIndex]]])==0,{ii,1,2},{aa,1,4},{jj,1,2},{bb,1,4},{FieldListIndex,1,Length[FieldListWithDerivatives]}]//Simplify//Flatten,True]//DeleteDuplicates;

Length[AlgebraDDEqs[tt1,tt2][LocalFieldNumber]]

ZSolutions[LocalFieldNumber] = Solve[AlgebraDDEqs[tt1,tt2][LocalFieldNumber],CentralChargeList[LocalFieldNumber]]//Flatten;

Final Results

DeleteCases[Table[Z[LocalFieldNumber,nn]==PrintSigmaProduct[Table[ZGen[LocalFieldNumber,nn][ii,jj],{ii,1,2},{jj,1,2}]/.ZSolutions[LocalFieldNumber]],{nn,1,LocalNumZCharges}],Z[_,_]==0]

ZTerms[LocalFieldNumber][ii_,jj_][aa_,bb_]:=2 I sigma[3][[ii,jj]]Sum[γdown[mu][[aa,bb]]D[Fh[t,x,y,z],Ϛ[mu]],{mu,0,3}]+2 I sigma[2][[ii,jj]]γ5down[[aa,bb]]Lap[Bb[t,x,y,z]]+ 2 I sigma[1][[ii,jj]]Sum[γdown[mu][[aa,bb]]D[Fb[t,x,y,z],Ϛ[mu]],{mu,0,3}]-2 I sigma[2][[ii,jj]]Sum[γ5γdown[mu][[aa,bb]]D[Gb[t,x,y,z],Ϛ[mu]],{mu,0,3}]

Z[LocalFieldNumber] = Table[ZTerms[LocalFieldNumber][ii,jj][aa,bb],{ii,1,2},{aa,1,4},{jj,1,2},{bb,1,4}];

AlgebraDDResult[tt1,tt2][LocalFieldNumber] =(Collect[(AlgebraDD[tt1,tt2][LocalField] == Z[LocalFieldNumber])/.LagrangianParameters/.SetW,FieldListWithDerivatives,Simplify])

Clear[LocalField,LocalFieldNumber,LocalNumZCharges]

{(Subscript[D, a]^i),(Subscript[(OverHat[D]), b]^j)}G= 2i ((((σ^3)))^ij)(Subscript[((γ^μ)), ab])(∂_μ)(OverHat[G])  -2i ((((σ^2)))^ij)(Subscript[((γ^5)), ab])□(Overscript[A, ⌣]) -2 i ((((σ^2)))^ij)(Subscript[((γ^5)(γ^μ)), ab])(∂_μ)(Overscript[F, ⌣]) -2 i ((((σ^1)))^ij)(Subscript[((γ^μ)), ab])(∂_μ)(Overscript[G, ⌣])  
FieldNumber = 4

Set Local Field

LocalField = G[t,x,y,z];
LocalFieldNumber = 4;
LocalNumZCharges =8;
CentralChargeList[LocalFieldNumber] = Table[z[LocalFieldNumber,FieldNumber,sigmaNumber],{FieldNumber,1,LocalNumZCharges},{sigmaNumber,0,3}]//Flatten

Solving for the Central Charges

ZAnsatzTerms[LocalFieldNumber][ii_,jj_][aa_,bb_]:=2 ZGen[LocalFieldNumber,1][ii,jj]γ5down[[aa,bb]]Lap[Ah[t,x,y,z]]+2I ZGen[LocalFieldNumber,2][ii,jj]Cmetric[[aa,bb]]Lap[Bh[t,x,y,z]]+ 2 ZGen[LocalFieldNumber,3][ii,jj]Sum[γ5γdown[mu][[aa,bb]]D[Fh[t,x,y,z],Ϛ[mu]],{mu,0,3}]+2 I ZGen[LocalFieldNumber,4][ii,jj]Sum[γdown[mu][[aa,bb]]D[Gh[t,x,y,z],Ϛ[mu]],{mu,0,3}]+2 ZGen[LocalFieldNumber,5][ii,jj]γ5down[[aa,bb]]Lap[Ab[t,x,y,z]]+2I ZGen[LocalFieldNumber,6][ii,jj]Cmetric[[aa,bb]]Lap[Bb[t,x,y,z]]+ 2 ZGen[LocalFieldNumber,7][ii,jj]Sum[γ5γdown[mu][[aa,bb]]D[Fb[t,x,y,z],Ϛ[mu]],{mu,0,3}]+2 I ZGen[LocalFieldNumber,8][ii,jj]Sum[γdown[mu][[aa,bb]]D[Gb[t,x,y,z],Ϛ[mu]],{mu,0,3}]

ZAnsatz[LocalFieldNumber] = Table[ZAnsatzTerms[LocalFieldNumber][ii,jj][aa,bb],{ii,1,2},{aa,1,4},{jj,1,2},{bb,1,4}];

AlgebraDDMinusClosure[tt1,tt2][LocalFieldNumber]=(Collect[(AlgebraDD[tt1,tt2][LocalField] - ZAnsatz[LocalFieldNumber])/.LagrangianParameters/.SetW,FieldListWithDerivatives,Simplify]);

AlgebraDDEqs[tt1,tt2][LocalFieldNumber] = DeleteCases[Table[(Coefficient[AlgebraDDMinusClosure[tt1,tt2][LocalFieldNumber][[ii,aa,jj,bb]],FieldListWithDerivatives[[FieldListIndex]]])==0,{ii,1,2},{aa,1,4},{jj,1,2},{bb,1,4},{FieldListIndex,1,Length[FieldListWithDerivatives]}]//Simplify//Flatten,True]//DeleteDuplicates;

Length[AlgebraDDEqs[tt1,tt2][LocalFieldNumber]]

ZSolutions[LocalFieldNumber] = Solve[AlgebraDDEqs[tt1,tt2][LocalFieldNumber],CentralChargeList[LocalFieldNumber]]//Flatten;

Final Results

DeleteCases[Table[Z[LocalFieldNumber,nn]==PrintSigmaProduct[Table[ZGen[LocalFieldNumber,nn][ii,jj],{ii,1,2},{jj,1,2}]/.ZSolutions[LocalFieldNumber]],{nn,1,LocalNumZCharges}],Z[_,_]==0]

ZTerms[LocalFieldNumber][ii_,jj_][aa_,bb_]:=2 I sigma[3][[ii,jj]]Sum[γdown[mu][[aa,bb]]D[Gh[t,x,y,z],Ϛ[mu]],{mu,0,3}]-2I sigma[2][[ii,jj]]γ5down[[aa,bb]]Lap[Ab[t,x,y,z]]-2 I sigma[2][[ii,jj]]Sum[γ5γdown[mu][[aa,bb]]D[Fb[t,x,y,z],Ϛ[mu]],{mu,0,3}]-2 I sigma[1][[ii,jj]]Sum[γdown[mu][[aa,bb]]D[Gb[t,x,y,z],Ϛ[mu]],{mu,0,3}]

Z[LocalFieldNumber] = Table[ZTerms[LocalFieldNumber][ii,jj][aa,bb],{ii,1,2},{aa,1,4},{jj,1,2},{bb,1,4}];

AlgebraDDResult[tt1,tt2][LocalFieldNumber] =(Collect[(AlgebraDD[tt1,tt2][LocalField] == Z[LocalFieldNumber])/.LagrangianParameters/.SetW,FieldListWithDerivatives,Simplify])

Clear[LocalField,LocalFieldNumber,LocalNumZCharges]

{(Subscript[D, a]^i),(Subscript[(OverHat[D]), b]^j)}d= 2i ((((σ^2)))^ij)(Subscript[((γ^5)), ab])□(OverHat[A])+ 2i ((((σ^2)))^ij)(Subscript[((γ^5)(γ^μ)), ab])(∂_μ)(OverHat[F])  +2 i ((((σ^1)))^ij)(Subscript[((γ^μ)), ab])(∂_μ)(OverHat[G]) +2 i ((((σ^3)))^ij)(Subscript[((γ^μ)), ab])(∂_μ)(Overscript[G, ⌣])   
FieldNumber = 5

Set Local Field

LocalField = d[t,x,y,z];
LocalFieldNumber = 5;
LocalNumZCharges =8;
CentralChargeList[LocalFieldNumber] = Table[z[LocalFieldNumber,FieldNumber,sigmaNumber],{FieldNumber,1,LocalNumZCharges},{sigmaNumber,0,3}]//Flatten

Solving for the Central Charges

ZAnsatzTerms[LocalFieldNumber][ii_,jj_][aa_,bb_]:=2 ZGen[LocalFieldNumber,1][ii,jj]γ5down[[aa,bb]]Lap[Ah[t,x,y,z]]+2I ZGen[LocalFieldNumber,2][ii,jj]Cmetric[[aa,bb]]Lap[Bh[t,x,y,z]]+ 2 ZGen[LocalFieldNumber,3][ii,jj]Sum[γ5γdown[mu][[aa,bb]]D[Fh[t,x,y,z],Ϛ[mu]],{mu,0,3}]+2 I ZGen[LocalFieldNumber,4][ii,jj]Sum[γdown[mu][[aa,bb]]D[Gh[t,x,y,z],Ϛ[mu]],{mu,0,3}]+2 ZGen[LocalFieldNumber,5][ii,jj]γ5down[[aa,bb]]Lap[Ab[t,x,y,z]]+2I ZGen[LocalFieldNumber,6][ii,jj]Cmetric[[aa,bb]]Lap[Bb[t,x,y,z]]+ 2 ZGen[LocalFieldNumber,7][ii,jj]Sum[γ5γdown[mu][[aa,bb]]D[Fb[t,x,y,z],Ϛ[mu]],{mu,0,3}]+2 I ZGen[LocalFieldNumber,8][ii,jj]Sum[γdown[mu][[aa,bb]]D[Gb[t,x,y,z],Ϛ[mu]],{mu,0,3}]

ZAnsatz[LocalFieldNumber] = Table[ZAnsatzTerms[LocalFieldNumber][ii,jj][aa,bb],{ii,1,2},{aa,1,4},{jj,1,2},{bb,1,4}];

AlgebraDDMinusClosure[tt1,tt2][LocalFieldNumber]=(Collect[(AlgebraDD[tt1,tt2][LocalField] - ZAnsatz[LocalFieldNumber])/.LagrangianParameters/.SetW,FieldListWithDerivatives,Simplify]);

AlgebraDDEqs[tt1,tt2][LocalFieldNumber] = DeleteCases[Table[(Coefficient[AlgebraDDMinusClosure[tt1,tt2][LocalFieldNumber][[ii,aa,jj,bb]],FieldListWithDerivatives[[FieldListIndex]]])==0,{ii,1,2},{aa,1,4},{jj,1,2},{bb,1,4},{FieldListIndex,1,Length[FieldListWithDerivatives]}]//Simplify//Flatten,True]//DeleteDuplicates;

Length[AlgebraDDEqs[tt1,tt2][LocalFieldNumber]]

ZSolutions[LocalFieldNumber] = Solve[AlgebraDDEqs[tt1,tt2][LocalFieldNumber],CentralChargeList[LocalFieldNumber]]//Flatten;

Final Results

DeleteCases[Table[Z[LocalFieldNumber,nn]==PrintSigmaProduct[Table[ZGen[LocalFieldNumber,nn][ii,jj],{ii,1,2},{jj,1,2}]/.ZSolutions[LocalFieldNumber]],{nn,1,LocalNumZCharges}],Z[_,_]==0]

ZTerms[LocalFieldNumber][ii_,jj_][aa_,bb_]:=2 I sigma[2][[ii,jj]]γ5down[[aa,bb]]Lap[Ah[t,x,y,z]]+ 2 I sigma[2][[ii,jj]]Sum[γ5γdown[mu][[aa,bb]]D[Fh[t,x,y,z],Ϛ[mu]],{mu,0,3}]+2 I sigma[1][[ii,jj]]Sum[γdown[mu][[aa,bb]]D[Gh[t,x,y,z],Ϛ[mu]],{mu,0,3}]+2 I sigma[3][[ii,jj]]Sum[γdown[mu][[aa,bb]]D[Gb[t,x,y,z],Ϛ[mu]],{mu,0,3}]

Z[LocalFieldNumber] = Table[ZTerms[LocalFieldNumber][ii,jj][aa,bb],{ii,1,2},{aa,1,4},{jj,1,2},{bb,1,4}];

AlgebraDDResult[tt1,tt2][LocalFieldNumber] =(Collect[(AlgebraDD[tt1,tt2][LocalField] == Z[LocalFieldNumber])/.LagrangianParameters/.SetW,FieldListWithDerivatives,Simplify])

Clear[LocalField,LocalFieldNumber,LocalNumZCharges]

{(Subscript[D, a]^i),(Subscript[(OverHat[D]), b]^j)}(Subscript[A, μ])=-i ((((σ^1)))^ij)(Subscript[([(Subscript[γ, μ]),(Subscript[γ, ν])]), ab])(∂^ν)(OverHat[A])  -2 i ((((σ^1)))^ij)(Subscript[((Subscript[γ, μ])), ab])(OverHat[F]) +2 i ((((σ^2)))^ij)(Subscript[((γ^5)(Subscript[γ, μ])), ab])(OverHat[G]) +i ((((σ^3)))^ij)(Subscript[([(Subscript[γ, μ]),(Subscript[γ, ν])]), ab])(∂^ν)(Overscript[A, ⌣])- (δ^ij)(Subscript[((γ^5)[(Subscript[γ, μ]),(Subscript[γ, ν])]), ab])(∂^ν)(Overscript[B, ⌣])+2 i ((((σ^3)))^ij)(Subscript[((Subscript[γ, μ])), ab])(Overscript[F, ⌣])-(∂_μ)(Subscript[(OverHat[λ]), ab]^ij)  ,
FieldNumber = 7
(Subscript[(OverHat[λ]), ab]^ij)=  2 i ((((σ^2)))^ij) (Subscript[((γ^5)), ab])(OverHat[B])

Set Local Field

LocalFieldNumber = 7;
LocalNumZCharges =12;
CentralChargeList[LocalFieldNumber] = Table[z[LocalFieldNumber,FieldNumber,sigmaNumber],{FieldNumber,1,LocalNumZCharges},{sigmaNumber,0,3}]//Flatten

Solving for the Central Charges

ZAnsatzTerms[LocalFieldNumber][mu_][ii_,jj_][aa_,bb_]:=2 I ZGen[LocalFieldNumber,1][ii,jj]Sum[CommuteGammastdowndown[mu,nu][[aa,bb]]UD[Ah[t,x,y,z],Ϛ[nu]],{nu,0,3}] + 2 ZGen[LocalFieldNumber,2][ii,jj]Sum[(γ5.CommuteGammastdowndown[mu,nu])[[aa,bb]]UD[Bh[t,x,y,z],Ϛ[nu]],{nu,0,3}]  +2 I ZGen[LocalFieldNumber,3][ii,jj]γstdowndown[mu][[aa,bb]]Fh[t,x,y,z]+ 2 ZGen[LocalFieldNumber,4][ii,jj]γ5γstdowndown[mu][[aa,bb]]Gh[t,x,y,z] + 2 I ZGen[LocalFieldNumber,5][ii,jj]Sum[CommuteGammastdowndown[mu,nu][[aa,bb]]UD[Ab[t,x,y,z],Ϛ[nu]],{nu,0,3}] + 2 ZGen[LocalFieldNumber,6][ii,jj]Sum[(γ5.CommuteGammastdowndown[mu,nu])[[aa,bb]]UD[Bb[t,x,y,z],Ϛ[nu]],{nu,0,3}]+2 I ZGen[LocalFieldNumber,7][ii,jj]γstdowndown[mu][[aa,bb]]Fb[t,x,y,z]+ 2 ZGen[LocalFieldNumber,8][ii,jj]γ5γstdowndown[mu][[aa,bb]]Gb[t,x,y,z]

Do[ZAnsatz[LocalFieldNumber][mu] = Table[ZAnsatzTerms[LocalFieldNumber][mu][ii,jj][aa,bb],{ii,1,2},{aa,1,4},{jj,1,2},{bb,1,4}],{mu,0,3}];

λAnsatzTerms[tt1,tt2][ii_,jj_][aa_,bb_][t,x,y,z] :=  2I ZGen[LocalFieldNumber,9][ii,jj]Cmetric[[aa,bb]] Ah[t,x,y,z]  + 2 ZGen[LocalFieldNumber,10][ii,jj]γ5down[[aa,bb]] Bh[t,x,y,z]+2I ZGen[LocalFieldNumber,11][ii,jj]Cmetric[[aa,bb]] Ab[t,x,y,z]+ 2 ZGen[LocalFieldNumber,12][ii,jj]γ5down[[aa,bb]] Bb[t,x,y,z]

λAnsatz[tt1,tt2][t,x,y,z] = Table[λAnsatzTerms[tt1,tt2][ii,jj][aa,bb][t,x,y,z],{ii,1,2},{aa,1,4},{jj,1,2},{bb,1,4}];

AlgebraDDMinusClosure[tt1,tt2][LocalFieldNumber]=Table[Collect[(AlgebraDD[tt1,tt2][Am[mu][t,x,y,z]] - (ZAnsatz[LocalFieldNumber][mu]-D[λAnsatz[tt1,tt2][t,x,y,z],Ϛ[mu]]))/.LagrangianParameters/.SetW,FieldListWithDerivatives,Simplify],{mu,0,3}];

Dimensions[AlgebraDDMinusClosure[tt1,tt2][LocalFieldNumber]]

AlgebraDDEqs[tt1,tt2][LocalFieldNumber] = DeleteCases[Table[(Coefficient[AlgebraDDMinusClosure[tt1,tt2][LocalFieldNumber][[mu,ii,aa,jj,bb]],FieldListWithDerivatives[[FieldListIndex]]])==0,{mu,1,4},{ii,1,2},{aa,1,4},{jj,1,2},{bb,1,4},{FieldListIndex,1,Length[FieldListWithDerivatives]}]//Simplify//Flatten,True]//DeleteDuplicates;

Length[AlgebraDDEqs[tt1,tt2][LocalFieldNumber]]

ZSolutions[LocalFieldNumber] = Solve[AlgebraDDEqs[tt1,tt2][LocalFieldNumber],CentralChargeList[LocalFieldNumber]]//Flatten

Final Results

DeleteCases[Table[Z[LocalFieldNumber,nn]==PrintSigmaProduct[Table[ZGen[LocalFieldNumber,nn][ii,jj],{ii,1,2},{jj,1,2}]/.ZSolutions[LocalFieldNumber]],{nn,1,LocalNumZCharges}],Z[_,_]==0]

ZTerms[LocalFieldNumber][mu_][ii_,jj_][aa_,bb_]:=-I sigma[1][[ii,jj]]Sum[CommuteGammastdowndown[mu,nu][[aa,bb]]UD[Ah[t,x,y,z],Ϛ[nu]],{nu,0,3}] + -2 I sigma[1][[ii,jj]]γstdowndown[mu][[aa,bb]]Fh[t,x,y,z]+ 2 I sigma[2][[ii,jj]]γ5γstdowndown[mu][[aa,bb]]Gh[t,x,y,z] + I sigma[3][[ii,jj]]Sum[CommuteGammastdowndown[mu,nu][[aa,bb]]UD[Ab[t,x,y,z],Ϛ[nu]],{nu,0,3}] + -KroneckerDelta[ii,jj]Sum[(γ5.CommuteGammastdowndown[mu,nu])[[aa,bb]]UD[Bb[t,x,y,z],Ϛ[nu]],{nu,0,3}]+2 I sigma[3][[ii,jj]]γstdowndown[mu][[aa,bb]]Fb[t,x,y,z]

Do[Z[LocalFieldNumber][mu] = Table[ZTerms[LocalFieldNumber][mu][ii,jj][aa,bb],{ii,1,2},{aa,1,4},{jj,1,2},{bb,1,4}],{mu,0,3}];

λTerms[LocalFieldNumber][ii_,jj_][aa_,bb_][t,x,y,z]:= 2 I sigma[2][[ii,jj]]γ5down[[aa,bb]] Bh[t,x,y,z]

λTable[tt1,tt2][t,x,y,z] = Table[λTerms[LocalFieldNumber][ii,jj][aa,bb][t,x,y,z],{ii,1,2},{aa,1,4},{jj,1,2},{bb,1,4}];

AlgebraDDResult[tt1,tt2][LocalFieldNumber] =Table[AlgebraDD[tt1,tt2][Am[mu][t,x,y,z]],{mu,0,3}] == Table[Z[LocalFieldNumber][mu]-D[λTable[tt1,tt2][t,x,y,z],Ϛ[mu]],{mu,0,3}]/.LagrangianParameters/.SetW//Simplify

Clear[LocalFieldNumber,LocalNumZCharges]

{(Subscript[D, a]^i),(Subscript[(OverHat[D]), b]^j)}(OverHat[A])= 2i ((((σ^3)))^ij)(Subscript[((γ^μ)), ab])(∂_μ)A  -2 i ((((σ^2)))^ij)(Subscript[((γ^5)), ab])d +i(1/2) ((((σ^1)))^ij) (Subscript[([(γ^μ),(γ^ν)]), ab])(Subscript[F, μν]) 
FieldNumber = 8

Set Local Field

LocalField = Ah[t,x,y,z];
LocalFieldNumber = 8;
LocalNumZCharges =7;
CentralChargeList[LocalFieldNumber] = Table[z[LocalFieldNumber,FieldNumber,sigmaNumber],{FieldNumber,1,LocalNumZCharges},{sigmaNumber,0,3}]//Flatten

Solving for the Central Charges

ZAnsatzTerms[LocalFieldNumber][ii_,jj_][aa_,bb_]:=2 I ZGen[LocalFieldNumber,1][ii,jj]Sum[γdown[mu][[aa,bb]]D[A[t,x,y,z],Ϛ[mu]],{mu,0,3}]+2 ZGen[LocalFieldNumber,2][ii,jj]Sum[γ5γdown[mu][[aa,bb]]D[B[t,x,y,z],Ϛ[mu]],{mu,0,3}] + 2 I ZGen[LocalFieldNumber,3][ii,jj]Cmetric[[aa,bb]]F[t,x,y,z]+ 2 ZGen[LocalFieldNumber,4][ii,jj]γ5down[[aa,bb]]G[t,x,y,z]+ 2 ZGen[LocalFieldNumber,5][ii,jj]γ5down[[aa,bb]]d[t,x,y,z]+
2I ZGen[LocalFieldNumber,6][ii,jj]Cmetric[[aa,bb]]Sum[UD[Am[mu][t,x,y,z],Ϛ[mu]],{mu,0,3}] + (1/2)I ZGen[LocalFieldNumber,7][ii,jj]Sum[CommuteGammadown[mu,nu][[aa,bb]]Fmn[mu,nu][t,x,y,z],{mu,0,3},{nu,0,3}]

ZAnsatz[LocalFieldNumber] = Table[ZAnsatzTerms[LocalFieldNumber][ii,jj][aa,bb],{ii,1,2},{aa,1,4},{jj,1,2},{bb,1,4}];

AlgebraDDMinusClosure[tt1,tt2][LocalFieldNumber]=(Collect[(AlgebraDD[tt1,tt2][LocalField] - ZAnsatz[LocalFieldNumber])/.LagrangianParameters/.SetW,FieldListWithDerivatives,Simplify]);

AlgebraDDEqs[tt1,tt2][LocalFieldNumber] = DeleteCases[Table[(Coefficient[AlgebraDDMinusClosure[tt1,tt2][LocalFieldNumber][[ii,aa,jj,bb]],FieldListWithDerivatives[[FieldListIndex]]])==0,{ii,1,2},{aa,1,4},{jj,1,2},{bb,1,4},{FieldListIndex,1,Length[FieldListWithDerivatives]}]//Simplify//Flatten,True]//DeleteDuplicates;

Length[AlgebraDDEqs[tt1,tt2][LocalFieldNumber]]

ZSolutions[LocalFieldNumber] = Solve[AlgebraDDEqs[tt1,tt2][LocalFieldNumber],CentralChargeList[LocalFieldNumber]]//Flatten;

Final Results

DeleteCases[Table[Z[LocalFieldNumber,nn]==PrintSigmaProduct[Table[ZGen[LocalFieldNumber,nn][ii,jj],{ii,1,2},{jj,1,2}]/.ZSolutions[LocalFieldNumber]],{nn,1,LocalNumZCharges}],Z[_,_]==0]

ZTerms[LocalFieldNumber][ii_,jj_][aa_,bb_]:=2 I sigma[3][[ii,jj]]Sum[γdown[mu][[aa,bb]]D[A[t,x,y,z],Ϛ[mu]],{mu,0,3}]-2I sigma[2][[ii,jj]]γ5down[[aa,bb]]d[t,x,y,z]+
(1/2)I sigma[1][[ii,jj]]Sum[CommuteGammadown[mu,nu][[aa,bb]]Fmn[mu,nu][t,x,y,z],{mu,0,3},{nu,0,3}]

Z[LocalFieldNumber] = Table[ZTerms[LocalFieldNumber][ii,jj][aa,bb],{ii,1,2},{aa,1,4},{jj,1,2},{bb,1,4}];

AlgebraDDResult[tt1,tt2][LocalFieldNumber] =(Collect[(AlgebraDD[tt1,tt2][LocalField] == Z[LocalFieldNumber])/.LagrangianParameters/.SetW,FieldListWithDerivatives,Simplify])

Clear[LocalField,LocalFieldNumber,LocalNumZCharges]

{(Subscript[D, a]^i),(Subscript[(OverHat[D]), b]^j)}(Overscript[A, ⌣])= 2i ((((σ^1)))^ij)(Subscript[((γ^μ)), ab])(∂_μ)A  +2 i ((((σ^2)))^ij)(Subscript[((γ^5)), ab])G -i(1/2) ((((σ^3)))^ij) (Subscript[([(γ^μ),(γ^ν)]), ab])(Subscript[F, μν]) 
FieldNumber = 9

Set Local Field

LocalField = Ab[t,x,y,z];
LocalFieldNumber = 9;
LocalNumZCharges =7;
CentralChargeList[LocalFieldNumber] = Table[z[LocalFieldNumber,FieldNumber,sigmaNumber],{FieldNumber,1,LocalNumZCharges},{sigmaNumber,0,3}]//Flatten

Solving for the Central Charges

ZAnsatzTerms[LocalFieldNumber][ii_,jj_][aa_,bb_]:=2 I ZGen[LocalFieldNumber,1][ii,jj]Sum[γdown[mu][[aa,bb]]D[A[t,x,y,z],Ϛ[mu]],{mu,0,3}]+2 ZGen[LocalFieldNumber,2][ii,jj]Sum[γ5γdown[mu][[aa,bb]]D[B[t,x,y,z],Ϛ[mu]],{mu,0,3}] + 2 I ZGen[LocalFieldNumber,3][ii,jj]Cmetric[[aa,bb]]F[t,x,y,z]+ 2 ZGen[LocalFieldNumber,4][ii,jj]γ5down[[aa,bb]]G[t,x,y,z]+ 2 ZGen[LocalFieldNumber,5][ii,jj]γ5down[[aa,bb]]d[t,x,y,z]+
2I ZGen[LocalFieldNumber,6][ii,jj]Cmetric[[aa,bb]]Sum[UD[Am[mu][t,x,y,z],Ϛ[mu]],{mu,0,3}] + (1/2)I ZGen[LocalFieldNumber,7][ii,jj]Sum[CommuteGammadown[mu,nu][[aa,bb]]Fmn[mu,nu][t,x,y,z],{mu,0,3},{nu,0,3}]

ZAnsatz[LocalFieldNumber] = Table[ZAnsatzTerms[LocalFieldNumber][ii,jj][aa,bb],{ii,1,2},{aa,1,4},{jj,1,2},{bb,1,4}];

AlgebraDDMinusClosure[tt1,tt2][LocalFieldNumber]=(Collect[(AlgebraDD[tt1,tt2][LocalField] - ZAnsatz[LocalFieldNumber])/.LagrangianParameters/.SetW,FieldListWithDerivatives,Simplify]);

AlgebraDDEqs[tt1,tt2][LocalFieldNumber] = DeleteCases[Table[(Coefficient[AlgebraDDMinusClosure[tt1,tt2][LocalFieldNumber][[ii,aa,jj,bb]],FieldListWithDerivatives[[FieldListIndex]]])==0,{ii,1,2},{aa,1,4},{jj,1,2},{bb,1,4},{FieldListIndex,1,Length[FieldListWithDerivatives]}]//Simplify//Flatten,True]//DeleteDuplicates;

Length[AlgebraDDEqs[tt1,tt2][LocalFieldNumber]]

ZSolutions[LocalFieldNumber] = Solve[AlgebraDDEqs[tt1,tt2][LocalFieldNumber],CentralChargeList[LocalFieldNumber]]//Flatten

Final Results

DeleteCases[Table[Z[LocalFieldNumber,nn]==PrintSigmaProduct[Table[ZGen[LocalFieldNumber,nn][ii,jj],{ii,1,2},{jj,1,2}]/.ZSolutions[LocalFieldNumber]],{nn,1,LocalNumZCharges}],Z[_,_]==0]

ZTerms[LocalFieldNumber][ii_,jj_][aa_,bb_]:=2 I sigma[1][[ii,jj]]Sum[γdown[mu][[aa,bb]]D[A[t,x,y,z],Ϛ[mu]],{mu,0,3}]+2 I sigma[2][[ii,jj]]γ5down[[aa,bb]]G[t,x,y,z]-(1/2)I sigma[3][[ii,jj]]Sum[CommuteGammadown[mu,nu][[aa,bb]]Fmn[mu,nu][t,x,y,z],{mu,0,3},{nu,0,3}]

Z[LocalFieldNumber] = Table[ZTerms[LocalFieldNumber][ii,jj][aa,bb],{ii,1,2},{aa,1,4},{jj,1,2},{bb,1,4}];

AlgebraDDResult[tt1,tt2][LocalFieldNumber] =(Collect[(AlgebraDD[tt1,tt2][LocalField] == Z[LocalFieldNumber])/.LagrangianParameters/.SetW,FieldListWithDerivatives,Simplify])

Clear[LocalField,LocalFieldNumber,LocalNumZCharges]

{(Subscript[D, a]^i),(Subscript[(OverHat[D]), b]^j)}(OverHat[B])= 2i (Subscript[(δ^ij)((γ^μ)), ab])(∂_μ)B  
FieldNumber = 10

Set Local Field

LocalField = Bh[t,x,y,z];
LocalFieldNumber = 10;
LocalNumZCharges =7;
CentralChargeList[LocalFieldNumber] = Table[z[LocalFieldNumber,FieldNumber,sigmaNumber],{FieldNumber,1,LocalNumZCharges},{sigmaNumber,0,3}]//Flatten

Solving for the Central Charges

ZAnsatzTerms[LocalFieldNumber][ii_,jj_][aa_,bb_]:=2 ZGen[LocalFieldNumber,1][ii,jj]Sum[γ5γdown[mu][[aa,bb]]D[A[t,x,y,z],Ϛ[mu]],{mu,0,3}]+2I ZGen[LocalFieldNumber,2][ii,jj]Sum[γdown[mu][[aa,bb]]D[B[t,x,y,z],Ϛ[mu]],{mu,0,3}] + 2 ZGen[LocalFieldNumber,3][ii,jj]γ5down[[aa,bb]]F[t,x,y,z]+ 2I ZGen[LocalFieldNumber,4][ii,jj]Cmetric[[aa,bb]]G[t,x,y,z]+ 2 I ZGen[LocalFieldNumber,5][ii,jj]Cmetric[[aa,bb]]d[t,x,y,z]+
2ZGen[LocalFieldNumber,6][ii,jj]γ5down[[aa,bb]]Sum[UD[Am[mu][t,x,y,z],Ϛ[mu]],{mu,0,3}] + (1/2) ZGen[LocalFieldNumber,7][ii,jj]Sum[(γ5.CommuteGammadown[mu,nu])[[aa,bb]]Fmn[mu,nu][t,x,y,z],{mu,0,3},{nu,0,3}]

ZAnsatz[LocalFieldNumber] = Table[ZAnsatzTerms[LocalFieldNumber][ii,jj][aa,bb],{ii,1,2},{aa,1,4},{jj,1,2},{bb,1,4}];

AlgebraDDMinusClosure[tt1,tt2][LocalFieldNumber]=(Collect[(AlgebraDD[tt1,tt2][LocalField] - ZAnsatz[LocalFieldNumber])/.LagrangianParameters/.SetW,FieldListWithDerivatives,Simplify]);

AlgebraDDEqs[tt1,tt2][LocalFieldNumber] = DeleteCases[Table[(Coefficient[AlgebraDDMinusClosure[tt1,tt2][LocalFieldNumber][[ii,aa,jj,bb]],FieldListWithDerivatives[[FieldListIndex]]])==0,{ii,1,2},{aa,1,4},{jj,1,2},{bb,1,4},{FieldListIndex,1,Length[FieldListWithDerivatives]}]//Simplify//Flatten,True]//DeleteDuplicates;

Length[AlgebraDDEqs[tt1,tt2][LocalFieldNumber]]

ZSolutions[LocalFieldNumber] = Solve[AlgebraDDEqs[tt1,tt2][LocalFieldNumber],CentralChargeList[LocalFieldNumber]]//Flatten

Final Results

DeleteCases[Table[Z[LocalFieldNumber,nn]==PrintSigmaProduct[Table[ZGen[LocalFieldNumber,nn][ii,jj],{ii,1,2},{jj,1,2}]/.ZSolutions[LocalFieldNumber]],{nn,1,LocalNumZCharges}],Z[_,_]==0]

ZTerms[LocalFieldNumber][ii_,jj_][aa_,bb_]:=2I KroneckerDelta[ii,jj]Sum[γdown[mu][[aa,bb]]D[B[t,x,y,z],Ϛ[mu]],{mu,0,3}]

Z[LocalFieldNumber] = Table[ZTerms[LocalFieldNumber][ii,jj][aa,bb],{ii,1,2},{aa,1,4},{jj,1,2},{bb,1,4}];

AlgebraDDResult[tt1,tt2][LocalFieldNumber] =(Collect[(AlgebraDD[tt1,tt2][LocalField] == Z[LocalFieldNumber])/.LagrangianParameters/.SetW,FieldListWithDerivatives,Simplify])

Clear[LocalField,LocalFieldNumber,LocalNumZCharges]

{(Subscript[D, a]^i),(Subscript[(OverHat[D]), b]^j)}(Overscript[B, ⌣])=- 2i ((((σ^2)))^ij)(Subscript[((γ^5)(γ^μ)), ab])(∂_μ)A -2 i ((((σ^2)))^ij)(Subscript[((γ^5)), ab])F+ (1/2)(δ^ij) (Subscript[((γ^5)[(γ^μ),(γ^ν)]), ab])(Subscript[F, μν]) 
FieldNumber = 11

Set Local Field

LocalField = Bb[t,x,y,z];
LocalFieldNumber = 11;
LocalNumZCharges =7;
CentralChargeList[LocalFieldNumber] = Table[z[LocalFieldNumber,FieldNumber,sigmaNumber],{FieldNumber,1,LocalNumZCharges},{sigmaNumber,0,3}]//Flatten

Solving for the Central Charges

ZAnsatzTerms[LocalFieldNumber][ii_,jj_][aa_,bb_]:=2 ZGen[LocalFieldNumber,1][ii,jj]Sum[γ5γdown[mu][[aa,bb]]D[A[t,x,y,z],Ϛ[mu]],{mu,0,3}]+2I ZGen[LocalFieldNumber,2][ii,jj]Sum[γdown[mu][[aa,bb]]D[B[t,x,y,z],Ϛ[mu]],{mu,0,3}] + 2 ZGen[LocalFieldNumber,3][ii,jj]γ5down[[aa,bb]]F[t,x,y,z]+ 2I ZGen[LocalFieldNumber,4][ii,jj]Cmetric[[aa,bb]]G[t,x,y,z]+ 2 I ZGen[LocalFieldNumber,5][ii,jj]Cmetric[[aa,bb]]d[t,x,y,z]+
2ZGen[LocalFieldNumber,6][ii,jj]γ5down[[aa,bb]]Sum[UD[Am[mu][t,x,y,z],Ϛ[mu]],{mu,0,3}] + (1/2) ZGen[LocalFieldNumber,7][ii,jj]Sum[(γ5.CommuteGammadown[mu,nu])[[aa,bb]]Fmn[mu,nu][t,x,y,z],{mu,0,3},{nu,0,3}]

ZAnsatz[LocalFieldNumber] = Table[ZAnsatzTerms[LocalFieldNumber][ii,jj][aa,bb],{ii,1,2},{aa,1,4},{jj,1,2},{bb,1,4}];

AlgebraDDMinusClosure[tt1,tt2][LocalFieldNumber]=(Collect[(AlgebraDD[tt1,tt2][LocalField] - ZAnsatz[LocalFieldNumber])/.LagrangianParameters/.SetW,FieldListWithDerivatives,Simplify]);

AlgebraDDEqs[tt1,tt2][LocalFieldNumber] = DeleteCases[Table[(Coefficient[AlgebraDDMinusClosure[tt1,tt2][LocalFieldNumber][[ii,aa,jj,bb]],FieldListWithDerivatives[[FieldListIndex]]])==0,{ii,1,2},{aa,1,4},{jj,1,2},{bb,1,4},{FieldListIndex,1,Length[FieldListWithDerivatives]}]//Simplify//Flatten,True]//DeleteDuplicates;

Length[AlgebraDDEqs[tt1,tt2][LocalFieldNumber]]

ZSolutions[LocalFieldNumber] = Solve[AlgebraDDEqs[tt1,tt2][LocalFieldNumber],CentralChargeList[LocalFieldNumber]]//Flatten

Final Results

DeleteCases[Table[Z[LocalFieldNumber,nn]==PrintSigmaProduct[Table[ZGen[LocalFieldNumber,nn][ii,jj],{ii,1,2},{jj,1,2}]/.ZSolutions[LocalFieldNumber]],{nn,1,LocalNumZCharges}],Z[_,_]==0]

ZTerms[LocalFieldNumber][ii_,jj_][aa_,bb_]:=-2 I sigma[2][[ii,jj]]Sum[γ5γdown[mu][[aa,bb]]D[A[t,x,y,z],Ϛ[mu]],{mu,0,3}]+ -2 I sigma[2][[ii,jj]]γ5down[[aa,bb]]F[t,x,y,z]+ (1/2) KroneckerDelta[ii,jj]Sum[(γ5.CommuteGammadown[mu,nu])[[aa,bb]]Fmn[mu,nu][t,x,y,z],{mu,0,3},{nu,0,3}]

Z[LocalFieldNumber] = Table[ZTerms[LocalFieldNumber][ii,jj][aa,bb],{ii,1,2},{aa,1,4},{jj,1,2},{bb,1,4}];

AlgebraDDResult[tt1,tt2][LocalFieldNumber] =(Collect[(AlgebraDD[tt1,tt2][LocalField] == Z[LocalFieldNumber])/.LagrangianParameters/.SetW,FieldListWithDerivatives,Simplify])

Clear[LocalField,LocalFieldNumber,LocalNumZCharges]

{(Subscript[D, a]^i),(Subscript[(OverHat[D]), b]^j)}(OverHat[F])=2i ((((σ^3)))^ij)(Subscript[((γ^μ)), ab])(∂_μ)F +2i ((((σ^2)))^ij)(Subscript[((γ^5)(γ^μ)), ab])(∂_μ)d-2i ((((σ^1)))^ij) (Subscript[((γ^μ)), ab])(∂^ν)(Subscript[F, μν]) 
FieldNumber = 12

Set Local Field

LocalField =Fh[t,x,y,z];
LocalFieldNumber = 12;
LocalNumZCharges =6;
CentralChargeList[LocalFieldNumber] = Table[z[LocalFieldNumber,FieldNumber,sigmaNumber],{FieldNumber,1,LocalNumZCharges},{sigmaNumber,0,3}]//Flatten

Solving for the Central Charges

ZAnsatzTerms[LocalFieldNumber][ii_,jj_][aa_,bb_]:=2 I ZGen[LocalFieldNumber,1][ii,jj]Cmetric[[aa,bb]]Lap[A[t,x,y,z]]+2ZGen[LocalFieldNumber,2][ii,jj]γ5down[[aa,bb]]Lap[B[t,x,y,z]] + 2 I ZGen[LocalFieldNumber,3][ii,jj]Sum[γdown[mu][[aa,bb]]D[F[t,x,y,z],Ϛ[mu]],{mu,0,3}]+ 2 ZGen[LocalFieldNumber,4][ii,jj]Sum[γ5γdown[mu][[aa,bb]]D[G[t,x,y,z],Ϛ[mu]],{mu,0,3}]+ 2 ZGen[LocalFieldNumber,5][ii,jj]Sum[γ5γdown[mu][[aa,bb]]D[d[t,x,y,z],Ϛ[mu]],{mu,0,3}]+  I ZGen[LocalFieldNumber,6][ii,jj]Sum[γdown[mu][[aa,bb]]UD[Fmn[mu,nu][t,x,y,z],Ϛ[nu]],{mu,0,3},{nu,0,3}]

ZAnsatz[LocalFieldNumber] = Table[ZAnsatzTerms[LocalFieldNumber][ii,jj][aa,bb],{ii,1,2},{aa,1,4},{jj,1,2},{bb,1,4}];

AlgebraDDMinusClosure[tt1,tt2][LocalFieldNumber]=(Collect[(AlgebraDD[tt1,tt2][LocalField] - ZAnsatz[LocalFieldNumber])/.LagrangianParameters/.SetW,FieldListWithDerivatives,Simplify]);

AlgebraDDEqs[tt1,tt2][LocalFieldNumber] = DeleteCases[Table[(Coefficient[AlgebraDDMinusClosure[tt1,tt2][LocalFieldNumber][[ii,aa,jj,bb]],FieldListWithDerivatives[[FieldListIndex]]])==0,{ii,1,2},{aa,1,4},{jj,1,2},{bb,1,4},{FieldListIndex,1,Length[FieldListWithDerivatives]}]//Simplify//Flatten,True]//DeleteDuplicates;

AlgebraDDEqs[tt1,tt2][LocalFieldNumber]

Length[AlgebraDDEqs[tt1,tt2][LocalFieldNumber]]

ZSolutions[LocalFieldNumber] = Solve[AlgebraDDEqs[tt1,tt2][LocalFieldNumber],CentralChargeList[LocalFieldNumber]]//Flatten

Final Results

DeleteCases[Table[Z[LocalFieldNumber,nn]==PrintSigmaProduct[Table[ZGen[LocalFieldNumber,nn][ii,jj],{ii,1,2},{jj,1,2}]/.ZSolutions[LocalFieldNumber]],{nn,1,LocalNumZCharges}],Z[_,_]==0]

ZTerms[LocalFieldNumber][ii_,jj_][aa_,bb_]:= 2 I sigma[3][[ii,jj]]Sum[γdown[mu][[aa,bb]]D[F[t,x,y,z],Ϛ[mu]],{mu,0,3}]+ 2 I sigma[2][[ii,jj]]Sum[γ5γdown[mu][[aa,bb]]D[d[t,x,y,z],Ϛ[mu]],{mu,0,3}]-2I sigma[1][[ii,jj]]Sum[γdown[mu][[aa,bb]]UD[Fmn[mu,nu][t,x,y,z],Ϛ[nu]],{mu,0,3},{nu,0,3}]

Z[LocalFieldNumber] = Table[ZTerms[LocalFieldNumber][ii,jj][aa,bb],{ii,1,2},{aa,1,4},{jj,1,2},{bb,1,4}];

AlgebraDDResult[tt1,tt2][LocalFieldNumber] =(Collect[(AlgebraDD[tt1,tt2][LocalField] == Z[LocalFieldNumber])/.LagrangianParameters/.SetW,FieldListWithDerivatives,Simplify])

Clear[LocalField,LocalFieldNumber,LocalNumZCharges]

{(Subscript[D, a]^i),(Subscript[(OverHat[D]), b]^j)}(Overscript[F, ⌣])=2i ((((σ^1)))^ij)(Subscript[((γ^μ)), ab])(∂_μ)F -2i ((((σ^2)))^ij)(Subscript[((γ^5)(γ^μ)), ab])(∂_μ)G+2i ((((σ^3)))^ij) (Subscript[((γ^μ)), ab])(∂^ν)(Subscript[F, μν]) 
FieldNumber = 13

Set Local Field

LocalField =Fb[t,x,y,z];
LocalFieldNumber = 13;
LocalNumZCharges =6;
CentralChargeList[LocalFieldNumber] = Table[z[LocalFieldNumber,FieldNumber,sigmaNumber],{FieldNumber,1,LocalNumZCharges},{sigmaNumber,0,3}]//Flatten

Solving for the Central Charges

ZAnsatzTerms[LocalFieldNumber][ii_,jj_][aa_,bb_]:=2 I ZGen[LocalFieldNumber,1][ii,jj]Cmetric[[aa,bb]]Lap[A[t,x,y,z]]+2ZGen[LocalFieldNumber,2][ii,jj]γ5down[[aa,bb]]Lap[B[t,x,y,z]] + 2 I ZGen[LocalFieldNumber,3][ii,jj]Sum[γdown[mu][[aa,bb]]D[F[t,x,y,z],Ϛ[mu]],{mu,0,3}]+ 2 ZGen[LocalFieldNumber,4][ii,jj]Sum[γ5γdown[mu][[aa,bb]]D[G[t,x,y,z],Ϛ[mu]],{mu,0,3}]+ 2 ZGen[LocalFieldNumber,5][ii,jj]Sum[γ5γdown[mu][[aa,bb]]D[d[t,x,y,z],Ϛ[mu]],{mu,0,3}]+  I ZGen[LocalFieldNumber,6][ii,jj]Sum[γdown[mu][[aa,bb]]UD[Fmn[mu,nu][t,x,y,z],Ϛ[nu]],{mu,0,3},{nu,0,3}]

ZAnsatz[LocalFieldNumber] = Table[ZAnsatzTerms[LocalFieldNumber][ii,jj][aa,bb],{ii,1,2},{aa,1,4},{jj,1,2},{bb,1,4}];

AlgebraDDMinusClosure[tt1,tt2][LocalFieldNumber]=(Collect[(AlgebraDD[tt1,tt2][LocalField] - ZAnsatz[LocalFieldNumber])/.LagrangianParameters/.SetW,FieldListWithDerivatives,Simplify]);

AlgebraDDEqs[tt1,tt2][LocalFieldNumber] = DeleteCases[Table[(Coefficient[AlgebraDDMinusClosure[tt1,tt2][LocalFieldNumber][[ii,aa,jj,bb]],FieldListWithDerivatives[[FieldListIndex]]])==0,{ii,1,2},{aa,1,4},{jj,1,2},{bb,1,4},{FieldListIndex,1,Length[FieldListWithDerivatives]}]//Simplify//Flatten,True]//DeleteDuplicates;

AlgebraDDEqs[tt1,tt2][LocalFieldNumber]

Length[AlgebraDDEqs[tt1,tt2][LocalFieldNumber]]

ZSolutions[LocalFieldNumber] = Solve[AlgebraDDEqs[tt1,tt2][LocalFieldNumber],CentralChargeList[LocalFieldNumber]]//Flatten

Final Results

DeleteCases[Table[Z[LocalFieldNumber,nn]==PrintSigmaProduct[Table[ZGen[LocalFieldNumber,nn][ii,jj],{ii,1,2},{jj,1,2}]/.ZSolutions[LocalFieldNumber]],{nn,1,LocalNumZCharges}],Z[_,_]==0]

ZTerms[LocalFieldNumber][ii_,jj_][aa_,bb_]:= 2 I sigma[1][[ii,jj]]Sum[γdown[mu][[aa,bb]]D[F[t,x,y,z],Ϛ[mu]],{mu,0,3}]-2 I sigma[2][[ii,jj]]Sum[γ5γdown[mu][[aa,bb]]D[G[t,x,y,z],Ϛ[mu]],{mu,0,3}]+ 2I sigma[3][[ii,jj]]Sum[γdown[mu][[aa,bb]]UD[Fmn[mu,nu][t,x,y,z],Ϛ[nu]],{mu,0,3},{nu,0,3}]

Z[LocalFieldNumber] = Table[ZTerms[LocalFieldNumber][ii,jj][aa,bb],{ii,1,2},{aa,1,4},{jj,1,2},{bb,1,4}];

AlgebraDDResult[tt1,tt2][LocalFieldNumber] =(Collect[(AlgebraDD[tt1,tt2][LocalField] == Z[LocalFieldNumber])/.LagrangianParameters/.SetW,FieldListWithDerivatives,Simplify])

Clear[LocalField,LocalFieldNumber,LocalNumZCharges]

{(Subscript[D, a]^i),(Subscript[(OverHat[D]), b]^j)}(OverHat[G])=2i ((((σ^3)))^ij)(Subscript[((γ^μ)), ab])(∂_μ)G +2i ((((σ^1)))^ij)(Subscript[((γ^μ)), ab])(∂_μ)d+2i ((((σ^2)))^ij) (Subscript[((γ^5)(γ^μ)), ab])(∂^ν)(Subscript[F, μν]) 
FieldNumber = 14

Set Local Field

LocalField =Gh[t,x,y,z];
LocalFieldNumber = 14;
LocalNumZCharges =6;
CentralChargeList[LocalFieldNumber] = Table[z[LocalFieldNumber,FieldNumber,sigmaNumber],{FieldNumber,1,LocalNumZCharges},{sigmaNumber,0,3}]//Flatten

Solving for the Central Charges

ZAnsatzTerms[LocalFieldNumber][ii_,jj_][aa_,bb_]:=2 ZGen[LocalFieldNumber,1][ii,jj]γ5down[[aa,bb]]Lap[A[t,x,y,z]]+2I ZGen[LocalFieldNumber,2][ii,jj]Cmetric[[aa,bb]]Lap[B[t,x,y,z]] + 2 ZGen[LocalFieldNumber,3][ii,jj]Sum[γ5γdown[mu][[aa,bb]]D[F[t,x,y,z],Ϛ[mu]],{mu,0,3}]+ 2 I ZGen[LocalFieldNumber,4][ii,jj]Sum[γdown[mu][[aa,bb]]D[G[t,x,y,z],Ϛ[mu]],{mu,0,3}]+ 2 I ZGen[LocalFieldNumber,5][ii,jj]Sum[γdown[mu][[aa,bb]]D[d[t,x,y,z],Ϛ[mu]],{mu,0,3}]+  2ZGen[LocalFieldNumber,6][ii,jj]Sum[γ5γdown[mu][[aa,bb]]UD[Fmn[mu,nu][t,x,y,z],Ϛ[nu]],{mu,0,3},{nu,0,3}]

ZAnsatz[LocalFieldNumber] = Table[ZAnsatzTerms[LocalFieldNumber][ii,jj][aa,bb],{ii,1,2},{aa,1,4},{jj,1,2},{bb,1,4}];

AlgebraDDMinusClosure[tt1,tt2][LocalFieldNumber]=(Collect[(AlgebraDD[tt1,tt2][LocalField] - ZAnsatz[LocalFieldNumber])/.LagrangianParameters/.SetW,FieldListWithDerivatives,Simplify]);

AlgebraDDEqs[tt1,tt2][LocalFieldNumber] = DeleteCases[Table[(Coefficient[AlgebraDDMinusClosure[tt1,tt2][LocalFieldNumber][[ii,aa,jj,bb]],FieldListWithDerivatives[[FieldListIndex]]])==0,{ii,1,2},{aa,1,4},{jj,1,2},{bb,1,4},{FieldListIndex,1,Length[FieldListWithDerivatives]}]//Simplify//Flatten,True]//DeleteDuplicates;

AlgebraDDEqs[tt1,tt2][LocalFieldNumber]

Length[AlgebraDDEqs[tt1,tt2][LocalFieldNumber]]

ZSolutions[LocalFieldNumber] = Solve[AlgebraDDEqs[tt1,tt2][LocalFieldNumber],CentralChargeList[LocalFieldNumber]]//Flatten

Final Results

DeleteCases[Table[Z[LocalFieldNumber,nn]==PrintSigmaProduct[Table[ZGen[LocalFieldNumber,nn][ii,jj],{ii,1,2},{jj,1,2}]/.ZSolutions[LocalFieldNumber]],{nn,1,LocalNumZCharges}],Z[_,_]==0]

ZTerms[LocalFieldNumber][ii_,jj_][aa_,bb_]:= 2 I sigma[3][[ii,jj]]Sum[γdown[mu][[aa,bb]]D[G[t,x,y,z],Ϛ[mu]],{mu,0,3}]+ 2 I sigma[1][[ii,jj]]Sum[γdown[mu][[aa,bb]]D[d[t,x,y,z],Ϛ[mu]],{mu,0,3}]+  2I sigma[2][[ii,jj]]Sum[γ5γdown[mu][[aa,bb]]UD[Fmn[mu,nu][t,x,y,z],Ϛ[nu]],{mu,0,3},{nu,0,3}]

Z[LocalFieldNumber] = Table[ZTerms[LocalFieldNumber][ii,jj][aa,bb],{ii,1,2},{aa,1,4},{jj,1,2},{bb,1,4}];

AlgebraDDResult[tt1,tt2][LocalFieldNumber] =(Collect[(AlgebraDD[tt1,tt2][LocalField] == Z[LocalFieldNumber])/.LagrangianParameters/.SetW,FieldListWithDerivatives,Simplify])

Clear[LocalField,LocalFieldNumber,LocalNumZCharges]

{(Subscript[D, a]^i),(Subscript[(OverHat[D]), b]^j)}(Overscript[G, ⌣])=-2i ((((σ^2)))^ij)(Subscript[((γ^5)), ab])□A+2((((σ^2)))^ij)(Subscript[C, ab])□B-2i ((((σ^2)))^ij)(Subscript[((γ^5)(γ^μ)), ab])(∂_μ)F-2i ((((σ^1)))^ij)(Subscript[((γ^μ)), ab])(∂_μ)G +2i ((((σ^3)))^ij)(Subscript[((γ^μ)), ab])(∂_μ)d 
FieldNumber = 15

Set Local Field

LocalField =Gb[t,x,y,z];
LocalFieldNumber = 15;
LocalNumZCharges =6;
CentralChargeList[LocalFieldNumber] = Table[z[LocalFieldNumber,FieldNumber,sigmaNumber],{FieldNumber,1,LocalNumZCharges},{sigmaNumber,0,3}]//Flatten

Solving for the Central Charges

ZAnsatzTerms[LocalFieldNumber][ii_,jj_][aa_,bb_]:=2 ZGen[LocalFieldNumber,1][ii,jj]γ5down[[aa,bb]]Lap[A[t,x,y,z]]+2I ZGen[LocalFieldNumber,2][ii,jj]Cmetric[[aa,bb]]Lap[B[t,x,y,z]] + 2 ZGen[LocalFieldNumber,3][ii,jj]Sum[γ5γdown[mu][[aa,bb]]D[F[t,x,y,z],Ϛ[mu]],{mu,0,3}]+ 2 I ZGen[LocalFieldNumber,4][ii,jj]Sum[γdown[mu][[aa,bb]]D[G[t,x,y,z],Ϛ[mu]],{mu,0,3}]+ 2 I ZGen[LocalFieldNumber,5][ii,jj]Sum[γdown[mu][[aa,bb]]D[d[t,x,y,z],Ϛ[mu]],{mu,0,3}]+  2ZGen[LocalFieldNumber,6][ii,jj]Sum[γ5γdown[mu][[aa,bb]]UD[Fmn[mu,nu][t,x,y,z],Ϛ[nu]],{mu,0,3},{nu,0,3}]

ZAnsatz[LocalFieldNumber] = Table[ZAnsatzTerms[LocalFieldNumber][ii,jj][aa,bb],{ii,1,2},{aa,1,4},{jj,1,2},{bb,1,4}];

AlgebraDDMinusClosure[tt1,tt2][LocalFieldNumber]=(Collect[(AlgebraDD[tt1,tt2][LocalField] - ZAnsatz[LocalFieldNumber])/.LagrangianParameters/.SetW,FieldListWithDerivatives,Simplify]);

AlgebraDDEqs[tt1,tt2][LocalFieldNumber] = DeleteCases[Table[(Coefficient[AlgebraDDMinusClosure[tt1,tt2][LocalFieldNumber][[ii,aa,jj,bb]],FieldListWithDerivatives[[FieldListIndex]]])==0,{ii,1,2},{aa,1,4},{jj,1,2},{bb,1,4},{FieldListIndex,1,Length[FieldListWithDerivatives]}]//Simplify//Flatten,True]//DeleteDuplicates;

AlgebraDDEqs[tt1,tt2][LocalFieldNumber]

Length[AlgebraDDEqs[tt1,tt2][LocalFieldNumber]]

ZSolutions[LocalFieldNumber] = Solve[AlgebraDDEqs[tt1,tt2][LocalFieldNumber],CentralChargeList[LocalFieldNumber]]//Flatten

Final Results

DeleteCases[Table[Z[LocalFieldNumber,nn]==PrintSigmaProduct[Table[ZGen[LocalFieldNumber,nn][ii,jj],{ii,1,2},{jj,1,2}]/.ZSolutions[LocalFieldNumber]],{nn,1,LocalNumZCharges}],Z[_,_]==0]

ZTerms[LocalFieldNumber][ii_,jj_][aa_,bb_]:= -2 I sigma[2][[ii,jj]]γ5down[[aa,bb]]Lap[A[t,x,y,z]]+2sigma[2][[ii,jj]]Cmetric[[aa,bb]]Lap[B[t,x,y,z]] - 2I sigma[2][[ii,jj]]Sum[γ5γdown[mu][[aa,bb]]D[F[t,x,y,z],Ϛ[mu]],{mu,0,3}]-2 I sigma[1][[ii,jj]]Sum[γdown[mu][[aa,bb]]D[G[t,x,y,z],Ϛ[mu]],{mu,0,3}]+ 2 I sigma[3][[ii,jj]]Sum[γdown[mu][[aa,bb]]D[d[t,x,y,z],Ϛ[mu]],{mu,0,3}]

Z[LocalFieldNumber] = Table[ZTerms[LocalFieldNumber][ii,jj][aa,bb],{ii,1,2},{aa,1,4},{jj,1,2},{bb,1,4}];

AlgebraDDResult[tt1,tt2][LocalFieldNumber] =(Collect[(AlgebraDD[tt1,tt2][LocalField] == Z[LocalFieldNumber])/.LagrangianParameters/.SetW,FieldListWithDerivatives,Simplify])

LocalField

Clear[LocalField,LocalFieldNumber,LocalNumZCharges]

Boson Closure summary of N=2 chiral-vector unhat-hat algebra (tt1 =1, tt2 = 0)

Table[AlgebraDDResult[tt1,tt2][FieldIndex],{FieldIndex,{1,2,3,4,5,7,8,9,10,11,12,13,14,15}}]==Table[True,{FieldIndex,{1,2,3,4,5,7,8,9,10,11,12,13,14,15}}]

tt1
tt2

Fermion Closure for N=2 chiral-vector unhat-hat algebra (tt1 =1, tt2 = 0)

{(Subscript[D, a]^i),(Subscript[(OverHat[D]), b]^j)}(Subscript[Ψ, c]^k)= 2 i ((z[lfn,1])^ijkl)(Subscript[((γ^μ)), ab])(∂_μ)(Subscript[(OverHat[Ψ]), c]^l)+2 i ((z[lfn,2])^ijkl)(Subscript[([(γ^μ),(γ^ν)]), ab])((Subscript[((Subscript[γ, ν])), c])^d)(∂_μ)(Subscript[(OverHat[Ψ]), d]^l)
+2 i ((z[lfn,3])^ijkl)(Subscript[((Subscript[γ, ν])), ab])((Subscript[([(γ^ν),(γ^μ)]), c])^d)(∂_μ)(Subscript[(OverHat[Ψ]), d]^l)+2 i ((z[lfn,4])^ijkl)(Subscript[((γ^5)[(γ^μ),(γ^ν)]), ab])((Subscript[((γ^5)(Subscript[γ, ν])), c])^d)(∂_μ)(Subscript[(OverHat[Ψ]), d]^l)
+2 i ((z[lfn,5])^ijkl)(Subscript[((γ^5)(γ^μ)), ab])((Subscript[((γ^5)), c])^d)(∂_μ)(Subscript[(OverHat[Ψ]), d]^l)+2 i ((z[lfn,6])^ijkl)(Subscript[((γ^5)(Subscript[γ, ν])), ab])((Subscript[((γ^5)[(γ^ν),(γ^μ)]), c])^d)(∂_μ)(Subscript[(OverHat[Ψ]), d]^l)
+2 i ((z[lfn,7])^ijkl)(Subscript[((γ^5)), ab])((Subscript[((γ^5)(γ^μ)), c])^d)(∂_μ)(Subscript[(OverHat[Ψ]), d]^l)+2 i ((z[lfn,8])^ijkl)((Subscript[(Subscript[C, ab])((γ^μ)), c])^d)(∂_μ)(Subscript[(OverHat[Ψ]), d]^l)
FieldNumber = 6

Set Local Field

LocalFieldNumber = 6;
NumZCharges = 8;
CentralChargeList[LocalFieldNumber] = Table[z[LocalFieldNumber,FieldNumber][ii,jj,kk,ll],{FieldNumber,1,NumZCharges},{ii,1,2},{jj,1,2},{kk,1,2},{ll,1,2}]//Flatten;

Length[CentralChargeList[LocalFieldNumber]]==128

Solving for Z:

ZAnsatzTerms[LocalFieldNumber][ii_,jj_,kk_][aa_,bb_,cc_]:=2 I Sum[z[LocalFieldNumber,1][ii,jj,kk,ll]γdown[mu][[aa,bb]]D[Ψh[ll][cc][t,x,y,z],Ϛ[mu]],{ll,1,2},{mu,0,3}]+2I Sum[z[LocalFieldNumber,2][ii,jj,kk,ll]CommuteGammadown[mu,nu][[aa,bb]]γstdown[nu][[cc,dd]]D[Ψh[ll][dd][t,x,y,z],Ϛ[mu]],{ll,1,2},{dd,1,4},{mu,0,3},{nu,0,3}]+2I Sum[z[LocalFieldNumber,3][ii,jj,kk,ll]γstdowndown[nu][[aa,bb]]CommuteGamma[nu,mu][[cc,dd]]D[Ψh[ll][dd][t,x,y,z],Ϛ[mu]],{ll,1,2},{dd,1,4},{mu,0,3},{nu,0,3}]+2I Sum[z[LocalFieldNumber,4][ii,jj,kk,ll](γ5.CommuteGammadown[mu,nu])[[aa,bb]]γ5γstdown[nu][[cc,dd]]D[Ψh[ll][dd][t,x,y,z],Ϛ[mu]],{ll,1,2},{dd,1,4},{mu,0,3},{nu,0,3}]+2 I Sum[z[LocalFieldNumber,5][ii,jj,kk,ll]γ5γdown[mu][[aa,bb]]γ5[[cc,dd]]D[Ψh[ll][dd][t,x,y,z],Ϛ[mu]],{ll,1,2},{mu,0,3},{dd,1,4}]+2I Sum[z[LocalFieldNumber,6][ii,jj,kk,ll]γ5γstdowndown[nu][[aa,bb]](γ5.CommuteGamma[nu,mu])[[cc,dd]]D[Ψh[ll][dd][t,x,y,z],Ϛ[mu]],{ll,1,2},{dd,1,4},{mu,0,3},{nu,0,3}]+2 I Sum[z[LocalFieldNumber,7][ii,jj,kk,ll]γ5down[[aa,bb]]γ5γ[mu][[cc,dd]]D[Ψh[ll][dd][t,x,y,z],Ϛ[mu]],{ll,1,2},{mu,0,3},{dd,1,4}]+2 I Sum[z[LocalFieldNumber,8][ii,jj,kk,ll]Cmetric[[aa,bb]]γ[mu][[cc,dd]]D[Ψh[ll][dd][t,x,y,z],Ϛ[mu]],{ll,1,2},{mu,0,3},{dd,1,4}]

Do[ZAnsatz[LocalFieldNumber][kk,cc] = Table[ZAnsatzTerms[LocalFieldNumber][ii,jj,kk][aa,bb,cc],{ii,1,2},{aa,1,4},{jj,1,2},{bb,1,4}],{kk,1,2},{cc,1,4}];

AlgebraDDMinusClosure[tt1,tt2][LocalFieldNumber]=Table[(Collect[(AlgebraDD[tt1,tt2][Ψ[kk][cc][t,x,y,z]] - ZAnsatz[LocalFieldNumber][kk,cc])/.LagrangianParameters/.SetW,FieldListWithDerivatives,Simplify]),{kk,1,2},{cc,1,4}];

AlgebraDDEqs[tt1,tt2][LocalFieldNumber] = DeleteCases[Table[(Coefficient[AlgebraDDMinusClosure[tt1,tt2][LocalFieldNumber][[kk,cc,ii,aa,jj,bb]],FieldListWithDerivatives[[FieldListIndex]]])==0,{kk,1,2},{cc,1,4},{ii,1,2},{aa,1,4},{jj,1,2},{bb,1,4},{FieldListIndex,1,Length[FieldListWithDerivatives]}]//Simplify//Flatten,True]//DeleteDuplicates;

ZSolutions[LocalFieldNumber] = Solve[AlgebraDDEqs[tt1,tt2][LocalFieldNumber],CentralChargeList[LocalFieldNumber]]//Flatten;

Dimensions[ZSolutions[LocalFieldNumber]]

NonZeroz[LocalFieldNumber] =DeleteCases[Table[z[LocalFieldNumber,nn][ii,jj,kk,ll]==(z[LocalFieldNumber,nn][ii,jj,kk,ll]/.ZSolutions[LocalFieldNumber]),{nn,1,NumZCharges},{ii,1,2},{jj,1,2},{kk,1,2},{ll,1,2}]//Flatten,z[_,_][_,_,_,_]==0]

zAnsatz[lfn_,nn_][ii_,jj_,kk_,ll_]:=cz[lfn,nn,1]KroneckerDelta[ii,jj]sigma[3][[kk,ll]]+cz[lfn,nn,2]KroneckerDelta[ii,kk]sigma[3][[jj,ll]] +cz[lfn,nn,3]KroneckerDelta[ii,ll]sigma[3][[jj,kk]]+cz[lfn,nn,4]KroneckerDelta[jj,kk]sigma[3][[ii,ll]]+cz[lfn,nn,5]KroneckerDelta[jj,ll]sigma[3][[ii,kk]]+cz[lfn,nn,6]KroneckerDelta[kk,ll]sigma[3][[ii,jj]]+I cz[lfn,nn,7]sigma[2][[ii,jj]]sigma[1][[kk,ll]]+I cz[lfn,nn,8]sigma[2][[ii,kk]]sigma[1][[jj,ll]]+I cz[lfn,nn,9]sigma[2][[ii,ll]]sigma[1][[jj,kk]]+I cz[lfn,nn,10]sigma[2][[jj,kk]]sigma[1][[ii,ll]]+I cz[lfn,nn,11]sigma[2][[jj,ll]]sigma[1][[ii,kk]]+I cz[lfn,nn,12]sigma[2][[kk,ll]]sigma[1][[ii,jj]]+cz[lfn,nn,13]KroneckerDelta[ii,jj]KroneckerDelta[kk,ll]+cz[lfn,nn,14]KroneckerDelta[ii,kk]KroneckerDelta[jj,ll] +cz[lfn,nn,15]KroneckerDelta[ii,ll]KroneckerDelta[jj,kk]+cz[lfn,nn,16]sigma[3][[ii,jj]]sigma[3][[kk,ll]]+cz[lfn,nn,17]sigma[3][[ii,kk]]sigma[3][[jj,ll]] +cz[lfn,nn,18]sigma[3][[ii,ll]]sigma[3][[jj,kk]]+cz[lfn,nn,19]sigma[2][[ii,jj]]sigma[2][[kk,ll]]+cz[lfn,nn,20]sigma[2][[ii,kk]]sigma[2][[jj,ll]] +cz[lfn,nn,21]sigma[2][[ii,ll]]sigma[2][[jj,kk]]+cz[lfn,nn,22]sigma[1][[ii,jj]]sigma[1][[kk,ll]]+cz[lfn,nn,23]sigma[1][[ii,kk]]sigma[1][[jj,ll]] +cz[lfn,nn,24]sigma[1][[ii,ll]]sigma[1][[jj,kk]]

czList[lfn_,nn_]:= Table[cz[lfn,nn,Index],{Index,1,24}]

zAnsatzEqs[LocalFieldNumber]=DeleteCases[Table[(zAnsatz[LocalFieldNumber,ChargeNum][ii,jj,kk,ll])==(z[LocalFieldNumber,ChargeNum][ii,jj,kk,ll]/.ZSolutions[LocalFieldNumber]),{ChargeNum,1,NumZCharges},{ii,1,2},{jj,1,2},{kk,1,2},{ll,1,2}]//Flatten,True];

Length[zAnsatzEqs[LocalFieldNumber]]

Zerocz[LocalFieldNumber] = Delete[Table[cz[LocalFieldNumber,ChargeNum,CoeffNum]->0,{ChargeNum,1,NumZCharges},{CoeffNum,1,24}],Flatten[Table[{ChargeNum,NonZeroCoeff},{ChargeNum,1,NumZCharges},{NonZeroCoeff,{18,22,23,24}}],1]]//Flatten

czSolution[LocalFieldNumber]=Join[Solve[zAnsatzEqs[LocalFieldNumber]/.Zerocz[LocalFieldNumber]//Flatten,Table[czList[LocalFieldNumber,ChargeNum],{ChargeNum,1,NumZCharges}]//Flatten],Zerocz[LocalFieldNumber]]//Flatten

((Z[lfn,nn])^ijkl) = cz[lfn,nn,18]((((σ^3)))^il)((((σ^3)))^jk)+cz[lfn,nn,22]((((σ^1)))^ij)((((σ^1)))^kl)+cz[lfn,nn,23]((((σ^1)))^ik)((((σ^1)))^jl)+cz[lfn,nn,24]((((σ^1)))^il)((((σ^1)))^jk)

z[LocalFieldNumber,nn_][ii_,jj_,kk_,ll_]:=cz[LocalFieldNumber,nn,18]sigma[3][[ii,ll]]sigma[3][[jj,kk]]+cz[LocalFieldNumber,nn,22]sigma[1][[ii,jj]]sigma[1][[kk,ll]]+cz[LocalFieldNumber,nn,23]sigma[1][[ii,kk]]sigma[1][[jj,ll]] +cz[LocalFieldNumber,nn,24]sigma[1][[ii,ll]]sigma[1][[jj,kk]]/.czSolution[LocalFieldNumber]

z[LocalFieldNumber,1][1,2,1,2]

ZTerms[LocalFieldNumber][ii_,jj_,kk_][aa_,bb_,cc_]:=2 I Sum[z[LocalFieldNumber,1][ii,jj,kk,ll]γdown[mu][[aa,bb]]D[Ψh[ll][cc][t,x,y,z],Ϛ[mu]],{ll,1,2},{mu,0,3}]+2I Sum[z[LocalFieldNumber,2][ii,jj,kk,ll]CommuteGammadown[mu,nu][[aa,bb]]γstdown[nu][[cc,dd]]D[Ψh[ll][dd][t,x,y,z],Ϛ[mu]],{ll,1,2},{dd,1,4},{mu,0,3},{nu,0,3}]+2I Sum[z[LocalFieldNumber,3][ii,jj,kk,ll]γstdowndown[nu][[aa,bb]]CommuteGamma[nu,mu][[cc,dd]]D[Ψh[ll][dd][t,x,y,z],Ϛ[mu]],{ll,1,2},{dd,1,4},{mu,0,3},{nu,0,3}]+2I Sum[z[LocalFieldNumber,4][ii,jj,kk,ll](γ5.CommuteGammadown[mu,nu])[[aa,bb]]γ5γstdown[nu][[cc,dd]]D[Ψh[ll][dd][t,x,y,z],Ϛ[mu]],{ll,1,2},{dd,1,4},{mu,0,3},{nu,0,3}]+2 I Sum[z[LocalFieldNumber,5][ii,jj,kk,ll]γ5γdown[mu][[aa,bb]]γ5[[cc,dd]]D[Ψh[ll][dd][t,x,y,z],Ϛ[mu]],{ll,1,2},{mu,0,3},{dd,1,4}]+2I Sum[z[LocalFieldNumber,6][ii,jj,kk,ll]γ5γstdowndown[nu][[aa,bb]](γ5.CommuteGamma[nu,mu])[[cc,dd]]D[Ψh[ll][dd][t,x,y,z],Ϛ[mu]],{ll,1,2},{dd,1,4},{mu,0,3},{nu,0,3}]+2 I Sum[z[LocalFieldNumber,7][ii,jj,kk,ll]γ5down[[aa,bb]]γ5γ[mu][[cc,dd]]D[Ψh[ll][dd][t,x,y,z],Ϛ[mu]],{ll,1,2},{mu,0,3},{dd,1,4}]+2 I Sum[z[LocalFieldNumber,8][ii,jj,kk,ll]Cmetric[[aa,bb]]γ[mu][[cc,dd]]D[Ψh[ll][dd][t,x,y,z],Ϛ[mu]],{ll,1,2},{mu,0,3},{dd,1,4}]

Do[Z[LocalFieldNumber][kk,cc] = Table[ZTerms[LocalFieldNumber][ii,jj,kk][aa,bb,cc],{ii,1,2},{aa,1,4},{jj,1,2},{bb,1,4}],{kk,1,2},{cc,1,4}];

AlgebraDDResult[tt1,tt2][LocalFieldNumber]=Table[AlgebraDD[tt1,tt2][Ψ[kk][cc][t,x,y,z]]/.LagrangianParameters/.SetW,{kk,1,2},{cc,1,4}] ==Table[Z[LocalFieldNumber][kk,cc],{kk,1,2},{cc,1,4}]//Simplify

ZPrintout[LocalFieldNumber]= Table[(Subscript[(OverHat[Z]), nn]^ijkl)== (cz[LocalFieldNumber,nn,18]((((σ^3)))^il)((((σ^("3"))))^jk) + cz[LocalFieldNumber,nn,22]((((σ^(HoldForm[1]))))^ij)((((σ^("1"))))^kl) + cz[LocalFieldNumber,nn,23]((((σ^(HoldForm[1]))))^ik)((((σ^("1"))))^jl)+cz[LocalFieldNumber,nn,24]((((σ^(HoldForm[1]))))^il)((((σ^("1"))))^jk)  )/.czSolution[LocalFieldNumber],{nn,1,NumZCharges}]
ZPrintoutTeXForm[LocalFieldNumber ] = TeXForm[%]

Clear[LocalFieldNumber,NumZCharges]

{(Subscript[D, a]^i),(Subscript[(OverHat[D]), b]^j)}(Subscript[(OverHat[Ψ]), c]^k)= 2 i ((z[lfn,1])^ijkl)(Subscript[((γ^μ)), ab])(∂_μ)(Subscript[Ψ, c]^l)+2 i ((z[lfn,2])^ijkl)(Subscript[([(γ^μ),(γ^ν)]), ab])((Subscript[((Subscript[γ, ν])), c])^d)(∂_μ)(Subscript[Ψ, d]^l)
+2 i ((z[lfn,3])^ijkl)(Subscript[((Subscript[γ, ν])), ab])((Subscript[([(γ^ν),(γ^μ)]), c])^d)(∂_μ)(Subscript[Ψ, d]^l)+2 i ((z[lfn,4])^ijkl)(Subscript[((γ^5)[(γ^μ),(γ^ν)]), ab])((Subscript[((γ^5)(Subscript[γ, ν])), c])^d)(∂_μ)(Subscript[Ψ, d]^l)
+2 i ((z[lfn,5])^ijkl)(Subscript[((γ^5)(γ^μ)), ab])((Subscript[((γ^5)), c])^d)(∂_μ)(Subscript[Ψ, d]^l)+2 i ((z[lfn,6])^ijkl)(Subscript[((γ^5)(Subscript[γ, ν])), ab])((Subscript[((γ^5)[(γ^ν),(γ^μ)]), c])^d)(∂_μ)(Subscript[Ψ, d]^l)
+2 i ((z[lfn,7])^ijkl)(Subscript[((γ^5)), ab])((Subscript[((γ^5)(γ^μ)), c])^d)(∂_μ)(Subscript[Ψ, d]^l)+2 i ((z[lfn,8])^ijkl)((Subscript[(Subscript[C, ab])((γ^μ)), c])^d)(∂_μ)(Subscript[Ψ, d]^l)
FieldNumber = 16

Set Local Field

LocalFieldNumber = 16;
NumZCharges = 8;
CentralChargeList[LocalFieldNumber] = Table[z[LocalFieldNumber,FieldNumber][ii,jj,kk,ll],{FieldNumber,1,NumZCharges},{ii,1,2},{jj,1,2},{kk,1,2},{ll,1,2}]//Flatten;

Length[CentralChargeList[LocalFieldNumber]]==128

Solving for Z:

ZAnsatzTerms[LocalFieldNumber][ii_,jj_,kk_][aa_,bb_,cc_]:=2 I Sum[z[LocalFieldNumber,1][ii,jj,kk,ll]γdown[mu][[aa,bb]]D[Ψ[ll][cc][t,x,y,z],Ϛ[mu]],{ll,1,2},{mu,0,3}]+2I Sum[z[LocalFieldNumber,2][ii,jj,kk,ll]CommuteGammadown[mu,nu][[aa,bb]]γstdown[nu][[cc,dd]]D[Ψ[ll][dd][t,x,y,z],Ϛ[mu]],{ll,1,2},{dd,1,4},{mu,0,3},{nu,0,3}]+2I Sum[z[LocalFieldNumber,3][ii,jj,kk,ll]γstdowndown[nu][[aa,bb]]CommuteGamma[nu,mu][[cc,dd]]D[Ψ[ll][dd][t,x,y,z],Ϛ[mu]],{ll,1,2},{dd,1,4},{mu,0,3},{nu,0,3}]+2I Sum[z[LocalFieldNumber,4][ii,jj,kk,ll](γ5.CommuteGammadown[mu,nu])[[aa,bb]]γ5γstdown[nu][[cc,dd]]D[Ψ[ll][dd][t,x,y,z],Ϛ[mu]],{ll,1,2},{dd,1,4},{mu,0,3},{nu,0,3}]+2 I Sum[z[LocalFieldNumber,5][ii,jj,kk,ll]γ5γdown[mu][[aa,bb]]γ5[[cc,dd]]D[Ψ[ll][dd][t,x,y,z],Ϛ[mu]],{ll,1,2},{mu,0,3},{dd,1,4}]+2I Sum[z[LocalFieldNumber,6][ii,jj,kk,ll]γ5γstdowndown[nu][[aa,bb]](γ5.CommuteGamma[nu,mu])[[cc,dd]]D[Ψ[ll][dd][t,x,y,z],Ϛ[mu]],{ll,1,2},{dd,1,4},{mu,0,3},{nu,0,3}]+2 I Sum[z[LocalFieldNumber,7][ii,jj,kk,ll]γ5down[[aa,bb]]γ5γ[mu][[cc,dd]]D[Ψ[ll][dd][t,x,y,z],Ϛ[mu]],{ll,1,2},{mu,0,3},{dd,1,4}]+2 I Sum[z[LocalFieldNumber,8][ii,jj,kk,ll]Cmetric[[aa,bb]]γ[mu][[cc,dd]]D[Ψ[ll][dd][t,x,y,z],Ϛ[mu]],{ll,1,2},{mu,0,3},{dd,1,4}]

Do[ZAnsatz[LocalFieldNumber][kk,cc] = Table[ZAnsatzTerms[LocalFieldNumber][ii,jj,kk][aa,bb,cc],{ii,1,2},{aa,1,4},{jj,1,2},{bb,1,4}],{kk,1,2},{cc,1,4}];

AlgebraDDMinusClosure[tt1,tt2][LocalFieldNumber]=Table[(Collect[(AlgebraDD[tt1,tt2][Ψh[kk][cc][t,x,y,z]] - ZAnsatz[LocalFieldNumber][kk,cc])/.LagrangianParameters/.SetW,FieldListWithDerivatives,Simplify]),{kk,1,2},{cc,1,4}];

AlgebraDDEqs[tt1,tt2][LocalFieldNumber] = DeleteCases[Table[(Coefficient[AlgebraDDMinusClosure[tt1,tt2][LocalFieldNumber][[kk,cc,ii,aa,jj,bb]],FieldListWithDerivatives[[FieldListIndex]]])==0,{kk,1,2},{cc,1,4},{ii,1,2},{aa,1,4},{jj,1,2},{bb,1,4},{FieldListIndex,1,Length[FieldListWithDerivatives]}]//Simplify//Flatten,True]//DeleteDuplicates;

ZSolutions[LocalFieldNumber] = Solve[AlgebraDDEqs[tt1,tt2][LocalFieldNumber],CentralChargeList[LocalFieldNumber]]//Flatten;

Dimensions[ZSolutions[LocalFieldNumber]]

NonZeroz[LocalFieldNumber] =DeleteCases[Table[z[LocalFieldNumber,nn][ii,jj,kk,ll]==(z[LocalFieldNumber,nn][ii,jj,kk,ll]/.ZSolutions[LocalFieldNumber]),{nn,1,NumZCharges},{ii,1,2},{jj,1,2},{kk,1,2},{ll,1,2}]//Flatten,z[_,_][_,_,_,_]==0]

zAnsatzEqs[LocalFieldNumber]=DeleteCases[Table[(zAnsatz[LocalFieldNumber,ChargeNum][ii,jj,kk,ll])==(z[LocalFieldNumber,ChargeNum][ii,jj,kk,ll]/.ZSolutions[LocalFieldNumber]),{ChargeNum,1,NumZCharges},{ii,1,2},{jj,1,2},{kk,1,2},{ll,1,2}]//Flatten,True];

Length[zAnsatzEqs[LocalFieldNumber]]

Zerocz[LocalFieldNumber] = Delete[Table[cz[LocalFieldNumber,ChargeNum,CoeffNum]->0,{ChargeNum,1,NumZCharges},{CoeffNum,1,24}],Flatten[Table[{ChargeNum,NonZeroCoeff},{ChargeNum,1,NumZCharges},{NonZeroCoeff,{18,22,23,24}}],1]]//Flatten

czSolution[LocalFieldNumber]=Join[Solve[zAnsatzEqs[LocalFieldNumber]/.Zerocz[LocalFieldNumber]//Flatten,Table[czList[LocalFieldNumber,ChargeNum],{ChargeNum,1,NumZCharges}]//Flatten],Zerocz[LocalFieldNumber]]//Flatten

((Z[lfn,nn])^ijkl) = cz[lfn,nn,18]((((σ^3)))^il)((((σ^3)))^jk)+cz[lfn,nn,22]((((σ^1)))^ij)((((σ^1)))^kl)+cz[lfn,nn,23]((((σ^1)))^ik)((((σ^1)))^jl)+cz[lfn,nn,24]((((σ^1)))^il)((((σ^1)))^jk)

z[LocalFieldNumber,nn_][ii_,jj_,kk_,ll_]:=cz[LocalFieldNumber,nn,18]sigma[3][[ii,ll]]sigma[3][[jj,kk]]+cz[LocalFieldNumber,nn,22]sigma[1][[ii,jj]]sigma[1][[kk,ll]]+cz[LocalFieldNumber,nn,23]sigma[1][[ii,kk]]sigma[1][[jj,ll]] +cz[LocalFieldNumber,nn,24]sigma[1][[ii,ll]]sigma[1][[jj,kk]]/.czSolution[LocalFieldNumber]

z[LocalFieldNumber,1][1,2,1,2]

ZTerms[LocalFieldNumber][ii_,jj_,kk_][aa_,bb_,cc_]:=2 I Sum[z[LocalFieldNumber,1][ii,jj,kk,ll]γdown[mu][[aa,bb]]D[Ψ[ll][cc][t,x,y,z],Ϛ[mu]],{ll,1,2},{mu,0,3}]+2I Sum[z[LocalFieldNumber,2][ii,jj,kk,ll]CommuteGammadown[mu,nu][[aa,bb]]γstdown[nu][[cc,dd]]D[Ψ[ll][dd][t,x,y,z],Ϛ[mu]],{ll,1,2},{dd,1,4},{mu,0,3},{nu,0,3}]+2I Sum[z[LocalFieldNumber,3][ii,jj,kk,ll]γstdowndown[nu][[aa,bb]]CommuteGamma[nu,mu][[cc,dd]]D[Ψ[ll][dd][t,x,y,z],Ϛ[mu]],{ll,1,2},{dd,1,4},{mu,0,3},{nu,0,3}]+2I Sum[z[LocalFieldNumber,4][ii,jj,kk,ll](γ5.CommuteGammadown[mu,nu])[[aa,bb]]γ5γstdown[nu][[cc,dd]]D[Ψ[ll][dd][t,x,y,z],Ϛ[mu]],{ll,1,2},{dd,1,4},{mu,0,3},{nu,0,3}]+2 I Sum[z[LocalFieldNumber,5][ii,jj,kk,ll]γ5γdown[mu][[aa,bb]]γ5[[cc,dd]]D[Ψ[ll][dd][t,x,y,z],Ϛ[mu]],{ll,1,2},{mu,0,3},{dd,1,4}]+2I Sum[z[LocalFieldNumber,6][ii,jj,kk,ll]γ5γstdowndown[nu][[aa,bb]](γ5.CommuteGamma[nu,mu])[[cc,dd]]D[Ψ[ll][dd][t,x,y,z],Ϛ[mu]],{ll,1,2},{dd,1,4},{mu,0,3},{nu,0,3}]+2 I Sum[z[LocalFieldNumber,7][ii,jj,kk,ll]γ5down[[aa,bb]]γ5γ[mu][[cc,dd]]D[Ψ[ll][dd][t,x,y,z],Ϛ[mu]],{ll,1,2},{mu,0,3},{dd,1,4}]+2 I Sum[z[LocalFieldNumber,8][ii,jj,kk,ll]Cmetric[[aa,bb]]γ[mu][[cc,dd]]D[Ψ[ll][dd][t,x,y,z],Ϛ[mu]],{ll,1,2},{mu,0,3},{dd,1,4}]

Do[Z[LocalFieldNumber][kk,cc] = Table[ZTerms[LocalFieldNumber][ii,jj,kk][aa,bb,cc],{ii,1,2},{aa,1,4},{jj,1,2},{bb,1,4}],{kk,1,2},{cc,1,4}];

AlgebraDDResult[tt1,tt2][LocalFieldNumber]=Table[AlgebraDD[tt1,tt2][Ψh[kk][cc][t,x,y,z]]/.LagrangianParameters/.SetW,{kk,1,2},{cc,1,4}] ==Table[Z[LocalFieldNumber][kk,cc],{kk,1,2},{cc,1,4}]//Simplify

ZPrintout[LocalFieldNumber]= Table[(Subscript[(OverHat[(OverHat[Z])]), nn]^ijkl)== (cz[LocalFieldNumber,nn,18]((((σ^3)))^il)((((σ^("3"))))^jk) + cz[LocalFieldNumber,nn,22]((((σ^(HoldForm[1]))))^ij)((((σ^("1"))))^kl) + cz[LocalFieldNumber,nn,23]((((σ^(HoldForm[1]))))^ik)((((σ^("1"))))^jl)+cz[LocalFieldNumber,nn,24]((((σ^(HoldForm[1]))))^il)((((σ^("1"))))^jk)  )/.czSolution[LocalFieldNumber],{nn,1,NumZCharges}]
ZPrintoutTeXForm[LocalFieldNumber ] = TeXForm[%]

Clear[LocalFieldNumber,NumZCharges]

unhat-hat fermion closure summary

AlgebraDDResult[tt1,tt2][6]==AlgebraDDResult[tt1,tt2][16]==True

tt1
tt2