We are using some of our own lattices and models, but the problem we are encountering can be identified using existing lattices ("inhomogeneous square lattice" and "inhomogeneous simple cubic lattice") and different fermion Hubbard models.
As one concrete example, with the following (simplified) parameter file, only for some of the runs does one get local density data. Two example outputs are shown. Note that the same thing occurs (of course) if we use two separate parameter files, each one listing a particular disorder seed. Seeds 11113, 11114 work, 11115 does not, etc. Also, the problem exists for N_total=6 and 4 ...
So, do constraints exist on valid choices of DISORDERSEED? (Also, since we would like to know the complexion of disorder associated with a given ground state wave function, we would in fact like to know the random number generator being used.)
Thanks. ----
LATTICE_LIBRARY = "/usr/local/ALPS/lib/xml/lattices.xml"; LATTICE = "inhomogeneous simple cubic lattice"; MODEL_LIBRARY= "/usr/local/ALPS/lib/xml/models.xml"; MODEL = "alternative fermion Hubbard"; CONSERVED_QUANTUMNUMBERS="Sz_total"; CONSERVED_QUANTUMNUMBERS="N_total"; MEASURE_LOCAL[Local charge density]=n; t = 1.0; U = 2.0; L = 2; Sz_total = 0; N_total = 8; TRANSLATION_SYMMETRY = false; { DISORDERSEED = 11111; mu=-6.0*(random()-0.5); } { DISORDERSEED = 11112; mu=-6.0*(random()-0.5); }
some of task1 output:
<EIGENVALUES number="1"> -13.57576898386332
</EIGENVALUES> <EIGENSTATES number="1"> <EIGENSTATE number="0"> <SCALAR_AVERAGE name="Energy"><MEAN>-13.57576898386332</MEAN></SCALAR_AVERAGE> <VECTOR_AVERAGE name="Local charge density"> <SCALAR_AVERAGE indexvalue="( 0,0,0 )"><MEAN>1.02945417967653</MEAN></SCALAR_AVERAGE> <SCALAR_AVERAGE indexvalue="( 0,0,1 )"><MEAN>0.9407378693270886</MEAN></SCALAR_AVERAGE> <SCALAR_AVERAGE indexvalue="( 0,1,0 )"><MEAN>0.8692266888738961</MEAN></SCALAR_AVERAGE> <SCALAR_AVERAGE indexvalue="( 0,1,1 )"><MEAN>0.2243604886901401</MEAN></SCALAR_AVERAGE> <SCALAR_AVERAGE indexvalue="( 1,0,0 )"><MEAN>0.9836579726343326</MEAN></SCALAR_AVERAGE> <SCALAR_AVERAGE indexvalue="( 1,0,1 )"><MEAN>0.9533153322138351</MEAN></SCALAR_AVERAGE> <SCALAR_AVERAGE indexvalue="( 1,1,0 )"><MEAN>1.256609332250588</MEAN></SCALAR_AVERAGE> <SCALAR_AVERAGE indexvalue="( 1,1,1 )"><MEAN>1.742638136334122</MEAN></SCALAR_AVERAGE> </VECTOR_AVERAGE> </EIGENSTATE> </EIGENSTATES>
some of task2 output:
<EIGENVALUES number="1"> -12.3429442715644
</EIGENVALUES> <EIGENSTATES number="1"> <EIGENSTATE number="0"> <SCALAR_AVERAGE name="Energy"><MEAN>-12.3429442715644</MEAN></SCALAR_AVERAGE> <VECTOR_AVERAGE name="Local charge density"> <SCALAR_AVERAGE indexvalue="( 0,0,0 )"><MEAN>0</MEAN></SCALAR_AVERAGE> <SCALAR_AVERAGE indexvalue="( 0,0,1 )"><MEAN>0</MEAN></SCALAR_AVERAGE> <SCALAR_AVERAGE indexvalue="( 0,1,0 )"><MEAN>0</MEAN></SCALAR_AVERAGE> <SCALAR_AVERAGE indexvalue="( 0,1,1 )"><MEAN>0</MEAN></SCALAR_AVERAGE> <SCALAR_AVERAGE indexvalue="( 1,0,0 )"><MEAN>0</MEAN></SCALAR_AVERAGE> <SCALAR_AVERAGE indexvalue="( 1,0,1 )"><MEAN>0</MEAN></SCALAR_AVERAGE> <SCALAR_AVERAGE indexvalue="( 1,1,0 )"><MEAN>0</MEAN></SCALAR_AVERAGE> <SCALAR_AVERAGE indexvalue="( 1,1,1 )"><MEAN>0</MEAN></SCALAR_AVERAGE> </VECTOR_AVERAGE> </EIGENSTATE> </EIGENSTATES>
Could it be just that the values of mu are such that the ground state has no particles?
Matthias
On Aug 19, 2009, at 1:23 PM, Bob Gooding wrote:
We are using some of our own lattices and models, but the problem we are encountering can be identified using existing lattices ("inhomogeneous square lattice" and "inhomogeneous simple cubic lattice") and different fermion Hubbard models.
As one concrete example, with the following (simplified) parameter file, only for some of the runs does one get local density data. Two example outputs are shown. Note that the same thing occurs (of course) if we use two separate parameter files, each one listing a particular disorder seed. Seeds 11113, 11114 work, 11115 does not, etc. Also, the problem exists for N_total=6 and 4 ...
So, do constraints exist on valid choices of DISORDERSEED? (Also, since we would like to know the complexion of disorder associated with a given ground state wave function, we would in fact like to know the random number generator being used.)
Thanks.
LATTICE_LIBRARY = "/usr/local/ALPS/lib/xml/lattices.xml"; LATTICE = "inhomogeneous simple cubic lattice"; MODEL_LIBRARY= "/usr/local/ALPS/lib/xml/models.xml"; MODEL = "alternative fermion Hubbard"; CONSERVED_QUANTUMNUMBERS="Sz_total"; CONSERVED_QUANTUMNUMBERS="N_total"; MEASURE_LOCAL[Local charge density]=n; t = 1.0; U = 2.0; L = 2; Sz_total = 0; N_total = 8; TRANSLATION_SYMMETRY = false; { DISORDERSEED = 11111; mu=-6.0*(random()-0.5); } { DISORDERSEED = 11112; mu=-6.0*(random()-0.5); }
some of task1 output:
<EIGENVALUES number="1"> -13.57576898386332
</EIGENVALUES> <EIGENSTATES number="1"> <EIGENSTATE number="0"> <SCALAR_AVERAGE name="Energy"><MEAN>-13.57576898386332</MEAN></SCALAR_AVERAGE> <VECTOR_AVERAGE name="Local charge density"> <SCALAR_AVERAGE indexvalue="( 0,0,0 )"><MEAN>1.02945417967653</MEAN></SCALAR_AVERAGE> <SCALAR_AVERAGE indexvalue="( 0,0,1 )"><MEAN>0.9407378693270886</MEAN></SCALAR_AVERAGE> <SCALAR_AVERAGE indexvalue="( 0,1,0 )"><MEAN>0.8692266888738961</MEAN></SCALAR_AVERAGE> <SCALAR_AVERAGE indexvalue="( 0,1,1 )"><MEAN>0.2243604886901401</MEAN></SCALAR_AVERAGE> <SCALAR_AVERAGE indexvalue="( 1,0,0 )"><MEAN>0.9836579726343326</MEAN></SCALAR_AVERAGE> <SCALAR_AVERAGE indexvalue="( 1,0,1 )"><MEAN>0.9533153322138351</MEAN></SCALAR_AVERAGE> <SCALAR_AVERAGE indexvalue="( 1,1,0 )"><MEAN>1.256609332250588</MEAN></SCALAR_AVERAGE> <SCALAR_AVERAGE indexvalue="( 1,1,1 )"><MEAN>1.742638136334122</MEAN></SCALAR_AVERAGE> </VECTOR_AVERAGE> </EIGENSTATE> </EIGENSTATES>
some of task2 output:
<EIGENVALUES number="1"> -12.3429442715644
</EIGENVALUES> <EIGENSTATES number="1"> <EIGENSTATE number="0"> <SCALAR_AVERAGE name="Energy"><MEAN>-12.3429442715644</MEAN></SCALAR_AVERAGE> <VECTOR_AVERAGE name="Local charge density"> <SCALAR_AVERAGE indexvalue="( 0,0,0 )"><MEAN>0</MEAN></ SCALAR_AVERAGE> <SCALAR_AVERAGE indexvalue="( 0,0,1 )"><MEAN>0</MEAN></ SCALAR_AVERAGE> <SCALAR_AVERAGE indexvalue="( 0,1,0 )"><MEAN>0</MEAN></ SCALAR_AVERAGE> <SCALAR_AVERAGE indexvalue="( 0,1,1 )"><MEAN>0</MEAN></ SCALAR_AVERAGE> <SCALAR_AVERAGE indexvalue="( 1,0,0 )"><MEAN>0</MEAN></ SCALAR_AVERAGE> <SCALAR_AVERAGE indexvalue="( 1,0,1 )"><MEAN>0</MEAN></ SCALAR_AVERAGE> <SCALAR_AVERAGE indexvalue="( 1,1,0 )"><MEAN>0</MEAN></ SCALAR_AVERAGE> <SCALAR_AVERAGE indexvalue="( 1,1,1 )"><MEAN>0</MEAN></ SCALAR_AVERAGE> </VECTOR_AVERAGE> </EIGENSTATE> </EIGENSTATES>
comp-phys-alps-users@lists.phys.ethz.ch