计算化学公社

标题: 求助：SDD基组网站看不懂基组不会用 [打印本页]

作者
Author: liaorongbao 时间: 2024-4-5 10:20
标题: 求助：SDD基组网站看不懂基组不会用
本帖最后由 liaorongbao 于 2024-4-5 10:49 编辑

打开这个SDD基组网站https://www.tc.uni-koeln.de/PP/clickpse.en.html。
我想找到Zn原子的Gaussian格式SDD基组。在Data format菜单后边选择Gaussian，然后在周期表中点击Zn原子。然后自动弹出一个窗口。
这个窗口显示Zn可选用的赝势基组一共有如下四个。
ECP10MDF
ECP10MHF
ECP10MWB
ECP28MWB
ECP28SDF
以上四种中我选择第一个，即ECP10MDF。这个符号的意思是Zn的内层的10个电子将会被赝势代替。那么价层还剩下20个价电子，依然需要用高斯基组表达。MDF中的DF和M的解释如下。
http://sobereva.com/60博文中说DF=Dirac-Fock全相对论，M的含义是考虑实际原子全部价电子的能量。相应地S是只考虑模型体系仅有的单个价电子的能量。我的这个Zn的模型，价电子有20个，而不是仅有1个，因此对应的字母是M。
点开第一个即ECP10MDF，弹出如下信息。
Pseudopotential ECP10MDF for Zn!
Zn 0
ECP10MDF 4 10
G-Komponente
1
2 1.000000 0.000000
S-G
2
2 34.174001 399.986399
2 14.456371 85.489750
P-G
4
2 39.888683 92.381077
2 39.655017 184.771176
2 15.290546 23.002541
2 14.903524 46.057427
D-G
4
2 43.708296 -13.690734
2 43.698536 -20.543980
2 15.150718 -1.316154
2 15.282441 -1.838715
F-G
2
2 8.160014 -0.370360
2 12.228422 -1.062943
这些信息包含了10个内核电子的赝势信息。但却不知道另外的20个价电子的基组信息在哪里？

作者
Author: paramecium86 时间: 2024-4-5 11:21
点开这个赝势之后往下继续翻一下会有 Available basis sets for Zn xxxxxxx 然后又各种配套基组可以选。

作者
Author: zjxitcc 时间: 2024-4-5 14:49
Gaussian内置的SDD赝势基组不满足你的需求？

作者
Author: sobereva 时间: 2024-4-5 17:31
目前的BSE上的SDD并不比Gaussian自带的更新、更全，没必要刻意用这上的

作者
Author: liaorongbao 时间: 2024-4-5 17:52

zjxitcc 发表于 2024-4-5 14:49
Gaussian内置的SDD赝势基组不满足你的需求？

之前听说直接写SDD有可能不如网上现找的新，但后边卢老师说了，高斯自带的也很新

作者
Author: 好多于 时间: 2024-5-29 09:25

sobereva 发表于 2024-4-5 17:31
目前的BSE上的SDD并不比Gaussian自带的更新、更全，没必要刻意用这上的

#p opt TD CAM-B3LYP/genECP iop(9/40=5) pop=full guess=read geom=check

Title Card Required

0 3

C H N 0
6-311+G**
****
Ni 0
SDD
****

Ni 0
SDD

请问卢老师，对于含Ni的体系，使用以上关键词会调用gaussian内置的SDD球谐型基函数是不，不用指定MDF了吧

作者
Author: sobereva 时间: 2024-5-30 06:03

好多于发表于 2024-5-29 09:25
#p opt TD CAM-B3LYP/genECP iop(9/40=5) pop=full guess=read geom=check

Title Card Required

用gen、genecp的时候总会用球谐型基函数
当前写法可以。不建议写pop=full，没实际意义，还导致输出文件极大。要看什么轨道的展开系数，应当Multiwfn载入fch文件，用主功能6的子功能4

作者
Author: 好多于 时间: 2024-5-30 09:46

sobereva 发表于 2024-5-30 06:03
用gen、genecp的时候总会用球谐型基函数
当前写法可以。不建议写pop=full，没实际意义，还导致输出文件 ...

好的谢谢！

欢迎光临计算化学公社 (http://bbs.keinsci.com/)