O atom binding energy

set term postscript eps color enhanced 14

set output "BGCb_dEO.eps"

set size 0.5,0.5

set size ratio 0.8

set xlabel "{/Times-Roman BGC (eV)}"

xmin=-7.1 ; xmax=-1.0

set xr[xmin:xmax]

set xtics -9,1,0 format "%.0f"

set mxtics 2

set xtics nomirror

set ylabel "{/Symbol D}{/Times-Itallic E}{/Times-Roman _{O} (eV)}"

ymin=1.6 ; ymax = 4

set yr[ymin:ymax]

set ytics 1.,0.5,5 format "%.1f"

set mytics 5

set key box lw 0.5

set key at -1.0,ymax-0.0  font "Helvetica,10"

set label at -3.45, 2.45 tc rgb "red" "{/=8 1{/Symbol \264}1 SrO}"        

set label at -4.50, 3.08 tc rgb "red" "{/=8 1{/Symbol \264}1 TiO_2}"      

set label at -5.63, 3.32 tc rgb "red" "{/=8 c(4{/Symbol \264}2) TiO_2 DL}"

set label at -5.96, 3.18 tc rgb "red" "{/=8 c(2{/Symbol \264}2) TiO_2 DL}"

set label at -6.60, 2.89 tc rgb "red" "{/=8 2{/Symbol \264}1 TiO_2 DL}"   

set label at -3.10, 2.65 tc rgb "black" "{/=8 c(2{/Symbol \264}2) BaO}"

set label at -3.55, 1.92 tc rgb "black" "{/=8 1{/Symbol \264}1 BaO DL}"

set label at -4.89, 2.90 tc rgb "black" "{/=8 c(2{/Symbol \264}2) TiO_2}"

set label at -6.70, 3.07 tc rgb "black" "{/=8 2{/Symbol \264}1 TiO_2 DL}"

set label at -5.90, 3.90 tc rgb "blue" "{/=8 1{/Symbol \264}1 TiO_2}"

set arrow from -1.5,2.2653945 to -6.8,3.3224424  nohead lt 1 lc 7 lw 2

set arrow from xmin,2.3 to xmax,2.3 nohead lt 5 lw 1 lc 7  

p 'Eads_sto.dat'    u ($8+$9)/2:10  w p pt 5  ps 1.0 lc 1 t 'STO' ,\

  'Eads_bto.dat'    u ($8+$9)/2:10  w p pt 7  ps 1.0 lc 7 t 'BTO' ,\

  'Eads_tio2.dat'   u ($8+$9)/2:10  w p pt 13 ps 1.0 lc 3 t 'Rutile'

##  2.3 w l lt 5 lw 2 lc 9 t ''  

##  f(x) lt 1 lc 7 lw 2 t ''

# f(x) = a*x+b

# a=-1.0 ; b= 4.0

#

# set term postscript eps color enhanced

# set out "BGCb_fit.eps"

#

# fit f(x) 'Eads_BGCb.dat' u ($8+$9)/2:10 via a,b

# plot 'Eads_BGCb.dat' u ($8+$9)/2:10 w p pt 1 ps 1 ,\

#       f(x) lt 1 lw 2

#

# ## a               = -0.199443        +/- 0.06447      (32.32%)

# ## b               = 1.96623          +/- 0.3297       (16.77%)