The central idea observes a recursive mapping of -body intramolecular interactions to -body terms that is consistent with the molecular topology. Iterative application of the line graph transformation is identified as a natural and elegant tool to accomplish the recursion. The procedure readily generalizes to arbitrary -body potentials. In particular, the method yields a complete characterization of -body interactions. The hierarchical structure of atomic index lists for each interaction order is compactly expressed as a directed acyclic graph. A pseudo-code description of the generating algorithm is given. With suitable data structures (e.g., edge lists or adjacency matrices), automatic enumeration and indexing of -body interactions can be implemented straightforwardly to handle large bio-molecular systems. Explicit examples are discussed, including a chemically relevant effective potential model of taurocholate bile salt.
|Number of pages||15|
|Journal||Journal of Mathematical Chemistry|
|Early online date||19 May 2015|
|Publication status||Published - Aug 2015|
- Computer simulation
- Molecular modelling
- Graph theory