$rij
$rundimensions
   natoms=2
   nbf(CAO)=38
   nbf(AO)=36
$energy    file=energy
$scfdump
$ricore    500
$drvopt
cartesian  on
   basis      off
   global     off
   hessian    on
   dipole     on
   nuclear polarizability
$forceupdate
ahlrichs numgeo=0  mingeo=3 maxgeo=4 modus=<g|dq> dynamic fail=0.3
   threig=0.005  reseig=0.005  thrbig=3.0  scale=1.00  damping=0.0
$dft
 functional pbe
 gridsize m5
$subtitle
SP,DFT-PBC/RI,charge 0,c1,pob-DZVP-rev2,
$periodic    3
$cell
  7.3335901515 7.3343969161 7.3338845743 59.9889124445 59.9990287166 60.0003299374
$scfiterlimit    100
$coord    file=coord
$jbas    file=auxbasis
$thime    5
$scfconv    8
$forceinit on
   diag=default
$scforbitalshift    automatic=.1
$grad    file=gradient
$interconversion    off
   qconv=1.d-7
   maxiter=25
$optimize
 internal off
 redundant off
 cartesian on
 global off
$symmetry    c1
$tmole
$kpoints
  nkpoints 9 9 9
  kptlines 10
  recipr  0  0  0  0.5  0.0  0.5  10
  recipr  0.5  0.0  0.5  0.5  0.25  0.75  10
  recipr  0.5  0.25  0.75  0.375  0.375  0.75  10
  recipr  0.375  0.375  0.75  0  0  0  10
  recipr  0  0  0  0.5  0.5  0.5  10
  recipr  0.5  0.5  0.5  0.625  0.25  0.625  10
  recipr  0.625  0.25  0.625  0.5  0.25  0.75  10
  recipr  0.5  0.25  0.75  0.5  0.5  0.5  10
  recipr  0.5  0.5  0.5  0.375  0.375  0.75  10
  recipr  0.625  0.25  0.625  0.5  0.0  0.5  10
$basis    file=basis
$scfmo    file=mos
$title    bandstructure
$atoms
si 1-2 \
   jbas = si universal                     \
   basis = si pob-DZVP-rev2
$marij
$coordinateupdate
dqmax=0.3
   interpolate  on
   statistics    5
$scfdamp     start=0.700 step=0.050 min=0.100
$operating system    unix
$riper
 pqmatdiag off
 lenonly on
$thize    0.10000000E-04
$scfintunit
unit=30       size=0        file=twoint
$forceapprox    file=forceapprox
$scfdiis
$xctype GGA (not to be modified here)
$last SCF energy change = -578.59249
$closed shells
 a       1-14                                   ( 2 )
$maxcor    500 MiB  per_core
$last step     riper
$end
