计算化学公社

标题: 利用Gaussian计算的赝势问题 [打印本页]

作者
Author: 让你变成回忆 时间: 2018-9-28 08:26
标题: 利用Gaussian计算的赝势问题
本帖最后由让你变成回忆于 2018-9-28 08:34 编辑

在sob老师的博文《使用Gaussian+PySOC在TDDFT下计算旋轨耦合矩阵元》中，有一个关于Os配合物的实例。实例中TDDFT的计算的关键词为：#p b3lyp/def2TZVP scrf(solvent=CH2Cl2) TD(nstates=10,50-50) 6d 10f nosymm gfinput
但是在输出文件中打印的基组信息中却没有包含Os原子的ECP部分，仅仅只有全电子基组部分。但是我在EMSL的网站上，找到Os的ECP部分在def2TZVP下的定义却是如下：
OS 0OS-ECP 3 60f-ul potential 12 2.76707510 18.90945701 s-ul potential 32 13.87575390 885.40571914 2 6.93787690 25.96704014 2 2.76707510 -18.90945701 p-ul potential 32 10.19379260 320.08390185 2 5.09689620 26.14876493 2 2.76707510 -18.90945701 d-ul potential 32 7.09923846 115.04484313 2 3.54961923 13.62257457 2 2.76707510 -18.90945701
而在Gaussian的输出文件中却没有该字段。另外，我计算了另外一个含有Te原子的体系，基组同样为def2TZVP，由于在Gaussian的输出文件中打印出了ECP信息，因此Pysoc程序无法正确读取基组字段，导致计算失效。
而在sob老师的Os的实例中，计算成功的原因就在于Gaussian的输出文件中不包含ECP字段。

（1）为啥在EMSL上能查到Os的ECP，但是在Gaussian的输出文件中却没有找到？使用的Gaussian版本为G16.A03；测试了一下用G09.E01版本，输出文件中就没有了ECP字段。

作者
Author: sobereva 时间: 2018-9-28 08:49
用genecp的话，gfinput会把ECP和基组信息都打印出来。而直接用def2TZVP关键词的话，那就不一定了（即便此基组对当前的一些元素是赝势基组，自动会加对应的赝势）

作者
Author: 让你变成回忆 时间: 2018-9-28 08:59
本帖最后由让你变成回忆于 2018-9-28 14:21 编辑

sobereva 发表于 2018-9-28 08:49
用genecp的话，gfinput会把ECP和基组信息都打印出来。而直接用def2TZVP关键词的话，那就不一定了（即便此基 ...

sob老师，我刚刚测试了一下。
如果使用G16的话（版本为A03），都会打印出ECP信息；但是如果是G09（版本为E01），就不会打印ECP信息。
所以如果利用PySOC来计算SOC的话，当体系包含有重原子并且考虑ECP的时候，最好就用G09？

作者
Author: sobereva 时间: 2018-9-28 09:03

让你变成回忆发表于 2018-9-28 08:59
sob老师，我刚刚测试了一下。
如果使用G16的话（版本为A03），都会打印出ECP信息；但是如果是G09（版本 ...

这点我没特别留意过，反正要是pysoc因为这个兼容有问题，也就是改几行语句的事，让之能正确识别格式即可

作者
Author: 让你变成回忆 时间: 2018-9-28 09:07
本帖最后由让你变成回忆于 2018-9-28 09:09 编辑

sobereva 发表于 2018-9-28 09:03
这点我没特别留意过，反正要是pysoc因为这个兼容有问题，也就是改几行语句的事，让之能正确识别格式即可

我仔细读过Pysoc的代码，它是在输出文件中通过AO basis 和primitive 这两个关键词之间的语句来定位基函数信息的。因此如果打印出了ECP，基函数信息就包含了ECP，而在pysoc转换的时候，仅仅只有S=1, P=3, SP=4, D=6, F=10这样的对应关系。不知道如果我读基函数的时候，不读ECP（即使包含有ECP），而直接读基组信息部分，这样传给molsoc计算时会不会有影响。因为不了解molsoc计算输出文件的具体格式。

作者
Author: sobereva 时间: 2018-9-28 09:29

让你变成回忆发表于 2018-9-28 09:07
我仔细读过Pysoc的代码，它是在输出文件中通过AO basis 和primitive 这两个关键词之间的语句来定位基函数 ...

应该没影响。反正赝势部分信息也不会被利用

作者
Author: 让你变成回忆 时间: 2018-9-28 09:39

sobereva 发表于 2018-9-28 09:29
应该没影响。反正赝势部分信息也不会被利用

OK，明白了，谢谢sob老师指点。

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