Position Operator and Wannier Functions in Diamond
We plot the WF in diamond along the [111] direction (ie along a C-C bond)
and the corresponding position operator. There is essentially no difference
between the WF from a 4x4x4 grid and the WF from a 8x8x8 grid.