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