&GLOBAL
  PRINT_LEVEL LOW
  PROJECT mfi_geo3
  RUN_TYPE GEO_OPT
  SAVE_MEM
&END GLOBAL

&MOTION
  &GEO_OPT
    MAX_ITER 40
    MINIMIZER CG
    &CG
      &LINE_SEARCH
        TYPE 2PNT
        &2PNT
        &END 2PNT
      &END LINE_SEARCH
    &END CG
  &END GEO_OPT
&END MOTION

&FORCE_EVAL
  METHOD FIST
  &MM
    &FORCEFIELD
      &BEND
        ATOMS Si O Si
        K 0.0
        KIND HARMONIC
        THETA0 0.0
      &END BEND
      &BEND
        ATOMS O Si O
        K 0.0
        KIND HARMONIC
        THETA0 0.0
      &END BEND
      &BOND
        ATOMS Si O
        K 0.5
        KIND HARMONIC
        R0 [angstrom] 1.59
      &END BOND
      &CHARGE
        ATOM Si
        CHARGE 0.0
      &END CHARGE
      &CHARGE
        ATOM O
        CHARGE 0.0
      &END CHARGE
      &NONBONDED
        &GENPOT
          ATOMS Si Si
          FUNCTION 0.0
          VARIABLES r
        &END GENPOT
        &GENPOT
          ATOMS Si O
          FUNCTION 0.0
          VARIABLES r
        &END GENPOT
        &GENPOT
          ATOMS O O
          FUNCTION 0.0
          VARIABLES r
        &END GENPOT
      &END NONBONDED
      &SPLINE
        R0_NB 0.1
      &END SPLINE
    &END FORCEFIELD
    &POISSON
      &EWALD
        GMAX 25
      &END EWALD
    &END POISSON
    &PRINT
      #&NEIGHBOR_LISTS
      #&END NEIGHBOR_LISTS
      &FF_INFO DEBUG
      &END FF_INFO
      &SUBCELL
      &END SUBCELL
    &END PRINT
  &END MM
  &SUBSYS
    &CELL
      ABC 20.413734 19.904612 13.453377
    &END CELL
    &TOPOLOGY
      COORDINATE XYZ
      COORD_FILE_NAME ../sample_xyz/mfi.xyz
      &GENERATE
        CREATE_MOLECULES T
      &END GENERATE
    &END TOPOLOGY
  &END SUBSYS
&END FORCE_EVAL
