<cg>
  <bonded>
    <!-- name of the interaction -->
    <name>bond</name>
    <min>2.0</min>
    <max>3.0</max>
    <step>0.01</step>
    <inverse>
      <!-- 300*0.0019872036 lammps "real" units kcal/mol -->
      <kBT>0.59616108</kBT>
      <!-- target distribution -->
      <target>bond.dist.tgt</target>
      <lammps>
        <kBT>0.59616108</kBT>
        <table>table_b1.pot</table>
      </lammps>
    </inverse>
  </bonded>
  <bonded>
    <!-- name of the interaction -->
    <name>angle</name>
    <min>1.5</min>
    <max>3.1415</max>
    <step>0.01</step>
    <inverse>
      <!-- target distribution -->
      <target>angle.dist.tgt</target>
      <lammps>
      <kBT>0.59616108</kBT>
        <table>table_a1.pot</table>
      </lammps>
    </inverse>
  </bonded>
  <!-- example for a non-bonded interaction entry -->
  <non-bonded>
    <!-- name of the interaction -->
    <name>A-A</name>
    <!-- types involved in this interaction -->
    <type1>A</type1>
    <type2>A</type2>
    <!-- dimension + grid spacing of tables for calculations -->
    <min>2.0</min>
    <max>12.0</max>
    <step>0.05</step>
    <inverse>
      <kBT>0.59616108</kBT>
      <p_target>1.0</p_target>
      <!-- target distribution (rdf), just give gromas rdf.xvg -->
      <target>A-A.dist.tgt</target>
      <!-- update cycles -->
      <do_potential>1 0 0</do_potential>
      <!-- additional post processing of dU before added to potential -->
      <post_update>pressure</post_update>
      <post_update_options>
        <pressure>
          <type>simple</type>
          <do>1 0 0</do>
          <simple>
            <scale>0.003</scale>
          </simple>
        </pressure>
      </post_update_options>
      <!-- additional post processing of U after dU added to potential -->
      <post_add>convergence</post_add>
      <!-- name of the table for gromacs run -->
      <lammps>
        <table>table_A_A.pot</table>
      </lammps>
    </inverse>
  </non-bonded>
  
  <non-bonded>
    <!-- name of the interaction -->
    <name>B-B</name>
    <!-- types involved in this interaction -->
    <type1>B</type1>
    <type2>B</type2>
    <!-- dimension + grid spacing of tables for calculations -->
    <min>2.0</min>
    <max>12.0</max>
    <step>0.05</step>
    <inverse>
      <kBT>0.59616108</kBT>
      <p_target>1.0</p_target>
      <!-- target distribution (rdf), just give gromas rdf.xvg -->
      <target>B-B.dist.tgt</target>
      <!-- update cycles -->
      <do_potential>0 1 0</do_potential>
     <post_update>pressure</post_update>
      <post_update_options>
        <pressure>
          <type>simple</type>
          <do>0 1 0</do>
          <simple>
            <scale>0.003</scale>
          </simple>
        </pressure>
      </post_update_options>
      <!-- additional post processing of dU before added to potential -->
      <post_update></post_update>
      <!-- additional post processing of U after dU added to potential -->
      <post_add>convergence</post_add>
      <!-- inverse monte carlo specific stuff -->
      <!-- name of the table for gromacs run -->
      <lammps>
        <kBT>0.59616108</kBT>
        <table>table_B_B.pot</table>
      </lammps>
    </inverse>
  </non-bonded>

  <non-bonded>
    <!-- name of the interaction -->
    <name>A-B</name>
    <!-- types involved in this interaction -->
    <type1>A</type1>
    <type2>B</type2>
    <!-- dimension + grid spacing of tables for calculations -->
    <min>2.0</min>
    <max>12.0</max>
    <step>0.05</step>
    <inverse>
      <kBT>0.59616108</kBT>
      <p_target>1.0</p_target>
      <!-- target distribution (rdf), just give gromas rdf.xvg -->
      <target>A-B.dist.tgt</target>
      <!-- update cycles -->
      <do_potential>0 0 1</do_potential>
      <!-- additional post processing of dU before added to potential -->
      <post_update>pressure</post_update>
      <post_update_options>
        <pressure>
          <type>simple</type>
          <do>0 0 1</do>
          <simple>
            <scale>0.003</scale>
          </simple>
        </pressure>
      </post_update_options>
      <!-- additional post processing of U after dU added to potential -->
      <post_add>convergence</post_add>
      <!-- inverse monte carlo specific stuff -->
      <!-- name of the table for gromacs run -->
      <lammps>
        <kBT>0.59616108</kBT>
        <table>table_A_B.pot</table>
      </lammps>
    </inverse>
  </non-bonded>

  <!-- general options for inverse script -->
  <inverse>
    <kBT>0.59616108</kBT>
    <!-- use lammps as simulation program -->
    <program>lammps</program>
    <!-- lammps specific options -->
    <lammps>
      <command>mpirun -np 32 lmp</command>
      <script>hexane_cg.in</script>
      <topol>topol.xml</topol>
      <traj>traj.dump</traj>
    </lammps>
    <map>hexane_cg.xml</map>
    <kBT>0.59616108</kBT>
    <initial_configuration>maindir</initial_configuration>
    <!-- these files are copied for each new run -->
    <filelist>hexane_cg.data hexane_cg.in topol.xml hexane_cg.dump</filelist>
    <!-- do so many iterations -->
    <iterations_max>500</iterations_max>
    <convergence_check>
      <type>default</type>
      <limit>0.05</limit>
    </convergence_check>
    <!-- ibm: inverse biltzmann imc: inverse monte carlo -->
    <method>ibi</method>
    <!-- write log to this file -->
    <log_file>inverse.log</log_file>
    <!-- write restart step to this file -->
    <restart_file>restart_points.log</restart_file>
  </inverse>
</cg>
