Dear Prof. Matthias Troyer,

   Thank you for your quick response. I attach one of the results.  Since I find “Tau” is ‘nan’ (what is its meaning), I naively think the result may be incorrect. Maybe it is due to the choice of mu that there is no particle and I will check later.

  By the way, if I calculate many mu, how can I easily extract the results in the command way? Here the physical quantity is checkerboard charge order, can I get the charge order with other momentum?

Sincerely,

Huaiming

发件人: comp-phys-alps-users-bounces@lists.phys.ethz.ch [mailto:comp-phys-alps-users-bounces@lists.phys.ethz.ch] 代表 Matthias Troyer
发送时间: 2015年9月18日 12:28
收件人: comp-phys-alps-users@lists.phys.ethz.ch
主题: Re: [ALPS-users] dirloop_sse

Hi,

Can you explain what you mean by the result not being correct?

On Sep 18, 2015, at 05:06, Huaiming Guo <hmguo@buaa.edu.cn> wrote:

Dear all,

  I use dirloop_sse to calculate the hardcore-boson Hubbard model, but with next nearest neighbor hoppings. I revised the model and lattice as below or in the attached file. However the result is not correct. I want to know where the problem is. Thank you very much.

Sincerely,

Huaiming Guo

<HAMILTONIAN name="hardcore boson">

  <PARAMETER name="mu" default="0"/>

  <PARAMETER name="t" default="1"/>

  <PARAMETER name="V" default="0"/>

  <PARAMETER name="t'" default="0"/>

  <PARAMETER name="V'" default="0"/>

  <PARAMETER name="t''" default="0"/>

  <PARAMETER name="V''" default="0"/>

  <PARAMETER name="t0" default="t"/>

  <PARAMETER name="t1" default="t'"/>

  <PARAMETER name="V0" default="V"/>

  <PARAMETER name="V1" default="V'"/>

  <BASIS ref="hardcore boson"/>

  <SITETERM site="i">

    <PARAMETER name="mu#" default="mu"/>

    -mu#*n(i)

  </SITETERM>

  <BONDTERM source="i" target="j">

    <PARAMETER name="t#" default="0"/>

    <PARAMETER name="V#" default="0"/>

    -t#*(bdag(i)*b(j)+bdag(j)*b(i)) + V#*n(i)*n(j)

  </BONDTERM>

</HAMILTONIAN>

<UNITCELL name="nnn2d" dimension="2">

  <VERTEX/>

  <EDGE type="0"><SOURCE vertex="1" offset="0 0"/><TARGET vertex="1" offset="0 1"/></EDGE>

  <EDGE type="1"><SOURCE vertex="1" offset="0 0"/><TARGET vertex="1" offset="1 0"/></EDGE>

  <EDGE type="2"><SOURCE vertex="1" offset="0 0"/><TARGET vertex="1" offset="1 1"/></EDGE>

  <EDGE type="2"><SOURCE vertex="1" offset="0 0"/><TARGET vertex="1" offset="1 -1"/></EDGE>

</UNITCELL>

<LATTICEGRAPH name = "ladder" vt_name="LadderLattice">

  <FINITELATTICE>

    <LATTICE ref="square lattice"/>

    <PARAMETER name="W" default="2"/>

    <EXTENT dimension="1" size="L"/>

    <EXTENT dimension="2" size="W"/>

    <BOUNDARY dimension="1" type="periodic"/>

    <BOUNDARY dimension="2" type="open"/>

  </FINITELATTICE>

  <UNITCELL ref="nnn2d"/>

</LATTICEGRAPH>

<dd.dat>