Thanks Professor! I see the function in ao_utils.F and it makes sense now. <div class="gmail_quote"><div dir="auto" class="gmail_attr">On Thursday, December 15, 2022 at 5:02:36 AM UTC-8 jgh wrote: </div><blockquote class="gmail_quote" style="margin: 0 0 0 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">Hi the maxiumum exponent included in the soft basis is: zet = gauss_exponent(lshell, rc_hard, eps_fit, 1.0_dp) where rc_hard is the "hard" radius of the atom "HARD_EXP_RADIUS". It is the exponent that leads to a functional value of eps_fit at a distance of rc_hard from the atom. regards JH ________________________________________ 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 Nicholas Winner <<a href data-email-masked rel="nofollow">nwi...@berkeley.edu</a>> Sent: Wednesday, December 14, 2022 9:10 PM To: cp2k Subject: Re: [CP2K:18205] How to set multigrid cutoff for GAPW? Thanks Jurg and Marcella. I ran a few tests and it looks like with the defaults all the exponents are pretty small (e.g. <5). How does EPSFIT used to determine the exponent cutoff? I read through some cp2k tutorial slides and all the ones I found say that the exponents are mapped, but not what the criteria exactly is. On Tuesday, December 13, 2022 at 11:51:28 PM UTC-8 jgh wrote: Hi the soft basis set (printed out with MEDIUM) is used for the GPW grid part. Typically, the max exponent of the soft basis is in the range of 5-8. regards JH ________________________________________ 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 Nicholas Winner <<a href data-email-masked rel="nofollow">nwi...@berkeley.edu</a>> Sent: Tuesday, December 13, 2022 11:47 PM To: cp2k Subject: [CP2K:18191] How to set multigrid cutoff for GAPW? I know there is a well defined way for setting CUTOFF in the multigrid for GPW. Following Lippert's paper (<a href="https://link.springer.com/article/10.1007/s002140050523" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=en&q=https://link.springer.com/article/10.1007/s002140050523&source=gmail&ust=1671212951388000&usg=AOvVaw1BjxXLDdGW05qW2XoTdeng">https://link.springer.com/article/10.1007/s002140050523</a>), we can take the relative cutoff and multiply it by the maximum exponent in the basis set. That seems like the wrong strategy for the all electron basis sets, which can have exponents over 10,000. Is there a scheme for setting the cutoff in GAPW, and is there a reference like Lippert's paper that can help understand that scheme? Thanks, Nick -- You received this message because you are subscribed to the Google Groups "cp2k" group. 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>>. To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/cfbee009-eff0-4db2-9f31-630f1242958cn%40googlegroups.com" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=en&q=https://groups.google.com/d/msgid/cp2k/cfbee009-eff0-4db2-9f31-630f1242958cn%2540googlegroups.com&source=gmail&ust=1671212951389000&usg=AOvVaw3gvmYB5tTqtn3AFf_DMSw6">https://groups.google.com/d/msgid/cp2k/cfbee009-eff0-4db2-9f31-630f1242958cn%40googlegroups.com</a><<a href="https://groups.google.com/d/msgid/cp2k/cfbee009-eff0-4db2-9f31-630f1242958cn%40googlegroups.com?utm_medium=email&utm_source=footer" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=en&q=https://groups.google.com/d/msgid/cp2k/cfbee009-eff0-4db2-9f31-630f1242958cn%2540googlegroups.com?utm_medium%3Demail%26utm_source%3Dfooter&source=gmail&ust=1671212951389000&usg=AOvVaw0iRYh05GXsP_beMrFsJcNv">https://groups.google.com/d/msgid/cp2k/cfbee009-eff0-4db2-9f31-630f1242958cn%40googlegroups.com?utm_medium=email&utm_source=footer</a>>. -- You received this message because you are subscribed to the Google Groups "cp2k" group. 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>>. To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/5b86b3b3-4c84-4b5c-9334-68015840f362n%40googlegroups.com" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=en&q=https://groups.google.com/d/msgid/cp2k/5b86b3b3-4c84-4b5c-9334-68015840f362n%2540googlegroups.com&source=gmail&ust=1671212951389000&usg=AOvVaw06RRTgd9EassuxDxkAs_pp">https://groups.google.com/d/msgid/cp2k/5b86b3b3-4c84-4b5c-9334-68015840f362n%40googlegroups.com</a><<a href="https://groups.google.com/d/msgid/cp2k/5b86b3b3-4c84-4b5c-9334-68015840f362n%40googlegroups.com?utm_medium=email&utm_source=footer" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=en&q=https://groups.google.com/d/msgid/cp2k/5b86b3b3-4c84-4b5c-9334-68015840f362n%2540googlegroups.com?utm_medium%3Demail%26utm_source%3Dfooter&source=gmail&ust=1671212951389000&usg=AOvVaw05yhe3yBYbI4J14EkZ4P-u">https://groups.google.com/d/msgid/cp2k/5b86b3b3-4c84-4b5c-9334-68015840f362n%40googlegroups.com?utm_medium=email&utm_source=footer</a>>. </blockquote></div> -- You received this message because you are subscribed to the Google Groups "cp2k" group. 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>. To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/f2938afc-3b0e-497b-a6c2-621eee0f4d26n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/f2938afc-3b0e-497b-a6c2-621eee0f4d26n%40googlegroups.com</a>.