<div dir="ltr">Hi,<div><br></div><div>I am not 100% sure but it is possible that the Tersoff potential is not compatible with one non bonding potential.</div><div>You can solve the problem running a multiple force environment calculation, where you define two independent force environments, one for</div><div>Tersoff and one for LJ. The global force should be the simple sum of the two contributions. </div><div>Please check for reference the examples of multiple force environments that you can find among the cp2k/tests</div><div>Kind regards</div><div>Marcella</div><div><br><br>On Wednesday, February 13, 2019 at 11:05:43 AM UTC+1, p clabaut wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><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></blockquote></div></div>