利用Gaussian计算的赝势问题

让你变成回忆

在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字段。

让你变成回忆

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

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

sobereva

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

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

让你变成回忆

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计算输出文件的具体格式。

sobereva

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

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

让你变成回忆

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

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

sobereva

用genecp的话，gfinput会把ECP和基组信息都打印出来。而直接用def2TZVP关键词的话，那就不一定了（即便此基组对当前的一些元素是赝势基组，自动会加对应的赝势）