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>