<div dir="ltr"><div><div><div><div>Hi,<br><br></div>You can get it to run by adding a section in &FORCEFIELD:<br></div><div style="margin-left:40px"><span style="font-family:courier new,monospace">&SPLINE</span><br><span style="font-family:courier new,monospace"></span><span style="font-family:courier new,monospace">    EMAX_SPLINE 3.0</span><br><span style="font-family:courier new,monospace"></span><span style="font-family:courier new,monospace">&END SPLINE</span><br><br></div>However, starting from the first step, the system basically explodes and all you get is <span style="font-family:courier new,monospace">NaN</span>s.<br></div><div><br></div><div>Could it be that something is wrong in your TIP4P forcefield definition? If I strip down the system to a single H2O molecule, the simulation runs, but the 'He' dummy atom immediately flies away from the molecule to the other side of the box.<br><br></div><div>Ralph<br></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Nov 13, 2014 at 4:02 PM, Chris Lee <span dir="ltr"><<a href="mailto:ca.l...@gmail.com" target="_blank">ca.l...@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Dear CP2K users,<br><br>I'm trying to run classical simulations of portlandite and TIP4P water however I run into the following error:<br><br> WARNING| Particles:    5341   1919 at distance [au]:     5.11334331 less than:      5.93645406; increase EMAX_SPLINE.<br><br><br> ******************************************************************<br> *** ERROR in build_neighbor_lists (MODULE fist_neighbor_lists) ***<br> ******************************************************************<br><br> *** GEOMETRY wrong or EMAX_SPLINE too small! ***<br><br> *** Program stopped at line number 623 of MODULE fist_neighbor_lists ***<br><br> ===== Routine Calling Stack =====<br><br>            6 build_neighbor_lists<br>            5 build_fist_neighbor_lists<br>            4 list_control<br>            3 fist_force_control<br>            2 qs_mol_dyn_low<br>            1 CP2K<br> CP2K| Abnormal program termination, stopped by process number 0<br><br><br>From what I can tell the problem seems to be with the water molecules as separate calculations on only the portlandite slab in the box work fine, but the calculations on only water fail with a similar error. The pdb file was created from the output from another molecular dynamics program (DLPOLY). Is the error to do with some of the water molecules crossing the periodic boundary and is there any way around this? Or will I have to solvate the slab making sure that none of the water molecules cross the boundary?<br><br>I've attached the input files that I've used and the output file produced.<br><br>Any help would be greatly appreciated, thank you.<span class="HOEnZb"><font color="#888888"><br></font></span></div><span class="HOEnZb"><font color="#888888">

<p></p>

-- <br>
You received this message because you are subscribed to the Google Groups "cp2k" group.<br>
To unsubscribe from this group and stop receiving emails from it, send an email to <a href="mailto:cp2k+uns...@googlegroups.com" target="_blank">cp2k+uns...@googlegroups.com</a>.<br>
To post to this group, send email to <a href="mailto:cp...@googlegroups.com" target="_blank">cp...@googlegroups.com</a>.<br>
Visit this group at <a href="http://groups.google.com/group/cp2k" target="_blank">http://groups.google.com/group/cp2k</a>.<br>
For more options, visit <a href="https://groups.google.com/d/optout" target="_blank">https://groups.google.com/d/optout</a>.<br>
</font></span></blockquote></div><br></div>