<div>Dear all,</div><div><br /></div><div>I have encountered a very disturbing phenomena. </div><div><br /></div><div>My system is fairly simple: One Na+ in water (up to 800 water molecules). I generated the system using tleap. The plan is the following:</div><div><br /></div><div>Since full semi empirical MD with XTB is costly for 801 atoms, I  tried to calculate MD with both conventional QMMM and abfQMMM. I do acknowledge, that conventional QMMM will not work since the QM atoms in the QM box change. Therefore abfQMMM seems to be the solution. In any case, I tried conventional QMMM first. <br /></div><div><br /></div><div>1) conventional QMMM<br /></div><div>Following the normal procedure, I did first minimization without QMMM, and it crashed for most of the time. Somehow it depends on the box size and the binary I execute. With any executable ssmp, psmp, popt, sopt using CPU only, it crashed. The only successful run was with the local_cuda executable sopt. Sometimes it also did not work. </div><div><br /></div><div>I run the cp2k.pdbg and cp2k.sdbg. The error message when I run cp2k.sdbg was:</div><div><br /></div><div><b>At line 766 of file /usr/local/cp2k/src/force_fields_ext.F<br />Fortran runtime error: Pointer argument 'amb_info' is not associated<br /></b><br />Error termination. Backtrace:<br />#0  0x2183c84 in __force_fields_ext_MOD_read_force_field_amber<br /><span style="white-space: pre;">        </span>at /usr/local/cp2k/src/force_fields_ext.F:766<br />#1  0x1d8aa45 in __force_fields_MOD_force_field_control<br /><span style="white-space: pre;">      </span>at /usr/local/cp2k/src/force_fields.F:161<br />#2  0x19d5df9 in __fist_environment_MOD_fist_init<br /><span style="white-space: pre;">        </span>at /usr/local/cp2k/src/fist_environment.F:176<br />#3  0x184bc1b in __fist_main_MOD_fist_create_force_env<br /><span style="white-space: pre;">       </span>at /usr/local/cp2k/src/fist_main.F:98<br />#4  0x15ed06b in __f77_interface_MOD_create_force_env<br /><span style="white-space: pre;">        </span>at /usr/local/cp2k/src/f77_interface.F:786<br />#5  0x473b25 in cp2k_run<br /><span style="white-space: pre;">        </span>at /usr/local/cp2k/src/start/cp2k_runs.F:308<br />#6  0x475a19 in __cp2k_runs_MOD_run_input<br /><span style="white-space: pre;">     </span>at /usr/local/cp2k/src/start/cp2k_runs.F:975<br />#7  0x46b112 in cp2k<br /><span style="white-space: pre;">  </span>at /usr/local/cp2k/src/start/cp2k.F:379<br />#8  0x46b25d in main<br /><span style="white-space: pre;">       </span>at /usr/local/cp2k/src/start/cp2k.F:44</div><div><br /></div><div>The sporadic errornous and successful runs continue as well in the following NVT and NPT equilibration with QMMM. For me the error message seems to originate in the amber force field parser.<br /></div><div><br /></div><div>2) abfQMMM. No matter what executable I called, everytime it crashed and the same error message appears. I enlarged the box, still no success.</div><div><br /></div><div>Best regards,</div><div><br /></div><div>T. Oka<br /></div><div><br /></div><div><br /></div>

<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+unsubscribe@googlegroups.com">cp2k+unsubscribe@googlegroups.com</a>.<br />
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/a4862465-69e9-4482-a737-cefb11b9c2b4n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/a4862465-69e9-4482-a737-cefb11b9c2b4n%40googlegroups.com</a>.<br />