Dear all,

I have been doing some tests for the Bose Hubbard model in 1D with dirloop_sse and I find that the the output density and the value of the Green's function at (0,0) do not coincide, as I would expect.

I understand that when measuring Green's functions for the Bose-Hubbard model the output is that of <b^\dagger_0 b_j> with j going over all lattice points. b^\dagger_i is a creation operator and b_j is a destruction at sites i and j respectively. If that is correct, then, the following should be correct too and I should get <b^\dagger_0 b_0>=<n_0>=\rho. Where I denote the density by \rho.

When I do the calculation, using the following input:

<?xml version="1.0" encoding="UTF-8"?>
<?xml-stylesheet type="text/xsl" href="ALPS.xsl"?>
<SIMULATION xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xsi:noNamespaceSchemaLocation="http://xml.comp-phys.org/2003/8/QMCXML.xsd">
  <PARAMETERS>
<PARAMETER name="MEASURE[Green Function]">True</PARAMETER>
<PARAMETER name="MEASURE[Correlations]">True</PARAMETER>
<PARAMETER name="L">64</PARAMETER>
<PARAMETER name="LATTICE">chain lattice</PARAMETER>
<PARAMETER name="SEED">189185290</PARAMETER>
<PARAMETER name="U">10.0</PARAMETER>
<PARAMETER name="T">0.5</PARAMETER>
<PARAMETER name="MODEL">boson Hubbard</PARAMETER>
<PARAMETER name="Nmax">3</PARAMETER>
<PARAMETER name="THERMALIZATION">1000</PARAMETER>
<PARAMETER name="SWEEPS">50000</PARAMETER>
<PARAMETER name="NONLOCAL">0</PARAMETER>
<PARAMETER name="mu">0.4</PARAMETER>
<PARAMETER name="t">0.5</PARAMETER>
  </PARAMETERS>
</SIMULATION>

The output I get is:

<SCALAR_AVERAGE name="Density">
      <COUNT>55257</COUNT>
      <MEAN method="simple">0.6790808</MEAN>
      <ERROR converged="yes" method="binning">0.000306</ERROR>
      <VARIANCE method="simple">0.00258</VARIANCE>
      <AUTOCORR method="binning">0.505</AUTOCORR>

 <VECTOR_AVERAGE name="Green's Function" nvalues="64">
      <SCALAR_AVERAGE indexvalue="( 0 ) -- ( 0 )">
        <COUNT>55257</COUNT>
        <MEAN method="simple">1.178037</MEAN>
        <ERROR converged="yes" method="binning">0.000505</ERROR>
        <VARIANCE method="simple">0.0117</VARIANCE>
        <AUTOCORR method="binning">0.103</AUTOCORR>

Which is very different. Am I interpreting the output correctly or am I just not running it long enough? 

Also, what is the meaning of the last three fields in the xml output field? n,n^2 and n^4? Because my initial guess would be to interpret n as the number of particles but it is not consistent with the density and the lattice size.

Thank you very much,

Francisco Cordobes



.