Hi Matthias,

I've been using python scripts such as:

import pyalps
import matplotlib.pyplot as plt
import pyalps.plot

#prepare the input parameters
parms = []
for t in [ 0.02, 0.04, 0.06, 0.1]:
    parms.append(
        { 
          'LATTICE'        : "honeycomb lattice", 
          'MODEL'          : "boson Hubbard",
          'T'              : 0.1,
          'L'              : 24 ,
          't'              : t ,
          'mu'             : 0.75,
          'U'              : 1.0,
          'NONLOCAL'       : 0 ,
          'Nmax'           : 2 ,
          'THERMALIZATION' : 10000,
          'SWEEPS'         : 500000
        }
    )

#write the input file and run the simulation
input_file = pyalps.writeInputFiles('parm5a',parms)
res = pyalps.runApplication('worm',input_file,Tmin=5)

#load the magnetization and collect it as function of field h
data = pyalps.loadMeasurements(pyalps.getResultFiles(prefix='parm5a'),'Stiffness')
rhos = pyalps.collectXY(data,x='t',y='Stiffness')

#make plot
plt.figure()
pyalps.plot.plot(rhos)
plt.xlabel('Hopping $t/U$')
plt.ylabel('Superfluid density $\\rho _s$')
plt.show()

With various combinations of parameters and the inclusion of V (in the hopes of at least reproducing results obtained in other studies).

Kind Regards,
Kristina

On 3 June 2015 at 15:16, Matthias Troyer <troyer@phys.ethz.ch> wrote:
Hi Kristina,

what code are you using and what are your input files? We cannot help without that information.

Matthias Troyer




> On Jun 2, 2015, at 21:07, Kristina Fidanovski <21292685@student.uwa.edu.au> wrote:
>
> Dear Developers,
>
> I'm trying to look into the extended boson Hubbard model on a honeycomb lattice, but no matter what I do the density always comes out as zero. I don't seem to have the same problem when using a honeycomb lattice with any of the other examples in the tutorials. What could be the problem?
>
> All the best,
> Kristina