计算化学公社

标题: 赝势角动量截断与收敛等问题求助【已收敛但截断问题未知】 [打印本页]

作者
Author: hdhxx123 时间: 2021-4-7 20:51
标题: 赝势角动量截断与收敛等问题求助【已收敛但截断问题未知】
本帖最后由 hdhxx123 于 2021-4-11 16:56 编辑

老师们好，在博文“赝势的函数形式以及在量子化学程序中定义的方式
http://sobereva.com/188”中提到可以对赝势进行截断，但是在GAMESS中不知该如何写赝势定义来截断，故求助，提前感谢帮助！
以下是La的SDD在GAMESS-US的定义，报错后不知该如何修改：
1-ECP GEN 28 5
1    ----- h-ul potential -----
0.00000000E+00 2 1.00000000E+00
1    ----- s-h potential -----
5.85201953E+02 2 1.94414180E+01
1    ----- p-h potential -----
3.30109510E+02 2 1.60163530E+01
1    ----- d-h potential -----
1.86058232E+02 2 1.51282590E+01
1    ----- f-h potential -----
-4.94333520E+01 2 2.31038750E+01
1    ----- g-h potential -----
-2.01230200E+01 2 1.56390200E+01

作者
Author: sobereva 时间: 2021-4-8 03:21
(, 下载次数 Times of downloads: 24)
(, 下载次数 Times of downloads: 22)

作者
Author: hdhxx123 时间: 2021-4-8 07:50

sobereva 发表于 2021-4-8 03:21

非常感谢sob！

作者
Author: hdhxx123 时间: 2021-4-8 18:37
本帖最后由 hdhxx123 于 2021-4-11 18:17 编辑

关于赝势的截断Sob老师已经说的十分清楚了，但是很不幸的是我碰到了一些问题导致即使用了这种截断仍然无法进行正常任务的情况。体系仍旧是http://bbs.keinsci.com/thread-22442-1-1.html中的体系，运行语句也是一样的，在此不再赘述。
问题如下：
可能是因为镧系的特殊情况，即使将La的h角动量截断，GAMESS仍然输出报错，目前试出来能正常运行的任务只有将赝势设置为SBKJC。对应了mol
失败的体系有以下三种情况：
1.截断后的SDD赝势。对应了mol_SDD

2.写了PP=READ,$ECP赝势全部设置为NONE的全电子基组SARC-DKH2。对应了mol_SARC_PP
3.删除PP=READ,同时删除了$ECP部分的全电子基组SARC-DKH2。对应了mol_SARC_NOPP

这部分是个人理解失误，目前GKS-EDA不支持使用相对论哈密顿。收敛的问题是因为目前用的是SBKJC的赝势，但mol.inp文件中却没有同时把La的SDD基组改变为SBKJC赝势对应的基组导致的，希望大家在改变赝势时多多注意，记得同时变更基组！截断为何报错仍未知（若将$GUESS中GUESS=HUCKEL改为GUESS=HCORE即能解决赝势的报错问题），但建议大家不要在GAMESS用高于g（即L>4)的赝势基组，可能会导致不收敛等杂问题

其中1、3报错一致，2的报错我没有理解，以上的文件都在“输入输出文件.7z”中，以上的报错有可能是个人的输入文件书写方式有误导致，供各位大佬指正讨论，感谢各位。

另：目前mol.out文件为收敛无望后体系手动ctrl-c中断的情况，该任务是能正常运行的，十分希望各位大佬能帮助解决不收敛的问题！再次提前感谢各位

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