<div dir="ltr">Hi Christopher<br><br>if I am not mistaken, all nonbonded pair potentials are tabulated using splines in the MM pre-processing phase and these are then used as lookup tables during the simulation. I think this is a quite common approach which keeps the computational effort at the same level independent of the algebraic complexity of the actual pair potential. A control of the mapping (splining) accuracy is possible by CP2K input parameters (e.g. EPS_SPLINE). That might answer your first two questions at least partially. As far as I know, it is not possible to avoid that tabulation (splining) of the nonbonded pair potentials currently without a code change. The Poisson solver follows a different approaches depending on the chosen method.<br><br>It would be helpful if you would describe your problems with the splines more precisely. The splining of the pre-defined pair potentials works usually fine. You can activate SPLINE_DATA and SPLINE_INFO in the section <big class="uctt"><font size="2"><a href="https://manual.cp2k.org/cp2k-3_0-branch/CP2K_INPUT.html">CP2K_INPUT</a> /
               <a href="https://manual.cp2k.org/cp2k-3_0-branch/CP2K_INPUT/FORCE_EVAL.html">FORCE_EVAL</a> /
               <a href="https://manual.cp2k.org/cp2k-3_0-branch/CP2K_INPUT/FORCE_EVAL/MM.html">MM</a> /
               <a href="https://manual.cp2k.org/cp2k-3_0-branch/CP2K_INPUT/FORCE_EVAL/MM/PRINT.html">PRINT</a> /
               <a href="https://manual.cp2k.org/cp2k-3_0-branch/CP2K_INPUT/FORCE_EVAL/MM/PRINT/FF_INFO.html">FF_INFO</a></font><br></big>This will print a data file for each nonbonded pair potential interaction (atomic kind pair) using splines.<br><br>HTH<br><br>Matthias<br><br><br>On Wednesday, 4 May 2016 23:07:34 UTC+2, Christopher Mauney  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>I'm running into a lot of problems with splines.  The documentation on these is frustratingly imprecise and vague.<div><br></div><div>Why does CP2K use splines by default in MM calculations?</div><div>Are we modelling potentials as a spline? (That is, are we fitting a spline based on finding the values of the potential on some underlying grid?)</div><div>Can we turn off this behavior and just calculate the potential values directly?<br></div><div>What kind of spline is being used? Can we modify this?</div><div>Does the Poisson solver use splines the same way?</div><div><br></div><div>Any help here is appreciated.  I do have specific problems with input files, but this is an aspect of CP2K I've had trouble with on several projects and would prefer a general response.</div><div><br></div><div>-Christopher Mauney</div><div><br></div></div></blockquote></div>