Dear Matthias,<div><br /></div><div>Thanks for pointing me to the right NONBONDED implementation, it seems to work well with
<a href="https://manual.cp2k.org/trunk/CP2K_INPUT/FORCE_EVAL/MM/FORCEFIELD/NONBONDED/BUCKMORSE.html" target="_blank" rel="nofollow" style="color: rgb(26, 115, 232); font-size: 14.6667px;">BUCKMORSE</a>. I also found my error with the other WILLIAMS parameters, I messed up the exponent for the units of C (kcalmol*angstrom^-6 instead of kcalmol*angstrom^6), so now they are working as they should! I was wondering though, how important it is to make sure your cell dimensions are larger than RCUT? I tried redoing the calculation with a 1x1x1 cell instead of 8x8x3 and the results seem to be more or less the same...</div><div><br /></div><div>Kind regards,</div><div>Léon</div><div><br /></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Thursday 14 March 2024 at 20:52:55 UTC+1 Krack Matthias wrote:<br/></div><blockquote class="gmail_quote" style="margin: 0 0 0 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">





<div lang="en-CH" link="blue" vlink="purple" style="word-wrap:break-word">
<div>
<p class="MsoNormal"><span lang="DE-CH" style="font-size:11.0pt">Hi Léon<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="DE-CH" style="font-size:11.0pt"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt">I think there is no need to use GENPOT, because that nonbonded force field contribution is implemented in CP2K as
<a href="https://manual.cp2k.org/trunk/CP2K_INPUT/FORCE_EVAL/MM/FORCEFIELD/NONBONDED/BUCKMORSE.html" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=en-GB&q=https://manual.cp2k.org/trunk/CP2K_INPUT/FORCE_EVAL/MM/FORCEFIELD/NONBONDED/BUCKMORSE.html&source=gmail&ust=1710583154052000&usg=AOvVaw03dNWRm5JutwLAp6h7IcFG">
BUCKMORSE</a>. There are often typos in papers with force field parameters and it is advisable to cross-check the parameter tables for the same force field in different papers for consistency.
<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt">Small differences to other codes can have many reasons: splining quality, single- vs double-precision, SPME mesh size, type and threshold values of the employed optimizer
 etc.<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt">Best<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt">Matthias<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt"><u></u> <u></u></span></p>
<div>
<div></div></div></div></div><div lang="en-CH" link="blue" vlink="purple" style="word-wrap:break-word"><div><div><div>
<div style="border:none;border-top:solid #b5c4df 1.0pt;padding:3.0pt 0cm 0cm 0cm">
<p class="MsoNormal" style="margin-right:0cm;margin-bottom:12.0pt;margin-left:36.0pt">
<b><span style="color:black">From: </span></b><span style="color:black"><a href data-email-masked rel="nofollow">cp...@googlegroups.com</a> <<a href data-email-masked rel="nofollow">cp...@googlegroups.com</a>> on behalf of Léon Luntadila Lufungula <<a href data-email-masked rel="nofollow">Leon.luntad...@uantwerpen.be</a>><br>
<b>Date: </b>Thursday, 14 March 2024 at 11:48<br>
<b>To: </b>cp2k <<a href data-email-masked rel="nofollow">cp...@googlegroups.com</a>><br>
<b>Subject: </b>Re: [CP2K:20027] Matsui-Akaogi Force Field for bulk TiO2<u></u><u></u></span></p>
</div>
<p class="MsoNormal" style="margin-left:36.0pt">Hi Matthias,<u></u><u></u></p>
<div>
<p class="MsoNormal" style="margin-left:36.0pt"><u></u> <u></u></p>
</div>
</div></div></div></div><div lang="en-CH" link="blue" vlink="purple" style="word-wrap:break-word"><div><div><div><div><p class="MsoNormal" style="margin-left:36.0pt">Thanks for pointing this out for me, I've adjusted the input accordingly and now the calculation seems to work fine. I've increased GMAX to 2 points per angstrom and also used MULTIPLE_UNIT_CELL 8 8 3 (each cell
 vector should thus be around 30 Angstrom in length), because I read that the cell edge should be at least twice the size of RCUT (i.e., > 24 Angstrom). I tried several versions of the Matsui-Akaogi Force Field to try and get decent values for the cell parameters
 after a cell optimization as was shown in their <a href="https://www.tandfonline.com/doi/abs/10.1080/08927029108022432" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=en-GB&q=https://www.tandfonline.com/doi/abs/10.1080/08927029108022432&source=gmail&ust=1710583154052000&usg=AOvVaw08mefzkva2EbUJ3xIC3KOA">
original paper</a>:<br>
</p></div></div></div></div></div><div lang="en-CH" link="blue" vlink="purple" style="word-wrap:break-word"><div><div><div><div><p class="MsoNormal" style="margin-left:36.0pt"><img border="0" width="1326" height="822" style="width:13.8125in;min-height:8.5625in" src="https://groups.google.com/group/cp2k/attach/2f0ec2fc530cb/MA%20cell%20params.png?part=0.3&view=1"><u></u><u></u></p>
</div></div></div></div></div><div lang="en-CH" link="blue" vlink="purple" style="word-wrap:break-word"><div><div><div>
<div>
<p class="MsoNormal" style="margin-left:36.0pt"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal" style="margin-left:36.0pt">However, no matter which definition I try, I always get cell parameters which are quite a bit larger than those reported and I can't seem to get the original FF definition to work with GENPOT. The parameter sets
 I tried are attached and labeled as follows:<u></u><u></u></p>
