It seems that we must put a potential on boundary to obtain the ground state using DMRG.

How i can add a term as V/2*(n(1)+n(L)) to fermion model in models.xml

Thanks