

In [12]:

    
%pylab inline
from __future__ import division









    



Populating the interactive namespace from numpy and matplotlib

The following is the imput file 6 \ 1000 \ 5.E-1 \ 2.82691E-2 \ -0.25 \ 2.5 \ 0 \ 0 \ 0 \ 0 \ 9999999 \ 0.2 \ 0.0 \ 1.0 \ 1.0 \ 0.0 \ 0.0 \ 0.0 \ 0.0 \ Got OneChQS .. SetSingleSite done! Fsq = 2 chi_m1 = 0.121207 eN(0) = 0 sqrt(ALambda)*chi_m1 = 0.125319



In [13]:

    
Nsitesmax = 5
Ncutoff = 1000
U = 0.5
Gamma = 0.0282691
ed = -0.25
Lambda = 2.5
ChemPot = 0
UpdateBefCut = 0
calcdens = 0
Oliveira_z = 1



In [14]:

    
Gamma









    Out[14]:





0.0282691



In [15]:

    
x0=(1+Lambda**(-1))/2
x1=2/(1+Lambda**(-1))
x2=2/(1+Lambda**(-1))*(Lambda**(-1))
x3=0.5*U*2/(1+Lambda**(-1))*(Lambda**(-1))#this is the energy E0
x4=(ed+0.5*U)/(0.5*(1+Lambda**(-1))*Lambda)
x5=0.5*U*2/(1+Lambda**(-1))
x6=((x1**2)*2*Gamma/pi)**(1/2)
print(x0,x1,x2,x3,x4,x5,x6)
print x6*Lambda**(-1/2) ,x6*Lambda**(-1/2)*(1+Lambda**(-1))  #This is the one!!!!









    



(0.7, 1.4285714285714286, 0.5714285714285715, 0.14285714285714288, 0.0, 0.35714285714285715, 0.19164522924380614)
0.121207085423



In [28]:

    
b1 = 1-Lambda**(-Nsitesmax -1)
b2 = 1-Lambda**(-2*Nsitesmax -1)
b3 = 1-Lambda**(-2*Nsitesmax -3)
b4 =b1/(b2*b3)
print(b1,b2,b3,b4)









    



(0.995904, 0.99995805696, 0.9999932891136, 0.9959524567171915)



In [29]:

    
b1 = 1-Lambda**(-1-1)
b2 = 1-Lambda**(-2*1 -1)
b3 = 1-Lambda**(-2*1 -3)
b4 =b1/(b2*b3)
print(b1,b2,b3,b4)









    



(0.84, 0.9359999999999999, 0.98976, 0.9067207175839572)



In [30]:

    
b1 = 1-Lambda**(-1)
b2 = 1-Lambda**( -1)
b3 = 1-Lambda**(-3)
b4 =b1/(b2*b3)
print(b1,b2,b3,b4)









    



(0.6, 0.6, 0.9359999999999999, 1.0683760683760684)



In [15]:



In [15]:



In [ ]: