计算化学公社

标题: Gaussian读取ORCA波函数的命令求助 [打印本页]

作者
Author: 好兄弟 时间: 2024-7-3 15:19
标题: Gaussian读取ORCA波函数的命令求助
尊敬的各位老师：
目前Gaussian优化出来的激发态波函数不稳定，想用ORCA的波函数进行计算，请问Gaussian读取ORCA的文件命令应该如何写？
另外Gaussian算激发态的参考态波函数不稳定怎么办？还有其他的解决办法吗？
恳请各位老师不吝赐教，感恩。

作者
Author: zjxitcc 时间: 2024-7-3 15:33
本帖最后由 zjxitcc 于 2024-7-3 15:44 编辑

假设你做的是TDDFT计算，那“优化出来的激发态波函数不稳定”这话是错的，是激发态计算中参考波函数不稳定，也即基态波函数不稳定。不是“激发态波函数不稳定”。解决办法：读取不稳定波函数的轨道，加上关键词nosymm stable=opt将波函数优化至稳定即可，关键词示例

#p UM062X chkbasis nosymm guess=read geom=allcheck stable=opt 其他参数如隐式溶剂模型等

复制代码

这里UM062X是一个示例，如果你用B3LYP就改成UB3LYP。底下不用写坐标，基组，电荷，自旋多重度，都用关键词进行读取了。注意，这里基态与激发态必须使用严格一样的泛函和基组，不能说随便搞一两个泛函/基组读一读。

“还有其他的解决办法吗”此话不成立，因为你连最常用办法也没试过（或者试过，但做了其他不合理操作，导致没有获得想要的结果）。

Gaussian读取ORCA的收敛的轨道，不是在gjf文件里写几个关键词就可以实现。需要使用MOKIT的mkl2fch小程序或molden2fch小程序，从ORCA的mkl或molden文件转化产生fch文件。不过这个问题里你用不到这个。

作者
Author: 好兄弟 时间: 2024-7-4 09:24

zjxitcc 发表于 2024-7-3 15:33
假设你做的是TDDFT计算，那“优化出来的激发态波函数不稳定”这话是错的，是激发态计算中参考波函数不稳定 ...

感谢老师的解答，用UKS计算的T1激发态并用stable测试T1，显示波函数不稳定“The wavefunction has an internal instability.”。
由于原来优化后的t1结构的.chk文件未找到（转为了fchk，原文件可能误删了），于是将原优化后的坐标粘贴在输入文件上，并加关键词“nosymm stable=opt”结果显示仍然内坐标不稳定。
采用内坐标法进行相关更改"opt=z-matrix"也暂未解决此问题。

作者
Author: zjxitcc 时间: 2024-7-4 11:01
本帖最后由 zjxitcc 于 2024-7-4 11:03 编辑

好兄弟发表于 2024-7-4 09:24
感谢老师的解答，用UKS计算的T1激发态并用stable测试T1，显示波函数不稳定“The wavefunction has an int ...

不要自己创造名词，什么“内坐标不稳定”，治好了我的低血压。您在1L和3L说的内容完全不搭边，“采用内坐标法进行相关更改'opt=z-matrix'也暂未解决此问题”已误入歧途，我在2L的建议也没看进去，自顾自乱试。需要立即上传输入、输出文件进行紧急抢救。

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