Dear Matthias

   thanks for your answer.

   I find the problem in my parameter file that I forgot to remove  "NUMBER_EIGENVALUES=1".

   So I only get the same ground state.

   

   I have run fulldiag_evaluate with 

   "fulldiag_evaluate --T_MIN 1 --T_MAX 1 --DELTA_T 0.01 1.task1.out.xml"

   and got the  following files

   ----------------------------------------------------

   1.task1.plot.energy.xml

   1.task1.plot.free_energy.xml 

   1.task1.plot.entropy.xml

   .....

   -----------------------------------------------------------------

   but I need the density distribution of the system,

   how can I get this result ... or it must be done by myself using all the eigenstates?

   thanks.

chenahai