计算化学公社»论坛首页 › 理论与计算化学 (Theoretical and Computational Chemistry) › 量子化学 (Quantum Chemistry) › Molpro在计算hf时占据轨道的能量为正值是什么原因

Molpro在计算hf时占据轨道的能量为正值是什么原因

查看数: 4081 | 评论数: 5 | 收藏 Add to favorites 0

关灯 | 提示：支持键盘翻页<-左右->

帖子模式

xinting

发布时间: 2021-11-2 10:09

正文摘要:

大家好，我想用Molpro对YbAg双原子分子进行CI计算，目前先进行的是HF计算，选取的活性空间是（6,3,3,2）。根据基组网站刚开始Yb选择的是带DK的基组而Ag选择的是赝势基组，结果输入文件在基组处出错。然后我选择Yb=cc ...

回复 Reply

xinting 发表于 Post on 2021-11-2 21:47:32

感谢老师们的解答，我找到问题的所在了

xinting 发表于 Post on 2021-11-2 15:35:22

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

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

xinting 发表于 Post on 2021-11-2 15:33:55

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计算结果中就出现了占据轨道能量为正的情况

sobereva 发表于 Post on 2021-11-2 11:45:27

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

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

biogon 发表于 Post on 2021-11-2 10:39:21

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

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register