On 16 Mar 2011, at 10:23, Miguel Dias Costa wrote:
Hello all.
As a test case for Alps 2, we are running fulldiag for a 1D Hubbard model with only three lattice sites,
res = pyalps.runApplication('fulldiag',input_file,writexml=True) data = pyalps.loadEigenstateMeasurements(pyalps.getResultFiles(prefix='parm1'))
and are measuring
'MEASURE_CORRELATIONS[Sx Sx]' : "Sx:Sx", 'MEASURE_CORRELATIONS[Sz Sz]' : "Sz:Sz", 'MEASURE_CORRELATIONS[NupNdw]' : "Nup:Ndown", 'MEASURE_CORRELATIONS[NupNup]' : "Nup:Nup",
and everything seems to check out fine, except of <Sx Sx>. We see that the site operator Sx is correctly defined from Splus and Sminus, but is this the way to use it in the measurements?
Cheers, Miguel
You will need to define SPlus and Sminus for the Hubbard model. Dod you do that?
Matthias