“…We shall develop both the one-(ǫ, t) and two-body(X, U, V, J z = J xy = P ) electron interaction parameters, evaluated as in Ref. 1, by assuming a non-degenerate orbital described by Wannier functions built from atomic orbitals of Gaussian shape. We associate to each site a Gaussian atomic-like orbital φ i (r − R i ), with the ions centered at the positions R i ≡ (±a/2 + u i , 0, 0) (i = 1, 2) .…”