# gnuplot file # #
N=6 linear chain, Heisenberg model

set title "E0/s^2 vs. 1/s"
plot [0:] "-" u (2/$3):($5/($3/2)**2/6) t "N=6" w p,\
 "-" u (2/$3):($5/($3/2)**2/6) t "N=8 J2=0"   w p,\
 "-" u (2/$3):($5/($3/2)**2/6) t "N=8 J2=0.4" w p,\
-1-0.433*x,-1.333-0.56*x,-0.93-0.6*x
# N=6 Heisenberg chain with bigger spin s (2s=24 means s=12)
# even 2s
N=6 2s= 24 E0= -894.53099938 k=0+ 16c=3m, k=X fail_1GB n1=5375005 1c=8.1+46+Xm
N=6 2s= 12 E0= -231.31441734 n1=204763 1c=41s
N=6 2s= 10 E0= -162.77800492
N=6 2s= 8  E0= -106.24129442
N=6 2s= 6  E0= -61.70393148
N=6 2s= 4  E0= -29.16469320
N=6 2s= 2  E0= -8.61742318

# odd 2s
N=6 2s= 9  E0= -133.00970091
N=6 2s= 7  E0= -82.47273462
N=6 2s= 5  E0= -43.93471053
N=6 2s= 3  E0= -17.39284514
N=6 2s= 1  E0= -2.80277564
e
# temporary check N=8  (-1.6-0.66*x)*5/8=-1-0.412*(5./8)*x
N=8 2s= 8  E0= -140.87143178 k=0+ n1=73558 1c=62s k=X.n1=2306025 SH=30%*1GB 7h
N=8 2s= 6  E0= -81.69797814  k=0+ 1c=30s
N=8 2s= 4  E0= -38.52186914  k=0+ 1c=6s
N=8 2s= 2  E0= -11.33695608  k=0+
e
# J2=0.4 classical? (-1.09-0.78*x)*8./5
#N=8 2s=24 96+96 k=0+ ns_recursive=2h-timeout l1=1.8GB++
N=8 2s= 8  E0= -103.74060211 k=0+ J2=0.4
N=8 2s= 6  E0= -60.90637854  k=0+ J2=0.4
N=8 2s= 4  E0= -29.45266533  k=0+ J2=0.4
N=8 2s= 2  E0=  -9.31935449  k=0+ J2=0.4
e
pause -1
#