Dear developers,
I'm using the dirloop_sse code to simulate the hardcore boson Hubbard model, I met some problems there:
(1) An error occors in many simulations, that is, Internal error: size of prob table is zero.
(2) In a ladder model, the stiffness I get from the sse code is much smaller (almost zero) than using other codes (not ALPS code). While in a square lattice, the sse code gives the same stiffness with other codes; BTW, does the definition of stiffness in the ALPS codes miss a factor of 2 in the denominator?
(3) When I input some values of chemical potential (mu), I can get the corresponding particle density (rho). When mu is close to zero, the data of rho are not good, while for other mu, the data are nice.
PS: I attach my parameter file below, do I make something wrong in the parameters?
LATTICE="ladder";
L=64;
MODEL="hardcore boson";
V = 0.0;
V' = 0.0;
t = 1.0;
t' = 1.0;
T = 0.03125;
SWEEPS=10000000;
THERMALIZATION=1000000;
{ mu=-0.1; }
Best,
TaoYing
--
Tao Ying,
Ph.D candidate,
Physics department
Harbin Institute of Technology