计算化学公社

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

[ORCA] orca画出自然轨道用于CASSCF的初猜出现的一些问题

[复制链接 Copy URL]

284

帖子

0

威望

949

eV
积分
1233

Level 4 (黑子)

跳转到指定楼层 Go to specific reply
楼主
本帖最后由 BinWang 于 2025-12-27 19:37 编辑

各位老师:
      大家好!我用orca计算一个双自由基,并打算用CASSCF方法计算能量,出现如下问题:
       根据Sob老师的帖子http://sobereva.com/264,三重态下先做UHF计算之后用guess=(naturalorbitals,read,save,only)生成的自然轨道,相应的alph-HOMO和alpa-HOMO-1轨道恰好就是两个双自由基所在的轨道。和自旋轨道吻合的非常好。并且该方法生成的自然轨道相同原子序号的alpha和beta轨道形状完全一样,非常适用于CASSCF计算。一旦去掉only,则和正则轨道没啥区别,很难用来选择CASSCF的活化轨道。
      我的问题是在orca里面如何实现这一点?我用
%scf                        
UHFNO true
end
结果发现输出的波函数加这一行与不加没啥区别,并不是物理意义很明显的自然轨道。这个该怎么解决?我准备做包含x2c相对论效应的CASSCF计算,所以筛选出合理的活性轨道至关重要,我的体系很大,计算出自然轨道非常重要。


1万

帖子

0

威望

1万

eV
积分
22876

Level 6 (一方通行)

2#
发表于 Post on 2025-12-27 21:54:08 | 只看该作者 Only view this author
做UNO计算,然后读.uno文件里的轨道。.gbw文件里仍然是正则轨道
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php,墙内下载ORCA:https://www.faccts.de/customer)开发团队成员
Google Scholar: https://scholar.google.com/citations?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

284

帖子

0

威望

949

eV
积分
1233

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2025-12-28 15:32:48 | 只看该作者 Only view this author
wzkchem5 发表于 2025-12-27 21:54
做UNO计算,然后读.uno文件里的轨道。.gbw文件里仍然是正则轨道

非常感谢王老师

1

帖子

0

威望

55

eV
积分
56

Level 2 能力者

4#
发表于 Post on 2026-6-27 20:56:52 | 只看该作者 Only view this author
wzkchem5 发表于 2025-12-27 21:54
做UNO计算,然后读.uno文件里的轨道。.gbw文件里仍然是正则轨道

是直接用uno的。uno文件就行了嘛王老师,我也在做类似的计算,用uno.gbw的结果也不合理

4370

帖子

4

威望

1万

eV
积分
14551

Level 6 (一方通行)

MOKIT开发者

5#
发表于 Post on 2026-6-27 21:45:32 | 只看该作者 Only view this author
mourning 发表于 2026-6-27 20:56
是直接用uno的。uno文件就行了嘛王老师,我也在做类似的计算,用uno.gbw的结果也不合理

你可以展示你的体系,具体问题具体提问
自动做多参考态计算的程序MOKIT

1万

帖子

0

威望

1万

eV
积分
22876

Level 6 (一方通行)

6#
发表于 Post on 2026-6-29 17:01:23 | 只看该作者 Only view this author
mourning 发表于 2026-6-27 20:56
是直接用uno的。uno文件就行了嘛王老师,我也在做类似的计算,用uno.gbw的结果也不合理

不合理的具体表现是什么?你的问题和楼主的问题未必是一个原因
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php,墙内下载ORCA:https://www.faccts.de/customer)开发团队成员
Google Scholar: https://scholar.google.com/citations?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

本版积分规则 Credits rule

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

GMT+8, 2026-7-24 15:13 , Processed in 0.534586 second(s), 20 queries , Gzip On.

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