set term post eps enh color 12 
set out "Hubbard_Heisenberg_NNN_v2.eps"
set size 0.5,0.5
set multiplot
set samples 1000
set rmargin 1
set lmargin 6

PS=0.8
PS1=0.7
PS2=0.6
set style line 2 lt 2 lc rgb "#00BB00" lw 3 pt 4  ps 0.5
set style line 5 lt 5 lc rgb "#00BBBB" lw 3
set style line 7 lt 1 lc rgb "#FFFF7E" lw 15
set style line 8 lt 1 lc rgb "#00FFFF" lw 15

set xlabel "s" offset 0,0.5
set ylabel "<S_i {/Symbol \327} S_{i+{/Symbol d}}>" offset 2
set key samplen 2 spacing 1.3 left Left reverse at screen 0.07,0.23
# reverse
set xtics 0.2
set mxtics 2
set ytics 0.2
set mytics 2


set arrow 1 from log(2),-0.05 to log(2),0.2 nohead lt 4 lc -1 front
#set label 1 "{/Symbol d} = (1,0)" at 0.1,-0.4
#set label 2 "{/Symbol d} = (1,1)" at 0.1,0.3
set label 2 "{/=14{/Symbol d} = NNN}" at screen 0.16,0.44

set label 3 "{/=9 Weiss MF}" at 0.1,0.96 rotate by -34

set xzeroaxis

plot [0:1.05][-0.05:1]\
"Data/Heisenberg/Heisenberg_WeissMF_T_E_s.dat" u 3:($3<0.1?2*$2*0.905:1/0) t"" w l lt 1,\
x<0.6?0.76-1.48*x**1.5  + 0.32*x**3:1/0  t"" w l lt 2 lc rgb "#00BB00",\
x<0.3? 0.63 -1.77 *x**2+x**4*1.1 :1/0 t"" lt 4 lw 1,\
"Data/DMFT/results_U150mu0_ent_corr.dt0" u 18:($7**2) t"DMFT" w lp lt 1 lw 3 pt 6 ps PS1,\
"Data/DQMC_3D/sisj110u15.dat" u 2:($2>0.2?3*$3:1/0) w lp t"d = 3" lt 3 lw 3 pt 12 ps PS1,\
"Data/DQMC_2D/U12.25L6_beta_s_sisj.dat" u ($2):($2>0.2?3*$4:1/0) w lp t "d = 2" ls 2,\
"< paste Data/DQMC_1D/DQMC_1d_U867_N16_32_64.av Data/DQMC_1D/BA_DQMCgrid_beta_s_E_D.dat" u 16:(3*$14) t"d = 1" w lp lt 4 lw 3 pt 1 ps PS1,\
"Data/Heisenberg/Heisenberg_WeissMF_T_E_s.dat" u 3:(-2*$2) t"" w l lt -1,\
"Data/Heisenberg/Heisenberg_16x16x16_RTS.dat" u 9:(12*$7) t"" w l lt 1 lw 2 lc 3,\
x<0.1?0.857-2.16*3/4*x**1.333:1/0 t"" w l lt 1 lw 2 lc 3,\
x<0.05?0.812-2.34*2/3*x**1.5:1/0 t"" w l lt 1 lw 2 lc 2,\
"Data/Heisenberg/Heisenberg_Scalettar_12x12.dat" u 9:($9>0.05?12*$7:1/0) t"" w l lt 1 lw 2 lc 2,\
"Data/Heisenberg/corr_xxx_AF_s.dat" u 10:(12*$3) t"" w l lt 1 lw 2 lc 4

set size 0.25,0.3
set origin 0.235,0.18
set xlabel ""
set ylabel ""
set format "{/=10 %g}"
set xtics 0.1
set ytics 0.1
set arrow 1 from log(2),-0.01 to log(2),0.1 nohead lt 4 lc -1 front


plot [0.55:0.75][-0.01:0.25]\
"Data/Heisenberg/Heisenberg_WeissMF_T_E_s.dat" u 3:($3<0.1?2*$2*0.905:1/0) t"" w l lt 1,\
x<0.6?0.76-1.48*x**1.5  + 0.32*x**3:1/0  t"" w l lt 2 lc rgb "#00BB00",\
x<0.3? 0.63 -1.77 *x**2+x**4*1.1 :1/0 t"" lt 4 lw 1,\
"Data/DMFT/results_U150mu0_ent_corr.dt0" u 18:($7**2) t"DMFT" w lp lt 1 lw 3 pt 6 ps PS1,\
"Data/DQMC_3D/sisj110u15.dat" u 2:($2>0.2?3*$3:1/0) w lp t"d = 3" lt 3 lw 3 pt 12 ps PS1,\
"Data/DQMC_2D/U12.25L6_beta_s_sisj.dat" u ($2):($2>0.2?3*$4:1/0) w lp t "d = 2" ls 2,\
"< paste Data/DQMC_1D/DQMC_1d_U867_N16_32_64.av Data/DQMC_1D/BA_DQMCgrid_beta_s_E_D.dat" u 16:(3*$14) t"d = 1" w lp lt 4 lw 3 pt 1 ps PS1,\
"Data/Heisenberg/Heisenberg_WeissMF_T_E_s.dat" u 3:(-2*$2) t"" w l lt -1,\
"Data/Heisenberg/Heisenberg_16x16x16_RTS.dat" u 9:(12*$7) t"" w l lt 1 lw 2 lc 3,\
x<0.05?0.857-x**1.333:1/0 t"" w l lt 1 lw 2 lc 3,\
x<0.05?0.812-x**1.5:1/0 t"" w l lt 1 lw 2 lc 2,\
"Data/Heisenberg/Heisenberg_Scalettar_12x12.dat" u 9:($9>0.05?12*$7:1/0) t"" w l lt 1 lw 2 lc 2,\
"Data/Heisenberg/corr_xxx_AF_s.dat" u 10:(12*$3) t"" w l lt 1 lw 2 lc 4


