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