On the one hand you may object that it is an additional burden to write Matlab functions for the individual components of the force field. On the other hand, it gives you additional flexibility and you can perform simulations not directly possible with other standard molecular simulation packages.
The potential energy (aka force field) of a general system of molecules of different species consists of three types of contributions
Intramolecular interaction, i.e., interactions within each molecule
Intermolecular interaction, i.e., interactions between molecules of the same species)
Intermolecular interaction, i.e., interactions between molecules of different species)
Efficient implementation of gradient and Hessian crucially relies on this function
please read this guide
please read this guide
Wiki: Developers.ForceFields.inter
Wiki: Developers.ForceFields.intra
Wiki: Developers.Main
Wiki: Users.Main