<html xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40"> <head> <meta http-equiv="Content-Type" content="text/html; charset=Windows-1252"> <meta name="Generator" content="Microsoft Word 15 (filtered medium)"> <style></style> </head> <body lang="en-CH" link="blue" vlink="purple" style="word-wrap:break-word"> <div class="WordSection1"> <p class="MsoNormal"><span lang="DE-CH" style="font-size:11.0pt;mso-fareast-language:EN-US">Hi Léon<o:p></o:p></span></p> <p class="MsoNormal"><span lang="DE-CH" style="font-size:11.0pt;mso-fareast-language:EN-US"><o:p> </o:p></span></p> <p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;mso-fareast-language:EN-US">your force field definition is incorrect. The unit for the C parameters should be [kcalmol*angstrom^-6] instead of [kcalmol*angstrom^-1] and the B parameter for Ti-O should be something like 5.15 rather than 0.00515. The simulation cell is also too small for the given RCUT value of 12, but that shouldn’t matter.<o:p></o:p></span></p> <p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;mso-fareast-language:EN-US"><o:p> </o:p></span></p> <p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;mso-fareast-language:EN-US">HTH<o:p></o:p></span></p> <p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;mso-fareast-language:EN-US"><o:p> </o:p></span></p> <p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;mso-fareast-language:EN-US">Matthias<o:p></o:p></span></p> <p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;mso-fareast-language:EN-US"><o:p> </o:p></span></p> <div id="mail-editor-reference-message-container"> <div> <div style="border:none;border-top:solid #B5C4DF 1.0pt;padding:3.0pt 0cm 0cm 0cm"> <p class="MsoNormal" style="mso-margin-top-alt:0cm;margin-right:0cm;margin-bottom:12.0pt;margin-left:36.0pt"> <b><span style="color:black">From: </span></b><span style="color:black">cp2k@googlegroups.com <cp2k@googlegroups.com> on behalf of Léon Luntadila Lufungula <Leon.luntadilalufungula@uantwerpen.be><br> <b>Date: </b>Wednesday, 13 March 2024 at 12:20<br> <b>To: </b>cp2k <cp2k@googlegroups.com><br> <b>Subject: </b>[CP2K:20022] Matsui-Akaogi Force Field for bulk TiO2<o:p></o:p></span></p> </div> <p class="MsoNormal" style="margin-left:36.0pt">Dear all,<o:p></o:p></p> <div> <p class="MsoNormal" style="margin-left:36.0pt"><o:p> </o:p></p> </div> <div> <p class="MsoNormal" style="margin-left:36.0pt">I want to test a force field for TiO2 so that I can use it for pre-equilibrating my system for AIMD calculations. However, I have no experience in using force fields and I'm having trouble implementing the FF calculations. I took the FF parameters from a paper from <a href="https://pubs.acs.org/doi/10.1021/jp050963q"> Naicker et al</a> and tried to calculate the energy of the anatase crystal structure, before going on to a cell optimization of the structure to compare to DFT. <o:p></o:p></p> </div> <div> <p class="MsoNormal" style="margin-left:36.0pt"><o:p> </o:p></p> </div> <div> <p class="MsoNormal" style="margin-left:36.0pt">However, when I try to run the calculation I get an error indicating "GEOMETRY wrong or EMAX_SPLINE too small!", however, the geometry is definitely correct. So I tried increasing EMAX_SPLINE to 1E2 and then got "SPLINE_INFO| Number of points: 2346797 obtained accuracy 0.953210E-07. MM SPLINE: no convergence on required accuracy (adjust EPS_SPLINE and rerun)". I then tried to increase EPS_SPLINE to 1E-6 and the energy calculation ran without a problem. I tried to run a CELL_OPT with these settings but the optimization seems to be going nowhere... (see attached input and output files) Can anyone tell me what I'm doing wrong? Is it perhaps the definition of the forcefield?<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="margin-left:36.0pt"><o:p> </o:p></p> </div> <div> <p class="MsoNormal" style="margin-left:36.0pt">Kind regards,<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="margin-left:36.0pt">Léon<o:p></o:p></p> </div> <p class="MsoNormal" style="margin-left:36.0pt">-- <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/e185a18f-fcbc-477a-ba8e-95c0191b4ae4n%40googlegroups.com?utm_medium=email&utm_source=footer"> https://groups.google.com/d/msgid/cp2k/e185a18f-fcbc-477a-ba8e-95c0191b4ae4n%40googlegroups.com</a>.<o:p></o:p></p> </div> </div> </div> </body> </html> <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/ZRAP278MB0827BFAFF5DE0330FCF5EB48F42A2%40ZRAP278MB0827.CHEP278.PROD.OUTLOOK.COM?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/ZRAP278MB0827BFAFF5DE0330FCF5EB48F42A2%40ZRAP278MB0827.CHEP278.PROD.OUTLOOK.COM</a>.<br />