计算化学公社

标题: 请教如何在态平均下让ORCA输出指定态能量 [打印本页]

作者
Author: Hugo_314cat 时间: 2025-11-18 11:48
标题: 请教如何在态平均下让ORCA输出指定态能量
本帖最后由 Hugo_314cat 于 2025-11-18 11:50 编辑

老师们好，

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

谢谢各位老师

作者
Author: zjxitcc 时间: 2025-11-18 12:07
它不是像你想的那样work，而是你需要去ORCA输出文件里看每个根的能量，挑你要的读取。在FINAL SINGLE POINT ENERGY上方有。

作者
Author: Hugo_314cat 时间: 2025-11-18 13:05

zjxitcc 发表于 2025-11-18 12:07
它不是像你想的那样work，而是你需要去ORCA输出文件里看每个根的能量，挑你要的读取。在FINAL SINGLE POINT ...

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

作者
Author: zjxitcc 时间: 2025-11-19 09:59
我记得ORCA的QD-NEVPT2只有能量，没有解析导数，难道你要用数值导数进行结构优化？那不得算到猴年马月

作者
Author: Hugo_314cat 时间: 2025-11-19 14:12

zjxitcc 发表于 2025-11-19 09:59
我记得ORCA的QD-NEVPT2只有能量，没有解析导数，难道你要用数值导数进行结构优化？那不得算到猴年马月

是数值导数，优化的是很小的模型体系。我做过速度测试，结合2-zeta基组时的优化速度还是可接受的。

作者
Author: zjxitcc 时间: 2026-1-5 10:04
本帖最后由 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

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