def2ch[2];

H1 = epsilon number[d[]] + U hubbard[d[]];
H2 = epsilon number[a[]] + U hubbard[a[]];
H12 = t12 hop[d[], a[]];
Himp = H1 + H2 + H12;

Hc = gammaPolCh[1] hop[f[0], d[]] + gammaPolCh[2] hop[f[1], a[]];

H = H0 + Hc + Himp;

Hselfd = Himp;
selfopd = ( Chop @ Expand @ komutator[Hselfd /. params, d[#1, #2]] )&;

Hselfa = Himp;
selfopa = ( Chop @ Expand @ komutator[Hselfa /. params, a[#1, #2]] )&;

