CP2K中利用纯泛函计算能量相关疑问

Yjc

最近在学习sob老师cp2k培训班中能量计算相关问题-纯泛函计算部分时有一些自己的理解和疑问。

在cp2k中，纯泛函结合赝势基组的计算任务常采用GPW（Gaussian and Plane Wave Method），而需要进行全电子基组进行运算时常采用GAPW（Gaussian and Augmented Plane Wave Method）来进行；在学习过程中我有以下两个疑问：

①cp2k自带的全电子势文件没有包括第四周期以后的元素，也就是对于较重元素没有定义全电子势，需要自己手写；影响全电子势的主要因素是r_loc，请问这个代表的是什么意思？为什么在定义其他元素的时候开发者都建议其设为0.5呢？对于第五周期的银原子，其全电子势的文件是否应该像下面这么写？
Ag ALLELECTRON ALL
  9  18  20  
   0.5   0

②在GAPW计算时支持全电子基组和赝势基组一起使用，讲义中提到NMR或者X光光谱的计算中我们会将感兴趣的原子使用全电子基组，其他原子选择赝势基组；那么对于具体例子而言哪些种类的原子是在这两种任务中我们比较感兴趣的呢（能否举一些例子）？此外讲义上给出白云石（CaMgC2O6）的计算例子，将C和O用全电子基组，将Ca和Mg使用赝势基组，这又是基于什么考虑？

感谢各位老师的解答！

sobereva

Yjc 发表于 2023-7-10 17:19
谢谢社长的指点，我基本上明白了；可您说的那篇JCP原文我在官网和Google上都没找到诶，能否请您发下网址 ...

J. Chem. Phys., 152, 194103 (2020)

Yjc

sobereva 发表于 2023-7-10 13:47
1 用来划分hard和soft部分。细节看CP2K的2020年JCP原文里的算法介绍，特别是16式
只要是一个大致合适的值 ...

谢谢社长的指点，我基本上明白了；可您说的那篇JCP原文我在官网和Google上都没找到诶，能否请您发下网址，谢谢了！

sobereva

1 用来划分hard和soft部分。细节看CP2K的2020年JCP原文里的算法介绍，特别是16式
只要是一个大致合适的值，能对两部分划分即可，只是数值计算上的处理，不需要太讲究

是

2 你自己做实际研究的时候自然会明白什么原子是自己感兴趣的。诸如算NMR碳谱显然碳得用全电子基组。

CaMgC2O6是随便举个GAPW的例子，让学员知道全电子基组和赝势基组混用的时候输入文件的格式是什么样的。