</div>
</div></div></div></div><div lang="en-CH" link="blue" vlink="purple" style="word-wrap:break-word"><div><div><div><div><p class="MsoNormal" style="margin-left:36.0pt">- <a href="https://www.tandfonline.com/doi/abs/10.1080/08927029108022432" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=en-GB&q=https://www.tandfonline.com/doi/abs/10.1080/08927029108022432&source=gmail&ust=1710583154052000&usg=AOvVaw08mefzkva2EbUJ3xIC3KOA">
Matsui-Akaogi (MA)</a>: original definition with GENPOT [CPASSERT failed]<br>
</p></div></div></div></div></div><div lang="en-CH" link="blue" vlink="purple" style="word-wrap:break-word"><div><div><div><div><p class="MsoNormal" style="margin-left:36.0pt"><img border="0" width="1337" height="327" style="width:13.927in;min-height:3.4062in" src="https://groups.google.com/group/cp2k/attach/2f0ec2fc530cb/MA%20definition.png?part=0.2&view=1"><u></u><u></u></p>
</div>
<div>
<p class="MsoNormal" style="margin-left:36.0pt"><img border="0" width="1366" height="227" style="width:14.2291in;min-height:2.3645in" src="https://groups.google.com/group/cp2k/attach/2f0ec2fc530cb/MA%20params.png?part=0.1&view=1"></p></div></div></div></div></div><div lang="en-CH" link="blue" vlink="purple" style="word-wrap:break-word"><div><div><div><div><p class="MsoNormal" style="margin-left:36.0pt"><br>
- <a href="https://pubs.acs.org/doi/10.1021/jp050963q" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=en-GB&q=https://pubs.acs.org/doi/10.1021/jp050963q&source=gmail&ust=1710583154052000&usg=AOvVaw0e1aTjw9FkkAi2NxeUjNdW">Naicker (N)</a>: seemingly unaltered WILLIAMS version of the original MA definition, although I'm not 100% sure how to convert the original function into a WILLIAMS type definition. But f(Bi+Bj) does coincide
 with rho_ij and Cij coincides with (Ci+Cj) [a = 3.8865, c = 10.0187]<u></u><u></u></p></div></div></div></div></div><div lang="en-CH" link="blue" vlink="purple" style="word-wrap:break-word"><div><div><div>
<div>
<p class="MsoNormal" style="margin-left:36.0pt">- <a href="https://pubs.acs.org/doi/10.1021/jp037197c" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=en-GB&q=https://pubs.acs.org/doi/10.1021/jp037197c&source=gmail&ust=1710583154052000&usg=AOvVaw1TX0CaZLRAMBbsF5KOpM8I">
Predota (P)</a>: slightly altered version to better simulate TiO2 surfaces [a = 3.8865, c = 10.0187]<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal" style="margin-left:36.0pt">- <a href="https://pubs.acs.org/doi/10.1021/jp207272e" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=en-GB&q=https://pubs.acs.org/doi/10.1021/jp207272e&source=gmail&ust=1710583154052000&usg=AOvVaw1H43-vEJrqk3Bv2_AYFXPG">
Alimommadi (A)</a>: parameters should be the same as Predota but are given in eV instead of kcalmol, but give slightly different result [a = 3.886, c = 10.0177]<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal" style="margin-left:36.0pt"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal" style="margin-left:36.0pt">Do you know why this could be? Is there still something wrong with my settings perhaps or should I just accept that the result might not be ideal? Also, I don't see what I'm doing wrong with the GENPOT definition
 of the original MA FF, so any help would be greatly appreciated!<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal" style="margin-left:36.0pt"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal" style="margin-left:36.0pt">Kind regards,<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal" style="margin-left:36.0pt">Léon<u></u><u></u></p>
