def1ch[1];

H1 = delta (number[d]-1) + U/2 pow[number[d]-1, 2] + B spinz[d];
Hc = gammaPol hop[f[0], d[]];

H = H0 + H1 + Hc;
