请教对MOLOPT基组格式的理解

cj4566 · 发表于 Post on 2021-9-8 17:54:38

本帖最后由 cj4566 于 2021-9-8 17:56 编辑

这两天在研究怎样用cp2k的OPTIMIZE_BASIS模块优化基组，看了网上的教程，一头雾水，有几个问题想请教大家。

1. 为什么这里H元素的主量子数是2？而且发现所有元素的MOLOPT基组里，主量子数都是2。这里是指所有元素都只考虑最外两层轨道，而与元素的实际主量子数无关吗？
2. 对最大角量子数的数值同样感到困惑。按照官网的例题（https://www.cp2k.org/exercises:2015_cecam_tutorial:basis_set_optimisation_using_optimize_basis），这里的最大角量子数好像跟考虑的基组函数的数量也有关系？
3. 对于Cu元素（3s2 3p6 3d10 4s），现有的DZVP基组如下：
Cu DZVP-MOLOPT-SR-GTH DZVP-MOLOPT-SR-GTH-q11
1
2 0 3 6 2 2 2 1
   5.804051150731  0.020918100390  0.045720931893 -0.004381592772 -0.021109803873  0.275442696345 -0.101811263028 -0.016523760157
   2.947777593081 -0.106208582202 -0.094026024883  0.017185613995  0.020960301873  0.351705110927 -0.207670618594  0.055142365254
   1.271621207972  0.307397740339 -0.110623536813 -0.089805629814  0.233442747472  0.331635969640 -0.224161904198 -0.286656089760
   0.517173767860  0.240805274553 -0.742218346329  0.054415126660  0.369266430953  0.259386540456  0.176105738988 -0.502349311598
   0.198006620331 -0.798718095004  2.208107372713  0.446326740476 -1.405067129701  0.151105835782  0.210534119173 -0.508940020577
   0.061684232135 -0.738671023869 -1.720016262377  0.468516012555  1.042169799071  0.030634833418  0.902456275117  0.682110135764
第三行最后的“2 2 2 1”是指2个s函数、2个p函数、2个d函数和1个f函数吗？这里的q11表明该基组只考虑了11个价电子。但是我现在想考察3s和3p轨道电子对DOS的贡献，因此需要生成考虑了19个价电子的基组。用OPTIMIZE_BASIS进行基组优化的时候，上述一系列参数该怎么设置呢？谢谢！

biogon · 发表于 Post on 2021-9-8 18:14:34

需要考虑19个价电子,为啥不用全电子基组

cj4566 · 发表于 Post on 2021-9-8 18:31:54

biogon 发表于 2021-9-8 18:14
需要考虑19个价电子,为啥不用全电子基组

如果Cu用全电子基组，是不是体系的所有元素都应该用全电子基组呢？我不知道混合基组在输入文件里该怎么描述，有类似的例子吗？谢谢！

biogon · 发表于 Post on 2021-9-9 10:59:43

cj4566 发表于 2021-9-8 18:31
如果Cu用全电子基组，是不是体系的所有元素都应该用全电子基组呢？我不知道混合基组在输入文件里该怎么描 ...

全都用呗

cj4566 · 发表于 Post on 2021-9-9 12:09:13

biogon 发表于 2021-9-9 10:59
全都用呗

谢谢啊。又遇到问题了，我尝试用全电子基组的时候，发现Cs元素可用的基组找不到。在BSE网站下载Cs的def2-TZVP-RIFIT，复制到EMSL_BASIS_SETS文件里，计算的时候竟然报错，显示找不到基组。我觉得另一个方法是复制def2-TZVP参考文献里的基组，但是格式有点不一样，还不知道怎么改。

左图是EMSL_BASIS_SETS文件里的格式，右图是Phys. Chem. Chem. Phys. 7, 3297 (2005)文献里的格式，您知道怎么转换吗？谢谢了！

cj4566 · 发表于 Post on 2021-9-9 21:59:36

这个问题解决了，打扰了

wolfli369 · 发表于 Post on 2021-9-9 22:13:41

cj4566 发表于 2021-9-9 21:59
这个问题解决了，打扰了

请问怎么解决的，能指导一下吗

biogon · 发表于 Post on 2021-9-10 11:05:33

cj4566 发表于 2021-9-9 12:09
谢谢啊。又遇到问题了，我尝试用全电子基组的时候，发现Cs元素可用的基组找不到。在BSE网站下载Cs的def2- ...

建议用pob-TZVP-rev而不是别的什么基组

cj4566 · 发表于 Post on 2021-9-23 10:41:27

biogon 发表于 2021-9-10 11:05
建议用pob-TZVP-rev而不是别的什么基组

好的，谢谢！

cj4566 · 发表于 Post on 2021-9-23 10:43:00

wolfli369 发表于 2021-9-9 22:13
请问怎么解决的，能指导一下吗

在google group找人替我优化的

wolfli369 · 发表于 Post on 2021-9-23 17:09:13

cj4566 发表于 2021-9-23 10:43
在google group找人替我优化的

好的，谢谢解答

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[CP2K] 请教对MOLOPT基组格式的理解

浏览过的版块