<html xmlns:v="urn:schemas-microsoft-com:vml" 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=utf-8"> <meta name="Generator" content="Microsoft Word 15 (filtered medium)"> <style></style> </head> <body lang="EN-US" link="blue" vlink="purple"> <div class="WordSection1"> <p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D">Dear Natalie<o:p></o:p></span></p> <p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D"><o:p> </o:p></span></p> <p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D">I will try to answer your questions:<o:p></o:p></span></p> <p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D"><o:p> </o:p></span></p> <p class="MsoListParagraph" style="text-indent:-18.0pt;mso-list:l0 level1 lfo1"><![if !supportLists]><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D"><span style="mso-list:Ignore">1)<span style="font:7.0pt "Times New Roman""> </span></span></span><![endif]><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D">Yes, I think that’s right when using an all-electron potential.<o:p></o:p></span></p> <p class="MsoListParagraph" style="text-indent:-18.0pt;mso-list:l0 level1 lfo1"><![if !supportLists]><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D"><span style="mso-list:Ignore">2)<span style="font:7.0pt "Times New Roman""> </span></span></span><![endif]><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D">EPSFIT is the key input parameter to define the split of the basis set into a soft and a hard part. The default value of 1.0E-4 is too large for most systems. I would rather consider values of 1.0E-6 for EPSFIT to ensure that the soft basis includes sufficiently large exponents, e.g. if the largest exponent in the soft basis set for oxygen is only 4 this will give poor results.<o:p></o:p></span></p> <p class="MsoListParagraph" style="text-indent:-18.0pt;mso-list:l0 level1 lfo1"><![if !supportLists]><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D"><span style="mso-list:Ignore">3)<span style="font:7.0pt "Times New Roman""> </span></span></span><![endif]><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D">GAPW splits the basis set into a hard (frozen) and a soft part. This can be applied for all-electron calculations but also for hard pseudopotentials. There is no additional pseudopotential.<o:p></o:p></span></p> <p class="MsoListParagraph" style="text-indent:-18.0pt;mso-list:l0 level1 lfo1"><![if !supportLists]><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D"><span style="mso-list:Ignore">4)<span style="font:7.0pt "Times New Roman""> </span></span></span><![endif]><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D">In a first step, you can use the default values for these input parameters. More import is EPSRHO0. Consider smaller values than the default one for better accuracy. The electron counts printed (use PRINT_LEVEL medium) during and after the SCF procedure can serve as an accuracy check. In addition, the parameters defining the local atomic (Lebedev) grids (see &KIND section) might need larger values in some cases, e.g. for heavier elements.<o:p></o:p></span></p> <p class="MsoListParagraph" style="text-indent:-18.0pt;mso-list:l0 level1 lfo1"><![if !supportLists]><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D"><span style="mso-list:Ignore">5)<span style="font:7.0pt "Times New Roman""> </span></span></span><![endif]><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D">See above (EPSFIT)<o:p></o:p></span></p> <p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D"><o:p> </o:p></span></p> <p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D">HTH<o:p></o:p></span></p> <p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D"><o:p> </o:p></span></p> <p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D">Matthias<o:p></o:p></span></p> <p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D"><o:p> </o:p></span></p> <p class="MsoNormal"><b><span style="font-size:11.0pt;font-family:"Calibri",sans-serif">From:</span></b><span style="font-size:11.0pt;font-family:"Calibri",sans-serif"> cp...@googlegroups.com <...@googlegroups.com> <b>On Behalf Of </b>Natalie Holzwarth<br> <b>Sent:</b> Dienstag, 14. April 2020 17:16<br> <b>To:</b> cp2k <...@googlegroups.com><br> <b>Subject:</b> [CP2K:13128] How to understand/tune parameters for GAPW calculations<o:p></o:p></span></p> <p class="MsoNormal"><o:p> </o:p></p> <div> <div> <p class="MsoNormal">Dear cp2k developers,<o:p></o:p></p> </div> <div> <p class="MsoNormal"> We are very much looking forward to possibly using CP2K particularly using the GAPW scheme to extend our calculations which up to now have been performed using ABINIT and QUANTUM ESPRESSO using PAW. We have tried a number of your tutorials and test inputs, but are have a hard time understanding how to properly pick the parameters. Is there perhaps a writeup somewhere that describes the detailed implementation or perhaps a list of the main subroutines that a user should study? Or perhaps someone might have time to answer some the following silly questions.<o:p></o:p></p> </div> <div> <p class="MsoNormal">1. These questions are based on the notion that we should use the all-electron bases since these would allow n^1 to be represented analogously to Blochl's PAW formulation. Is this the correct way to think about it?<o:p></o:p></p> </div> <div> <p class="MsoNormal">2. According the paper by Lippert et al (1999), the soft wavefunctions and densities are chosen by eliminating the Gaussian basis primitives with too large an exponent. Which of the input parameters determines the soft basis parameters?<o:p></o:p></p> </div> <div> <p class="MsoNormal">3. In addition to the soft density, Blochl's PAW scheme has a local pseudo potential. I am not sure what is the analogous local pseudo potential for GAPW and which parameter controls it. I have assumed that we should use "ALL" potentials in the input files. Is this correct? <o:p></o:p></p> </div> <div> <p class="MsoNormal">4. Some of the examples use the following parameters which seem to be important, but I am not sure how to choose them<o:p></o:p></p> </div> <div> <p class="MsoNormal"><o:p> </o:p></p> </div> <div> <p class="MsoNormal">MAX_RAD_LOCAL<o:p></o:p></p> </div> <div> <p class="MsoNormal">ALPHA0_H<o:p></o:p></p> </div> <div> <p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Arial",sans-serif;color:black">HARD_EXP_RADIUS</span><o:p></o:p></p> </div> <div> <p class="MsoNormal"><o:p> </o:p></p> </div> <div> <p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Arial",sans-serif;color:black">5. Is there a provision/recommendation to use the frozen core approximation and if so, what parameters control that?</span><o:p></o:p></p> </div> <div> <p class="MsoNormal"><o:p> </o:p></p> </div> <div> <p class="MsoNormal">Thanks in advance for any suggestions you can provide and much gratitude for sharing and maintaining the cp2k software package.<br> <br> <o:p></o:p></p> </div> <div> <p class="MsoNormal"><o:p> </o:p></p> </div> <div> <p class="MsoNormal">Sincerely, Natalie Holzwarth<o:p></o:p></p> </div> </div> <p class="MsoNormal">-- <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+unsu...@googlegroups.com">cp2k+unsu...@googlegroups.com</a>.<br> To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/e59f80e3-1dda-4bb4-b969-93b783399ad9%40googlegroups.com?utm_medium=email&utm_source=footer"> https://groups.google.com/d/msgid/cp2k/e59f80e3-1dda-4bb4-b969-93b783399ad9%40googlegroups.com</a>.<o:p></o:p></p> </div> </body> </html>