Hi Chen, Ahai,

You don't need to fix a temperature when you run fulldiag but should specify the temperature only later when running fulldiag_evaluate. Have you tried that?

Matthias Troyer

On 11 Feb 2012, at 05:40, Chen Ahai wrote:

Dear  all ALPS users 

   I want to calculate the density distribution of hubbard model with 
   external potential in finite temperature using the "fulldiag" command, 
   the main parameters are
   ------------------------------------------------------------------------
   .........
   t=1.0
   CONSERVED_QUANTUMNUMBERS="Sz,Nup,Ndown,N"
   MEASURE_LOCAL[n_density]=n
   MEASURE_LOCAL[n_up]=n_up
   MEASURE_LOCAL[n_down]=n_down
   {T=1.0}
   {T=10.0}
   ------------------------------------------------------------------------

   It can get different thermodynamic quantities (e.g. free energy, specific heat) 
   in different T, but the quantities I defined above keep the same in the two tasks.
   
   I can't find  the problem in this test. 
   Any help would be appreciated....
   Thanks.


   Chen Ahai
   Department of Physics, Zhejiang Normal University
   Jinhua, Zhejiang Province, 321004, China
   e-mail: chenahaiphysics@gmail.com;
              chenahai@yahoo.cn