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>