ALPS users and developers, I am trying to calculate the phase diagram of extended boson hubbard using dmrg. The lattice is the finite one-dimension lattice with open boundaries. The chemical potentials are: u-=E(N)-E(N-1), u+=E(N+1)-E(N). E(N) is the ground energy for a total particle number N. However I found that there is something unusual when calculating the chemical potential: E(N+1)-E(N). When V is not zero, the u+ calculated by DMRG, sparsediag and mps_optim is less than the expected value. For example when L=100, N=50, t=0 and V=0.4, the correct value of u+ is 0.8, however using dmrg(mps and sparsediag) the value is 0.4 .
I can obtain the correct value using u+=E(N+2)-E(N+1) when choosing L=41 and N=21 instead. However within L=101 and N=101, even choose u+=E(N+2)-E(N+1) the result is still 0.4. It seems that when V does not equal to zero the calculated chemical potential is related to specific total particle numbers.
Xiao Liang