</div>
<div>
<div>
<p class="MsoNormal" style="margin-left:36.0pt">On Wednesday 13 March 2024 at 18:26:08 UTC+1 Krack Matthias wrote:<u></u><u></u></p>
</div>
<blockquote style="border:none;border-left:solid #cccccc 1.0pt;padding:0cm 0cm 0cm 6.0pt;margin-left:4.8pt;margin-right:0cm">
<div>
<div>
<p class="MsoNormal" style="margin-left:36.0pt">
<span lang="DE-CH" style="font-size:11.0pt">Hi Léon</span><u></u><u></u></p>
<p class="MsoNormal" style="margin-left:36.0pt">
<span lang="DE-CH" style="font-size:11.0pt"> </span><u></u><u></u></p>
<p class="MsoNormal" style="margin-left:36.0pt">
<span lang="EN-US" style="font-size:11.0pt">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.</span><u></u><u></u></p>
<p class="MsoNormal" style="margin-left:36.0pt">
<span lang="EN-US" style="font-size:11.0pt"> </span><u></u><u></u></p>
<p class="MsoNormal" style="margin-left:36.0pt">
<span lang="EN-US" style="font-size:11.0pt">HTH</span><u></u><u></u></p>
<p class="MsoNormal" style="margin-left:36.0pt">
<span lang="EN-US" style="font-size:11.0pt"> </span><u></u><u></u></p>
<p class="MsoNormal" style="margin-left:36.0pt">
<span lang="EN-US" style="font-size:11.0pt">Matthias</span><u></u><u></u></p>
<p class="MsoNormal" style="margin-left:36.0pt">
<span lang="EN-US" style="font-size:11.0pt"> </span><u></u><u></u></p>
</div>
</div>
<div>
<div>
<div>
<div>
<div style="border:none;border-top:solid #b5c4df 1.0pt;padding:3.0pt 0cm 0cm 0cm">
<p class="MsoNormal" style="margin-bottom:12.0pt;margin-left:72.0pt">
<b><span style="color:black">From: </span></b><span style="color:black"><span><a href data-email-masked rel="nofollow">cp...@googlegroups.com</a></span> <<span><a href data-email-masked rel="nofollow">cp...@googlegroups.com</a></span>> on behalf of Léon Luntadila Lufungula <<span><a href data-email-masked rel="nofollow">Leon.luntad...@uantwerpen.be</a></span>><br>
<b>Date: </b>Wednesday, 13 March 2024 at 12:20<br>
<b>To: </b>cp2k <<span><a href data-email-masked rel="nofollow">cp...@googlegroups.com</a></span>><br>
<b>Subject: </b>[CP2K:20022] Matsui-Akaogi Force Field for bulk TiO2</span><u></u><u></u></p>
</div>
<p class="MsoNormal" style="margin-left:72.0pt">
Dear all,<u></u><u></u></p>
<div>
<p class="MsoNormal" style="margin-left:72.0pt">
 <u></u><u></u></p>
</div>
<div>
<p class="MsoNormal" style="margin-left:72.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" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=en-GB&q=https://pubs.acs.org/doi/10.1021/jp050963q&source=gmail&ust=1710583154052000&usg=AOvVaw0e1aTjw9FkkAi2NxeUjNdW">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. <u></u><u></u></p>
</div>
<div>
<p class="MsoNormal" style="margin-left:72.0pt">
 <u></u><u></u></p>
</div>
<div>
<p class="MsoNormal" style="margin-left:72.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?<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal" style="margin-left:72.0pt">
 <u></u><u></u></p>
</div>
<div>
<p class="MsoNormal" style="margin-left:72.0pt">
Kind regards,<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal" style="margin-left:72.0pt">
Léon<u></u><u></u></p>
</div>
</div>
</div>
</div>
</div>
<div>
<div>
<div>
<div>
<p class="MsoNormal" style="margin-left:72.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
<span><a href data-email-masked rel="nofollow">cp2k+uns...@googlegroups.com</a></span>.<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" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=en-GB&q=https://groups.google.com/d/msgid/cp2k/e185a18f-fcbc-477a-ba8e-95c0191b4ae4n%2540googlegroups.com?utm_medium%3Demail%26utm_source%3Dfooter&source=gmail&ust=1710583154052000&usg=AOvVaw0jBjood8a5nkus8a9TsmxE">
https://groups.google.com/d/msgid/cp2k/e185a18f-fcbc-477a-ba8e-95c0191b4ae4n%40googlegroups.com</a>.<u></u><u></u></p>
</div>
</div>
</div>
</div>
</blockquote>
</div>
</div></div></div></div><div lang="en-CH" link="blue" vlink="purple" style="word-wrap:break-word"><div><div><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 data-email-masked rel="nofollow">cp2k+uns...@googlegroups.com</a>.<br></p></div></div></div></div><div lang="en-CH" link="blue" vlink="purple" style="word-wrap:break-word"><div><div><div><p class="MsoNormal" style="margin-left:36.0pt">
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/9e3e09dd-33ad-448c-bbff-bc163b4b46a0n%40googlegroups.com?utm_medium=email&utm_source=footer" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=en-GB&q=https://groups.google.com/d/msgid/cp2k/9e3e09dd-33ad-448c-bbff-bc163b4b46a0n%2540googlegroups.com?utm_medium%3Demail%26utm_source%3Dfooter&source=gmail&ust=1710583154052000&usg=AOvVaw0zThmgf-3pGB9Lyw2R3tU2">
https://groups.google.com/d/msgid/cp2k/9e3e09dd-33ad-448c-bbff-bc163b4b46a0n%40googlegroups.com</a>.<u></u><u></u></p>
</div>
</div>
</div>
</div>

</blockquote></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/73c2f20e-cfcb-4f3b-bc6b-dcb294acdd57n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/73c2f20e-cfcb-4f3b-bc6b-dcb294acdd57n%40googlegroups.com</a>.<br />