new imp_solver class done, need to adapt mps_iface.py
bug fixes
now interface to MPS_DMRG solver also working
removed unused code
several bugs related to mean-field calculation
exact nqp is 4*nsites