Dear Prof. Hutter,<br /><br />Thanks for pointing me into the right direction!<div><br /></div><div>Unfortunately, the defaults are not sufficient for my system as I am getting rather large values for the lines you suggested to check:<br /><br />(Inner SCF loop 1)<br /> Total charge density (r-space): 0.0000427891<br /> Total Rho_soft + Rho0_soft (g-space): 0.0001310539<br />(Inner SCF loop 2)<br /> Total charge density (r-space): 0.0000494000<br /> Total Rho_soft + Rho0_soft (g-space): 0.0001376711</div><div><br /></div><div>I have attached my input file and the short part that was written to the output if you would like to see it in more detail. My system consists of a TiO2 [anatase (101)] slab with organic molecules adsorbed onto the surface, so I suspect that it might be the Ti atoms for which the defaults are insufficient.<br /><div><br /></div><div>I have been looking into GAPW in the regtests, documentation and forum over the last couple of days and from what I understand, the most important parameters to change should be the following:<br />(FORCE_EVAL / SUBSYS / KIND)<br /><ul><li><b>LEBEDEV_GRID</b>: I'm a bit unsure how to get a correct value for this parameter for each atom kind. Does this involve a convergence test for each atom type individually or is there some rule of thumb?</li><li><b>RADIAL_GRID</b>: same as for LEBEDEV_GRID.</li></ul><div>(FORCE_EVAL /DFT / QS)</div></div></div><div><ul><li><b>EPS_DEFAULT</b>: This one is already at 1.0E-12, should I reduce it further?</li><li><b>EPS_FIT</b>: The default is 1.0E-4, but I've read in another post in a <a href="https://groups.google.com/g/cp2k/c/FUvjXq9RvWY/m/8VwHvHGUAgAJ">reply from M. Krack</a> that it is too large for most systems and should be reduced to 1.0E-6. Should I reduce it in my case?</li><li><b>EPS_ISO</b>: I've seen it mentioned a couple of times, but never if the value should be reduced or what to do with it exactly.</li><li><b>EPS_RHO0</b>: According to the <a href="https://groups.google.com/g/cp2k/c/FUvjXq9RvWY/m/8VwHvHGUAgAJ">reply of M. Krack</a> smaller values than the default give better accuracy.</li></ul><div>(FORCE_EVAL / DFT / MGRID)</div></div><div><ul><li><b>CUTOFF</b>: This is generally lower than for a GPW calculation as I understand it but I guess I should perform a convergence test to choose an appropriate value. Alternatively, I could also multiply the relative cutoff by the maximum exponent in the soft basis set (PRINT_LEVEL MEDIUM to get the values right?) as you and Marcella suggested in <a href="https://groups.google.com/g/cp2k/c/JX8bDi_w478/m/okVfMy6EAQAJ">another post</a>.</li><li><b>REL_CUTOFF</b>: I'm a bit unsure about how to choose this as all the convergence schemes I've read about first converge the CUTOFF with a large REL_CUTOFF and then reduce the value of REL_CUTOFF to get the best balance between accuracy and efficiency, but as explained above, the rule of thumb for the cutoff uses the value of REL_CUTOFF to determine the CUTOFF...</li></ul><div>Any suggestions on how to optimize these parameters would be greatly appreciated!</div></div><div><br /></div><div>I also have two other quick questions about your reply:</div><div><ol><li>I thought that GPW was used for pseudopotential calculations and GAPW was used for all electron calculations, but now you say that you can also do pseudopotential calculations with GAPW. Is there an advantage to using GAPW over GPW? I know that the disadvantage is that some features are not present for GAPW calculations, but is there a reason to use GPW if all your intended features are implemented for GAPW?</li><li>The BASIS_MOLOPT_UZH file contains SVP,- TZVPP- and QZVPP-MOLOPT-GGA-ae basis sets. Is there also a double zeta all electron (MOLOPT) basis set somewhere I can use? I have always been told that single zetas are not intended for production runs, but I suspect that TZVPP is already quite computationally expensive for a pre-optimization.</li></ol><div>My apologies for the barrage of questions, but I don't have anyone in my group with experience in CP2K so it's hard sometimes to figure these things out on my own...</div></div><div><br /></div><div>Thanks in advance and for your previous reply!</div><div><br /></div><div>Kind regards,</div><div>Léon</div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Tuesday, 14 November 2023 at 10:51:05 UTC+1 Jürg Hutter 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;">Hi
<br>
<br>For a pseudopotential calculation you can just change "METHOD GPW" to
<br>"METHOD GAPW".
<br>If you want to do all electron calculations, you also need to use corresponding
<br>basis sets and potentials, e.g.
<br> BASIS_SET SVP-MOLOPT-GGA-ae
<br> POTENTIAL ALL
<br>Use the default settings first, check the lines
<br> Total charge density (r-space): -0.0000000683
<br> Total Rho_soft + Rho0_soft (g-space): -0.0000000043
<br>in the output. If those values are larger than 0.00001 you need probably to
<br>change defaults.
<br>See the regtests and manual for more information.
<br>Come back here for if this basic setup fails.
<br>
<br>regards
<br>JH
<br>
<br>________________________________________
<br>From: <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>Sent: Wednesday, November 8, 2023 5:21 PM
<br>To: cp2k
<br>Subject: [CP2K:19481] Switching from GPW to GAPW
<br>
<br>Dear all,
<br>
<br>I want to perform an AIM analysis with Critic2 which (for as far as I know) requires an all-electron density (i.e. not the density obtained from a GPW calculation) and also plan to calculate NMR chemical shifts in the future which will also require the all-electron GAPW method. Currently, I have figured out how to do calculations with the GPW method but I'm still a bit unsure about how to go from a GPW to a GAPW calculation. Is there someone who could elaborate on what to look at when modifying my input file? I've looked into it a couple of times, but it's a lot of different parameters and I can't seem to figure out how it works...
<br>
<br>All the best,
<br>Léon
<br>
<br>
<br>--
<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><mailto:<a href data-email-masked rel="nofollow">cp2k+uns...@googlegroups.com</a>>.
<br>To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/670d2377-8c38-4444-9130-14ebd3d32d09n%40googlegroups.com" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=en-GB&q=https://groups.google.com/d/msgid/cp2k/670d2377-8c38-4444-9130-14ebd3d32d09n%2540googlegroups.com&source=gmail&ust=1700043596978000&usg=AOvVaw1Q8NmLFZF_6S5Jy1-ZHCDs">https://groups.google.com/d/msgid/cp2k/670d2377-8c38-4444-9130-14ebd3d32d09n%40googlegroups.com</a><<a href="https://groups.google.com/d/msgid/cp2k/670d2377-8c38-4444-9130-14ebd3d32d09n%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/670d2377-8c38-4444-9130-14ebd3d32d09n%2540googlegroups.com?utm_medium%3Demail%26utm_source%3Dfooter&source=gmail&ust=1700043596978000&usg=AOvVaw1MLVliAZkiW9Prmxoa-L__">https://groups.google.com/d/msgid/cp2k/670d2377-8c38-4444-9130-14ebd3d32d09n%40googlegroups.com?utm_medium=email&utm_source=footer</a>>.
<br></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/36846b0d-2a8d-4ce1-a35e-991b76854170n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/36846b0d-2a8d-4ce1-a35e-991b76854170n%40googlegroups.com</a>.<br />