Dear all,<div><br /></div><div>I have been trying my hand at performing GAPW calculations for the last couple of weeks and I have made due with the defaults for the first-, second- and thir-row elements I plan on using during my calculations. However, elements from the fourth row or higher (Ti, Zr and Hf in my case) are a bit harder to get right.</div><div><br /></div><div>I've scoured the mailing list and have done a lot of test calculations, but I'm still unsure how to optimize all of the parameters in a systematic way. Let me quickly sum up what I know about the different parameters and then ask what is still unclear.</div><div><br /></div><div><u><b>QS/EPSFIT and KIND/HARD_EXP_RADIUS:</b></u></div><div>Both of these parameters shift the border between the soft and hard part of the basis set and have an inverse relationship with the "cutoff" exponent. Larger values of either parameter lowers the cutoff exponent and will thus result in a soft basis set that is softer (i.e., lower exponents). I'm quessing HARD_EXP_RADIUS cannot be too large so as to keep the hard density confined. Matthias Krack suggested to lower EPSFIT to 1.0E-6 in <a href="https://groups.google.com/g/cp2k/c/FUvjXq9RvWY/m/2xSxfzuVAgAJ">another post</a> and Jürg Hutter stated <a href="https://groups.google.com/g/cp2k/c/JX8bDi_w478/m/tjLIF_UUAQAJ">here</a> that the max exponent of the soft basis should generally be in the range of 5-8. My questions here are:</div><div><ol><li>When using the defaults, I generally get max exponents between 1-5 for the soft basis. Should I make sure that the max exponents are in the range 5-8?</li><li>Should you start with lowering EPSFIT (global parameter) with the default HARD_EXP_RADIUS values and then optimize the HARD_EXP_RADIUS (kind-specific parameter) for each element in the system?</li><li>Since EPSFIT and HARD_EXP_RADIUS are both connected to the basis set, should you do this optimization for each basis set individually or can you do it for one and then apply it to other basis sets?<br /></li></ol><div><u><b>EPSISO and EPSRHO0:</b></u></div></div><div>Jürg Hutter stated <a href="https://groups.google.com/g/cp2k/c/n70_s_Tbt1A/m/LFlsSmA4BgAJ">here</a> that EPSISO is useful for heavy elements, but did not state whether to lower or increase it to get better results. Matthias Krack also wrote <a href="https://groups.google.com/g/cp2k/c/FUvjXq9RvWY/m/8VwHvHGUAgAJ">here</a> that reducing EPSRHO0 will lead to better accuracy. My questions here are:</div><div><ol><li>I only found these two suggestions of altering these parameters. Is it worth playing around with them or should I keep them at the default?</li><li>If I want to alter them, do I do this after EPSFIT and HARD_EXP_RADIUS and do I just use the electron counts (i.e., 


"Total charge density (r-space)" and "Total Rho_soft + Rho0_soft (g-space)") to check which values give the best results?</li></ol><div><u><b>LEBEDEV_GRID and RADIAL_GRID</b></u></div><div>It has been said many times over in the mailing list to increase these values for heavier elements as <a href="https://groups.google.com/g/cp2k/c/FUvjXq9RvWY/m/8VwHvHGUAgAJ">Matthias Krack</a> said as well in an earlier post. My question here are:</div><div><ol><li>Marcella Iannuzi stated in <a href="https://groups.google.com/g/cp2k/c/IH5k1SdMu5k/m/kZz4QxbupK8J">2014</a> that LEBEDEV_GRID is usually in the range 50-80 and RADIAL_GRID in the range 100-200. Is this still accurate? I've seen Jürg Hutter use a 100/100 combination for titanium in a <a href="https://groups.google.com/g/cp2k/c/6qjk0K5oRqU/m/YL83SHAIBgAJ">previous post of mine</a>.</li><li>How do I know if these values are too small? Do I start with them at 80/200 and then optimize EPSFIT/HARD_EXP_RADIUS and then reduce the LEBEDEV_GRID/RADIAL_GRID values until I see an increase in the electron counts?</li></ol></div><div><u><b>LMAXN0, LMAXN1 and ALPHA0_HARD:</b></u></div></div><div>I feel like I shouldn't mess around with these parameters but I have seen Jürg Hutter define a larger LMAXN0 for Ti in an <a href="https://groups.google.com/g/cp2k/c/6qjk0K5oRqU/m/YL83SHAIBgAJ">earlier post of mine</a> and seen ALPHA0_HARD set to 10 in several occasions, for example in the <a href="https://www.youtube.com/watch?v=L0hKLjvjIFU&t=1754s">GAPW class</a> on youtube from 2019.</div><div><ol><li>Should I leave these at the defaults or try and increase LMAXN0 and set ALPHA0_HARD to 10?</li></ol><div>My apologies for the long post with a lot of questions, but I feel like some of these details are still missing from the documentation and mailing list and I know that a proper answer would be immensely helpful for people like me who want to use GAPW for non-standard (heavy) elements and run into trouble with the defaults.</div><div><br /></div><div>Kind regards,</div><div>Léon </div></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 visit <a href="https://groups.google.com/d/msgid/cp2k/77608165-4fbc-431c-a73f-9b4a6623bc41n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/77608165-4fbc-431c-a73f-9b4a6623bc41n%40googlegroups.com</a>.<br />