计算化学公社

标题: Molpro在计算hf时占据轨道的能量为正值是什么原因 [打印本页]

作者
Author: xinting 时间: 2021-11-2 10:09
标题: Molpro在计算hf时占据轨道的能量为正值是什么原因
大家好，我想用Molpro对YbAg双原子分子进行CI计算，目前先进行的是HF计算，选取的活性空间是（6,3,3,2）。根据基组网站刚开始Yb选择的是带DK的基组而Ag选择的是赝势基组，结果输入文件在基组处出错。然后我选择Yb=cc-pwCVTZ-DK3和Ag=cc-pwCVTZ-DK基组，结果可以运行，但是在HF输出轨道能量的时候占据轨道的能量变成正值，非占据轨道能量为负值，如下图所示

(, 下载次数 Times of downloads: 17) (, 下载次数 Times of downloads: 27)
我的问题

1.PP基组和DK基组是不能一块是用的不？为什么？
2.为什么会存在占据轨道能量为正的情况，如何解决？ (, 下载次数 Times of downloads: 1)
感谢各位老师的回答

作者
Author: biogon 时间: 2021-11-2 10:39
算HF后面又说选活性空间干什么？
cc-pvnz-DK基组是为了全电子相对论计算用的，不能和cc-pvnz-pp混用，还有你用相对论哈密顿了吗？

作者
Author: sobereva 时间: 2021-11-2 11:45
要么用专为全电子相对论计算优化的全电子基组结合相对论哈密顿计算，要么用赝势结合赝势基组计算而不考虑相对论哈密顿。你的计算完全莫名其妙，如果没用相对论哈密顿就根本不应当用专为全电子相对论计算优化的DK后缀的基组。

看下文了解点相对论计算相关常识
使用Gaussian做镧系金属配合物的量子化学计算
http://sobereva.com/581（http://bbs.keinsci.com/thread-21274-1-1.html）

作者
Author: xinting 时间: 2021-11-2 15:33

biogon 发表于 2021-11-2 10:39
算HF后面又说选活性空间干什么？
cc-pvnz-DK基组是为了全电子相对论计算用的，不能和cc-pvnz-pp混用，还有 ...

感谢您的解答
那个活性空间在算HF时是不需要的信息，在这里没用
在相对论的基组计算的时候我使用了相对论哈密顿，设置DKROLL=1
具体的输入文件如下：
***ybag
geometry={angstrom;yb;ag,yb,r}
r=8.0
basis={yb=cc-pwCVTZ-DK3;ag=cc-pwCVTZ-DK}
DKROLL=1
dkho=3
hf
---
然后在HF计算结果中就出现了占据轨道能量为正的情况

作者
Author: xinting 时间: 2021-11-2 15:35

sobereva 发表于 2021-11-2 11:45
要么用专为全电子相对论计算优化的全电子基组结合相对论哈密顿计算，要么用赝势结合赝势基组计算而不考虑相 ...

好的，谢谢sob老师解答，我看看您附的链接

作者
Author: xinting 时间: 2021-11-2 21:47
感谢老师们的解答，我找到问题的所在了

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