计算化学公社

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

[ORCA] 请教如何在态平均下让ORCA输出指定态能量

[复制链接 Copy URL]

54

帖子

0

威望

219

eV
积分
273

Level 3 能力者

跳转到指定楼层 Go to specific reply
楼主
本帖最后由 Hugo_314cat 于 2025-11-18 11:50 编辑

老师们好,

我用ORCA程序以SA-CASSCF得到的波函数做QD-NEVPT2计算,输出的"FINAL SINGLE POINT ENERGY"也是基态和各个激发态混合之后的值。
我希望在保持SA-CASSCF这个框架下让"FINAL SINGLE POINT ENERGY"只输出QD-NEVPT2计算后得到的某个指定态的能量,输入文件应当如何写?

谢谢各位老师

四非院校大菜鸟

4279

帖子

4

威望

9473

eV
积分
13832

Level 6 (一方通行)

MOKIT开发者

2#
发表于 Post on 2025-11-18 12:07:05 | 只看该作者 Only view this author
它不是像你想的那样work,而是你需要去ORCA输出文件里看每个根的能量,挑你要的读取。在FINAL SINGLE POINT ENERGY上方有。
自动做多参考态计算的程序MOKIT

54

帖子

0

威望

219

eV
积分
273

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2025-11-18 13:05:27 | 只看该作者 Only view this author
zjxitcc 发表于 2025-11-18 12:07
它不是像你想的那样work,而是你需要去ORCA输出文件里看每个根的能量,挑你要的读取。在FINAL SINGLE POINT ...

我了解您的意思。但是如果他不在末尾输出我所需要的值,我就没法使用这个级别去进行结构优化。
四非院校大菜鸟

4279

帖子

4

威望

9473

eV
积分
13832

Level 6 (一方通行)

MOKIT开发者

4#
发表于 Post on 2025-11-19 09:59:50 | 只看该作者 Only view this author
我记得ORCA的QD-NEVPT2只有能量,没有解析导数,难道你要用数值导数进行结构优化?那不得算到猴年马月
自动做多参考态计算的程序MOKIT

54

帖子

0

威望

219

eV
积分
273

Level 3 能力者

5#
 楼主 Author| 发表于 Post on 2025-11-19 14:12:52 | 只看该作者 Only view this author
zjxitcc 发表于 2025-11-19 09:59
我记得ORCA的QD-NEVPT2只有能量,没有解析导数,难道你要用数值导数进行结构优化?那不得算到猴年马月

是数值导数,优化的是很小的模型体系。我做过速度测试,结合2-zeta基组时的优化速度还是可接受的。
四非院校大菜鸟

4279

帖子

4

威望

9473

eV
积分
13832

Level 6 (一方通行)

MOKIT开发者

6#
发表于 Post on 2026-1-5 10:04:18 | 只看该作者 Only view this author
本帖最后由 zjxitcc 于 2026-1-5 10:14 编辑
Hugo_314cat 发表于 2025-11-18 13:05
我了解您的意思。但是如果他不在末尾输出我所需要的值,我就没法使用这个级别去进行结构优化。

没必要一定要在末尾。你自己写几行脚本提取你要的能量,然后用Gaussian作为optimizer或者使用NEWUOA算法就可以了。

将Gaussian与Grimme的xtb程序联用搜索过渡态、产生IRC、做振动分析
http://sobereva.com/421

无需导数的局部极小化算法NEWUOA在Fortran中的使用简例
http://sobereva.com/536
自动做多参考态计算的程序MOKIT

本版积分规则 Credits rule

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

GMT+8, 2026-1-23 13:11 , Processed in 0.209178 second(s), 20 queries , Gzip On.

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