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
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
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
Which version of ALPS are you using? On 11.02.2012, at 07:13, Chen Ahai <chenahaiphysics@gmail.com> wrote:
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
ALPS 1.3.4 thanks. chenahai
ok, thanks very much chenahai
participants (2)
-
Chen Ahai
-
Matthias Troyer