<div dir="ltr"><div>Hello,</div><div><br></div><div>I'm trying to use a manybody potential between a metallic slab and water molecules in FIST but I would also need to use a simple Lennard-Jones potential to govern my water/water interactions.</div><div><br></div><div>Unfortunately, it would seems that the neighbour list building went badly between those modules because as long as my manybody potential (here, Tersoff) is present in the input, some of the Oxygen/Oxygen Lennard-Jones interactions are ignored.</div><div>Moreover, the verlet_skin parameters starts switching on or off some of these interactions.</div><div><br></div><div>For an example, I looked at two couples of water molecules away from my Pt slab and distant for each other.</div><div><br></div><div>- When I use <br></div><div>    &NEIGHBOR_LISTS<br>      VERLET_SKIN 0.00<br>    &END</div><div>with </div><div>    &TERSOFF<br>        ATOMS O Pt<br>        RCUT [angstrom] 10<br>     &END TERSOFF</div><div>and <br></div><div>    &LENNARD-JONES<br>          atoms O O<br>          EPSILON [kcalmol] 0.102<br>          SIGMA [angstrom] 3.188<br>          RCUT 10<br>     &END LENNARD-JONES</div><div><br></div><div>the oxygens of only one of my water couple are experiencing LJ interaction from each others while the others doesn'tfeel any forces at all <br></div><div><br></div><div><div>- When I use <br></div><div>    &NEIGHBOR_LISTS<br>      VERLET_SKIN 1.00<br>    &END</div><div>with </div><div>    &TERSOFF<br>        ATOMS O Pt<br>        RCUT [angstrom] 10<br>     &END TERSOFF</div><div>and <br></div><div>    &LENNARD-JONES<br>          atoms O O<br>          EPSILON [kcalmol] 0.102<br>          SIGMA [angstrom] 3.188<br>          RCUT 10<br>     &END LENNARD-JONES</div><div><br></div><div>with no others changes, it is the second water couple that undergoes LJ interaction but not the first one anymore.</div><div><br></div><div><div>- When I use <br></div><div>    &NEIGHBOR_LISTS<br>      VERLET_SKIN 1.00<br>    &END</div><div>and <br></div><div>    &LENNARD-JONES<br>          atoms O O<br>          EPSILON [kcalmol] 0.102<br>          SIGMA [angstrom] 3.188<br>          RCUT 10<br>     &END LENNARD-JONES</div><div><br></div><div>with no manybodypotential, both couples are resenting LJ interactions.</div><div><br></div><div><br></div><div><br></div><div>I have observed the same exact behavior with a Williams potential instead of a Lennard-Jones. I'm joining the example input in the message.</div><div><br></div><div>Could you tell me what is going and if it is intended or a bug in the programm?</div><div><br></div><div>Thanks you by advance !<br></div></div></div></div>