#!/bin/sh

echo "Geometry optimization"

$ADFBIN/band << eor
Title Graphene geometry optimization

Units
  Length Angstrom
end

! For Graphene we need to use a symmetryc grid
KSpace 5 &
   Grid Symmetric
End

ATOMS
  C   0.0 0.0 0.0
  C   1.23 0.7101408312 0.0
END

Lattice
   2.46 0.000000 0
   1.23 2.130422493 0
End

GeoOpt
  OptimizeLattice true
  Converge grad=1e-5
End

BasisDefaults
  BasisType DZ
end

end input
eor

mv RUNKF Graphene.runkf


echo "Phonon calculation"


$ADFBIN/band << eor
Title Graphene phonon calc

Units
  Length Angstrom
end

KSpace 3 &
   Grid Symmetric
End

Restart
  File Graphene.runkf
  Geometry
End

runtype
  phonons
end

phononConfig
  stepSize 0.0913
  superCell
    2 0
    0 2
  subend
end

ATOMS
  C   0.0 0.0 0.0
  C   1.23 0.7101408312 0.0
END

Lattice
   2.46 0.000000 0
   1.23 2.130422493 0
End

BasisDefaults
  BasisType DZ
end

end input
eor

NSCM=1
export NSCM
echo ""
echo "Begin TOC"

$ADFBIN/dmpkf -n 1 RUNKF --toc

echo "End TOC"