计算化学公社

标题: OpenMolcas的UHF/ANO-RCC-VDZ能量与其他程序不一致 [打印本页]

作者
Author: zjxitcc 时间: 2020-11-5 21:31
标题: OpenMolcas的UHF/ANO-RCC-VDZ能量与其他程序不一致
本帖最后由 zjxitcc 于 2020-11-6 09:04 编辑

本意是想检验一下各个程序用DKH2的时候Hatree-Fock能量是否一致，结果大部分一致，有一个不一致的

不知道是啥原因，请教大家

体系：1个Mn原子，中性，6重态
基组：ANO-RCC-VDZ，BSE基组库上拷贝的，输入文件用自定义基组。
以下各个程序都加了DKH2关键词，且用的是点核模型，全加了Noysmm。UHF能量如下（a.u.）

全用Cartesian型基函数（6D 10F），OpenMolcas异常
高斯：          -1157.374797
GAMESS：    -1157.374798
OpenMolcas： -1157.373056 （即使用内置ANO-RCC-VDZ基组，结果仍一样）
这里OpenMolcas的初始，我试过传入轨道，或者用程序默认的初猜，结果都是一样的，感觉不像波函数解的问题。

全用球谐型基函数（5D 7F），结果很一致
高斯：          -1157.374597
ORCA：       -1157.374598
OpenMolcas：-1157.374597

是不是还有哪个细节可能没有调整一致？附上6D 10F的输入文件
(, 下载次数 Times of downloads: 31)

作者
Author: beefly 时间: 2020-11-6 17:03
标量相对论计算不建议用Cartesian基函数，会出现数值不稳定，导致Molcas的Cartesian基函数能量比球函数能量更高。Gaussian、Gamess用两种基函数不应该得到相同的能量，应该Cartesian基函数的输入错了

作者
Author: zjxitcc 时间: 2020-11-6 17:08
本帖最后由 zjxitcc 于 2020-11-6 17:09 编辑

beefly 发表于 2020-11-6 17:03
标量相对论计算不建议用Cartesian基函数，会出现数值不稳定，导致Molcas的Cartesian基函数能量比球函数能量 ...

再请教一下“Gaussian、Gamess用两种基函数不应该得到相同的能量，应该Cartesian基函数的输入错了”是啥意思？

高斯用Cartesian型基函数的UHF能量是-1157.374797，高斯用球谐型基函数的UHF能量是-1157.374597，能量是前者低，后者高，看上去是符合直觉的。我这GAMESS只算了Cartesian型基函数，没算球谐型的。您是说这两个数据有误还是？

作者
Author: beefly 时间: 2020-11-7 12:04

zjxitcc 发表于 2020-11-6 17:08
再请教一下“Gaussian、Gamess用两种基函数不应该得到相同的能量，应该Cartesian基函数的输入错了”是啥 ...

两个能量只有小数点后第四位数字不一样，我之前没看出来

作者
Author: zjxitcc 时间: 2020-11-7 13:09

beefly 发表于 2020-11-7 12:04
两个能量只有小数点后第四位数字不一样，我之前没看出来

感谢。那看来全电子相对论计算只能用球谐型基函数算了

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