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
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
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
Hi Matthias,
Don't worry, I believe the problem has been solved. It looks like a problem with the parameters I was using. Thanks for your help anyway!
All the best, Kristina
On 3 June 2015 at 16:08, Kristina Fidanovski 21292685@student.uwa.edu.au wrote:
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
comp-phys-alps-users@lists.phys.ethz.ch