&GLOBAL
  PRINT_LEVEL HIGH
  PROJECT sod
  RUN_TYPE MD
&END GLOBAL

&MOTION
  &MD
    ENSEMBLE NVE
    STEPS 10
    TEMPERATURE 300
    TEMP_TOL 50
    TIMESTEP 1.0
  &END MD
&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 1
      &END CHARGE
      &CHARGE
        ATOM O
        CHARGE -0.5
      &END CHARGE
      &NONBONDED
        &GENPOT
          ATOMS Si Si
          FUNCTION 1*exp(-r/2.0)
          VARIABLES r
        &END GENPOT
        &GENPOT
          ATOMS Si O
          FUNCTION 1*exp(-r/2.0)
          VARIABLES r
        &END GENPOT
        &GENPOT
          ATOMS O O
          FUNCTION 1*exp(-r/2.0)
          VARIABLES r
        &END GENPOT
      &END NONBONDED
      &NONBONDED14
        &GENPOT
          ATOMS Si O
          FUNCTION 1*exp(-r/2.0)
          VARIABLES r
        &END GENPOT
      &END NONBONDED14
      &SPLINE
        EMAX_ACCURACY 1.0
        EMAX_SPLINE 1.0
        R0_NB [bohr] 3.5
      &END SPLINE
    &END FORCEFIELD
    &POISSON
      &EWALD
        EWALD_TYPE NONE
      &END EWALD
    &END POISSON
    &PRINT
      #&NEIGHBOR_LISTS
      #&END NEIGHBOR_LISTS
      &FF_INFO DEBUG
      &END FF_INFO
      &SUBCELL
      &END SUBCELL
    &END PRINT
  &END MM
  &PRINT
    &FORCES
    &END FORCES
  &END PRINT
  &SUBSYS
    &CELL
      ABC 8.965 8.965 8.965
    &END CELL
    &PRINT
      &CELL
        UNIT nm
      &END CELL
      #&INTERATOMIC_DISTANCES
      #&END INTERATOMIC_DISTANCES
      #&TOPOLOGY_INFO
      #    XYZ_INFO T
      #    GENERATE_INFO F
      #&END TOPOLOGY_INFO
      #&KINDS
      #&END KINDS
    &END PRINT
    &TOPOLOGY
      CONNECTIVITY UPSF
      CONN_FILE_NAME ../sample_psf/sod.psf
      COORDINATE XYZ
      COORD_FILE_NAME ../sample_xyz/sod.xyz
      #&GENERATE
      #    CREATE_MOLECULES T
      #&END GENERATE
    &END TOPOLOGY
  &END SUBSYS
&END FORCE_EVAL
