&GLOBAL
  PRINT_LEVEL LOW
  PROJECT test
  RUN_TYPE DEBUG
&END GLOBAL

&DEBUG
  DEBUG_FORCES T
  DEBUG_STRESS_TENSOR F
  EPS_NO_ERROR_CHECK 5.0E-4
&END DEBUG

&FORCE_EVAL
  METHOD QMMM
  &DFT
    &QS
      METHOD DFTB
      &DFTB
        DISPERSION F
        DO_EWALD F
        SELF_CONSISTENT T
        &PARAMETER
          PARAM_FILE_NAME scc_parameter
          PARAM_FILE_PATH DFTB/scc
        &END PARAMETER
      &END DFTB
    &END QS
    &SCF
      EPS_SCF 1.E-10
      SCF_GUESS MOPAC
    &END SCF
  &END DFT
  &MM
    &FORCEFIELD
      &BEND
        ATOMS H O H
        K 0.
        THETA0 1.8
      &END BEND
      &BOND
        ATOMS O H
        K 0.
        R0 1.8
      &END BOND
      &CHARGE
        ATOM O
        CHARGE -0.8476
      &END CHARGE
      &CHARGE
        ATOM H
        CHARGE 0.4238
      &END CHARGE
      &NONBONDED
        &LENNARD-JONES
          ATOMS O O
          EPSILON 78.198
          RCUT 11.4
          SIGMA 3.166
        &END LENNARD-JONES
        &LENNARD-JONES
          ATOMS O H
          EPSILON 0.0
          RCUT 11.4
          SIGMA 3.6705
        &END LENNARD-JONES
        &LENNARD-JONES
          ATOMS H H
          EPSILON 0.0
          RCUT 11.4
          SIGMA 3.30523
        &END LENNARD-JONES
      &END NONBONDED
    &END FORCEFIELD
    &POISSON
      &EWALD
        ALPHA .44
        EWALD_TYPE SPME
        GMAX 42
      &END EWALD
    &END POISSON
  &END MM
  &QMMM
    CENTER SETUP_ONLY
    ECOUPL POINT_CHARGE
    NOCOMPATIBILITY
    &CELL
      ABC 8.0 8.0 8.0
    &END CELL
    &MM_KIND H
      RADIUS 0.44
    &END MM_KIND
    &MM_KIND O
      RADIUS 0.78
    &END MM_KIND
    &PERIODIC
      &POISSON
        &EWALD
          ALPHA .44
          EWALD_TYPE SPME
          GMAX 21
        &END EWALD
      &END POISSON
    &END PERIODIC
    &QM_KIND H
      MM_INDEX 2 3
    &END QM_KIND
    &QM_KIND O
      MM_INDEX 1
    &END QM_KIND
  &END QMMM
  &SUBSYS
    &CELL
      ABC 24.955 24.955 24.955
    &END CELL
    &COORD
      O     0.000000     0.000000     0.000000    H2O1
      H     0.000000     0.000000     1.000000    H2O1
      H     0.942809     0.000000    -0.333333    H2O1
      O    -1.617979    -0.948062    -2.341650    H2O2
      H    -2.529195    -1.296822    -2.122437    H2O2
      H    -1.534288    -0.833088    -3.331486    H2O2
      O    -1.507990     2.117783     1.555094    H2O3
      H    -1.501128     2.645178     2.403050    H2O3
      H    -2.090603     1.352766     1.597519    H2O3
    &END COORD
    &TOPOLOGY
    &END TOPOLOGY
  &END SUBSYS
&END FORCE_EVAL
