INTERACTION subblock in QUILD block

One of the most important input-parts for multi-level jobs is the INTERACTION subblock of QUILD, where one should define how the different descriptions should be applied to the different regions. At the part where we explained the multi-level approaches, we already showed some examples of how to combine different methodologies. Below is another example input where all possible options are given.

QUILD
  INTERACTION
    TOTAL     description 1
    REPLACE   region 1 region 2    description 3 for description 2
    REPLACE   region 1             description 4 for description 3
    INTXN     region 1 region 2    description 3 for description 2
  SUBEND
END

There should always be a line with the description of the total system, as shown in the first line of the INTERACTION subblock. Then if you want to replace the interactions for one (or more) region(s), you could do so as indicated in the second and third line. Finally, if you want to replace the interaction between two regions, as we need for DNA where we replace the BP86 π-stacking by LDA π-stacking, the last line of the INTERACTION subblock should be used.

Note that in all cases it is not necessary at all to add the "region", "description" and "for" words in the INTERACTION subblock; they are ignored when reading the input. The program reads the line, uses the last two integers for the descriptions and the ones before for the regions. Therefore, a completely equivalent input would be:

QUILD
  INTERACTION
    TOTAL             1
    REPLACE   1 2   3 2
    REPLACE   1     4 3
    INTXN     1 2   3 2
  SUBEND
END

However, for better readability, it is to be advised to always use the additional text anyway.