&GLOBAL
  PRINT_LEVEL low
  PROJECT LiH
  RUN_TYPE ENERGY_FORCE
  WALLTIME 600
&END GLOBAL

&MOTION
  &GEO_OPT
    MAX_DR 0.001
    MAX_FORCE 0.00015
    MAX_ITER 2
    OPTIMIZER BFGS
    RMS_DR 0.0005
    RMS_FORCE 0.0001
  &END GEO_OPT
  &PRINT
    &STRESS
    &END STRESS
    &STRUCTURE_DATA
      DISTANCE 1 5
    &END STRUCTURE_DATA
  &END PRINT
&END MOTION

&FORCE_EVAL
  METHOD Quickstep
  STRESS_TENSOR analytical
  #  STRESS_TENSOR numerical
  &DFT
    &POISSON
      &EWALD
        ALPHA 0.5
        EWALD_TYPE EWALD
        GMAX 41 41 41
        &MULTIPOLES
          MAX_MULTIPOLE_EXPANSION QUADRUPOLE
        &END MULTIPOLES
      &END EWALD
    &END POISSON
    &QS
      METHOD PM6
      &SE
        INTEGRAL_SCREENING KDSO-D
        PERIODIC EWALD
        &COULOMB
          CUTOFF 12.0
        &END COULOMB
      &END SE
    &END QS
    &SCF
      EPS_SCF 1.0E-6
      IGNORE_CONVERGENCE_FAILURE
      MAX_SCF 10
      SCF_GUESS MOPAC
      &OT
        LINESEARCH 3PNT
        MINIMIZER CG
        PRECONDITIONER NONE
      &END OT
    &END SCF
  &END DFT
  &PRINT
    &STRESS_TENSOR
    &END STRESS_TENSOR
  &END PRINT
  &SUBSYS
    &CELL
      ABC 4.796302 4.796302 4.796302
    &END CELL
    &COORD
      SCALED
      Li  0.000 0.000 0.000
      Li  0.500 0.500 0.000
      Li  0.500 0.000 0.500
      Li  0.000 0.500 0.500
      H   0.000 0.500 0.000
      H   0.000 0.000 0.500
      H   0.500 0.000 0.000
      H   0.500 0.500 0.500
    &END COORD
    &PRINT
      &CELL
      &END CELL
    &END PRINT
  &END SUBSYS
&END FORCE_EVAL
