<div dir="auto"><div dir="auto">Hi everyone, <br></div><div dir="auto">I would like to use a force field in GROMACS coming from my Trained Neural Network. My Neural Network function represents the local potential energy of an atom i , E_{i}. </div><div dir="auto"><br></div><div dir="auto">E_{i} depends on, say, n descriptors (D_{i1},...,D_{in}), E_{i}(D_{i1},...,D_{in}) . These descriptors depend on coordinates of atom i and coordinates of the first, say, 20 nearest neighbour  atoms. </div><div dir="auto"><br></div><div dir="auto">I found the Tabulated interaction functions topics in the reference manual. It seems possible to use a custom pairwise interaction E(r_{ij})  but not more. </div><div dir="auto"><br></div><div dir="auto">I&#39;m looking for a routine where:</div><div dir="auto"><br></div><div dir="auto">I can compute at each step the descriptors for each molecule then evaluate all E_{i} and then derivate it to get forces. </div><div dir="auto"><br></div><div dir="auto">or </div><div dir="auto"><br></div><div dir="auto">I can tabulate all the descriptors and the corresponding energies and forces.</div><div dir="auto"><br></div><div dir="auto">Is there a way to do something like that in GROMACS? If I have to modify some source codes, how it could be difficult? </div><div dir="auto"><br></div><div dir="auto">FGM</div></div>