OpenMolcas的UHF/ANO-RCC-VDZ能量与其他程序不一致

zjxitcc

本意是想检验一下各个程序用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的输入文件
Mn_ANO-RCC-VDZ_6D10F_uhf.zip (31.26 KB, 下载次数 Times of downloads: 31)

2020-11-5 21:39 上传 Uploaded

点击下载
Click to download

beefly

标量相对论计算不建议用Cartesian基函数，会出现数值不稳定，导致Molcas的Cartesian基函数能量比球函数能量更高。Gaussian、Gamess用两种基函数不应该得到相同的能量，应该Cartesian基函数的输入错了

zjxitcc

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

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

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

beefly

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

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

zjxitcc

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

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