vasp_tools/e_fc_tersoff

   1 #!/bin/bash
   2
   3 ready=`grep General\ timing $1/OUTCAR | awk '{ print $1 }'`
   4
   5 if [ "$ready" = "General" ]; then
   6         echo "$i seems to be finisched ..."
   7
   8 e=`./e_form_tersoff plain_si_2333_nosym_sp plain_sic_2333_nosym_sp $1 | grep M2 | \
   9         awk '{ print $4 }'`
  10 echo "E_f = $e"
  11
  12 ec100=3.72027
  13 r=`echo $e $ec100 $ec100 | awk '{ print $1-$2-$3 }'`
  14 echo "C 100: $r"
  15
  16 ec110=4.16327
  17 r=`echo $e $ec100 $ec110 | awk '{ print $1-$2-$3 }'`
  18 echo "C 110: $r"
  19
  20 ecsub=1.9548
  21 r=`echo $e $ec100 $ecsub | awk '{ print $1-$2-$3 }'`
  22 echo "C sub: $r"
  23
  24 ecbc=4.66227
  25 r=`echo $e $ec100 $ecbc | awk '{ print $1-$2-$3 }'`
  26 echo "C bc: $r"
  27
  28 esvac=3.62905
  29 r=`echo $e $ec100 $esvac | awk '{ print $1-$2-$3 }'`
  30 echo "Si vac: $r"
  31
  32 es110=3.38999
  33 r=`echo $e $ec100 $es110 | awk '{ print $1-$2-$3 }'`
  34 echo "Si 110: $r"
  35
  36 estet=3.76899
  37 r=`echo $e $ec100 $estet | awk '{ print $1-$2-$3 }'`
  38 echo "Si tet: $r"
  39
  40 eshex=3.41999
  41 r=`echo $e $ec100 $eshex | awk '{ print $1-$2-$3 }'`
  42 echo "Si hex: $r"
  43
  44 r=`echo $e $es110 $ecsub | awk '{ print $1-$2-$3 }'`
  45 echo "Si 110 and C sub: $r"
  46
  47 r=`echo $e $estet $ecsub | awk '{ print $1-$2-$3 }'`
  48 echo "Si tet and C sub: $r"
  49
  50 else
  51         echo "$i not yet finished ..."
  52 fi
  53