&FORCE_EVAL
  METHOD QMMM
  &DFT
    BASIS_SET_FILE_NAME ../BASIS_SET
    POTENTIAL_FILE_NAME ../POTENTIAL
    &MGRID
      COMMENSURATE
      CUTOFF 50
    &END MGRID
    &QS
    &END QS
    &SCF
      SCF_GUESS atomic
    &END SCF
    &XC
      &XC_FUNCTIONAL pade
      &END XC_FUNCTIONAL
    &END XC
  &END DFT
  &MM
    &FORCEFIELD
      parm_file_name ../sample_pot/2H2O.pot
      parmtype CHM
      &spline
       rcut_nb 5.2917720830
      &end
    &END FORCEFIELD
    &POISSON
      &EWALD
        EWALD_TYPE spme
        ALPHA .44
        GMAX 64
        O_SPLINE 6
      &END EWALD
    &END POISSON
  &END MM
  &QMMM
    MM_POTENTIAL_FILE_NAME ../MM_POTENTIAL
    &CELL
      ABC 7.0 7.0 7.0
    &END CELL
    ECOUPL NONE
    NOCOMPATIBILITY
    &QM_KIND H
      MM_INDEX 5 6
    &END QM_KIND
    &QM_KIND O
      MM_INDEX 4
    &END QM_KIND
  # 
  # QM_KINDS
  # 
  &END QMMM
  &SUBSYS
    &CELL
      ABC 20.00 20.00 20.00
    &END CELL
    &KIND H
      BASIS_SET DZVP-GTH-BLYP
      POTENTIAL GTH-BLYP-q1
    &END KIND
    &KIND O
      BASIS_SET DZVP-GTH-BLYP
      POTENTIAL GTH-BLYP-q6
    &END KIND
    &TOPOLOGY
      CONN_FILE_NAME ../sample_psf/2H2O.psf
      CONNECTIVITY UPSF
      COORD_FILE_NAME ../sample_pdb/2H2O.pdb
      COORDINATE pdb
    &END TOPOLOGY
  &END SUBSYS
&END FORCE_EVAL
&GLOBAL
  PROJECT H2O-qmmm-none-1
  PRINT_LEVEL MEDIUM
  RUN_TYPE GEO_OPT
&END GLOBAL
&MOTION
  &GEO_OPT
    MAX_ITER 3
    OPTIMIZER BFGS
  &END GEO_OPT
&END MOTION
