<div dir="ltr"><div>I am trying to use FIST to do a box of water molecules with a force field (de Leeuw N. H., and Parker S. C. (1998a) Molecular-dynamics simulations of MgO surfaces in liquid water using a shell-model potential for water. Phys. Rev. B 58, 13901–13908.). </div><div><br></div><div>There is the coulombic subtraction in Morse potential for bond.</div><div><br></div><div>In CP2K the morse potential is in the form</div><div><br></div><div><em style="color: rgb(0, 0, 0); font-family: "Times New Roman"; font-size: medium;">K1*[(1-exp(-K2*(R-R0)))^2-1])</em><span style="color: rgb(0, 0, 0); font-family: "Times New Roman"; font-size: small;"> </span><br></div><div><span style="color: rgb(0, 0, 0); font-family: "Times New Roman"; font-size: small;"><br></span></div><div><span style="color: rgb(0, 0, 0); font-family: "Times New Roman"; font-size: small;">But I want to know how to add the the term of coulombic interaction.</span></div><div><br></div><div>I am not sure if it is possible in CP2K. </div><div><br></div><div>Please give me some suggestions. Thank you very much.</div></div>