计算化学公社

 找回密码 Forget password
 注册 Register
Views: 622|回复 Reply: 3
打印 Print 上一主题 Last thread 下一主题 Next thread

[ORCA] orca交换轨道不成功

[复制链接 Copy URL]

50

帖子

0

威望

251

eV
积分
301

Level 3 能力者

跳转到指定楼层 Go to specific reply
楼主
各位老师好,请问如何解决ORCA交换轨道不成功的问题(任务显示正常结束,但是输出文件中显示轨道没有交换)。以下是我的输入文件:

np_2_3i_pbe_alter.inp

668 Bytes, 下载次数 Times of downloads: 13

3806

帖子

4

威望

7996

eV
积分
11882

Level 6 (一方通行)

MOKIT开发者

2#
发表于 Post on 2024-7-1 11:11:47 | 只看该作者 Only view this author
如果你的任务是正常结束,且SCF迭代过程不是1圈收敛,是从一个新的能量 SCF迭代了几十圈(甚至更多)才收敛,那说明轨道调换成功了,只不过你调换产生的是SCF初猜,而SCF迭代过程有可能又优化回原来的轨道,迭代过程不是人为可控的。如果你想寻找不同的ROHF/RODFT解,应当尝试其他手段。如果有更多问题或想获得更多帮助,最好上传输入、输出文件。
自动做多参考态计算的程序MOKIT

50

帖子

0

威望

251

eV
积分
301

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2024-7-11 09:32:36 | 只看该作者 Only view this author
zjxitcc 发表于 2024-7-1 11:11
如果你的任务是正常结束,且SCF迭代过程不是1圈收敛,是从一个新的能量 SCF迭代了几十圈(甚至更多)才收敛 ...

谢谢老师的回复,确实是SCF迭代过程中又优化回原来的轨道了,但是我计算的金属的轨道的电子可能占据在d轨道上,也可能占据在f轨道上,所以我计算CAS的时候想包括所有的f轨道和一个d轨道,但是无论我如何交换,d轨道都没法交换到活性空间,请问如何解决呀?orca有没有类似于gaussian中guess=read, alter, only这样的关键词呢

3806

帖子

4

威望

7996

eV
积分
11882

Level 6 (一方通行)

MOKIT开发者

4#
发表于 Post on 2024-7-11 13:28:43 | 只看该作者 Only view this author
本帖最后由 zjxitcc 于 2024-7-11 13:31 编辑
CLJJ 发表于 2024-7-11 09:32
谢谢老师的回复,确实是SCF迭代过程中又优化回原来的轨道了,但是我计算的金属的轨道的电子可能占据在d轨 ...

尝试加上NoIter关键词,不过这个也只是碰运气。因为即使你做完了RODFT计算,成功交换了轨道(NoIter应该可以保证轨道不动),但你在CASSCF步骤读入轨道后,CASSCF本身也是轨道优化过程,能否保持你的初始轨道 又是另一码事。如果到时候你想要的d/f轨道又被旋转出去了、不在活性空间,那不是软件功能不全 导致的,而是用户的主观意愿与客观事实之间的矛盾导致的,关于这点可以看这个帖子的讨论http://bbs.keinsci.com/forum.php?mod=viewthread&tid=46834

另外,我注意到你似乎在使用RODFT轨道做CASSCF初始轨道,并通过手动调换轨道构造活性空间初猜。这十分原始和传统,非常不推荐。建议使用MOKIT自动做CASSCF计算,看看自动给出的活性轨道是否符合期望;如果不符合,可以用计算过程中给出的GVB轨道fch文件自行挑选轨道,GVB轨道呈空间局域,具有成键、反键轨道的特征,挑选起来很简单,不会有歧义(这个轨道像什么、可能是什么、怎么这么离域)。
自动做多参考态计算的程序MOKIT

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2024-11-24 10:52 , Processed in 0.245669 second(s), 24 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list