计算化学公社

标题: Multiwfn对PSI4做电子激发分析产生输入文件求助 [打印本页]

作者
Author: hnile 时间: 2022-3-16 21:29
标题: Multiwfn对PSI4做电子激发分析产生输入文件求助
求教！！！
我希望用mltiwfn做电子激发分析, 要求输入的第一个文件为基函数和分子轨道信息，但我是用psi4做的计算，不知道用psi4怎么输出这个文件？

作者
Author: sobereva 时间: 2022-3-17 03:12
帖子里别没事写那么多叹号，太浮躁

没有叫mltiwfn的程序

写上比如如下内容就会产生fchk文件
energy, wfn = energy('scf', return_wfn=True)
fchk_writer = psi4.FCHKWriter(wfn)
fchk_writer.write('h2o.fchk')

Multiwfn不直接支持从PSI4的输出文件里读入组态系数。若死活非要用PSI4（虽然我看不出有任何理由），就把输出的组态系数信息自己写成Multiwfn支持的通用格式，看Multiwfn手册3.21.A节里General cases: You can also use plain text file as the input file...下面的内容。

作者
Author: zjxitcc 时间: 2022-3-17 09:45
本帖最后由 zjxitcc 于 2022-3-17 09:47 编辑

注意你得说出你用的激发态方法，并不是所有激发态方法的组态系数都能被当前的Multiwfn所分析。常见的CIS/TDHF/TDDFT任务，只要格式合适，当然是没问题的。

另外，PSI4自身就支持产生高斯fch文件，2行即可，以HF/DFT计算为例
scfenergy, scf_wfn = energy('scf', return_wfn=True)
fchk(scf_wfn, 'a.fch')
注意fch文件里是没有激发组态系数的，轨道是基态HF/DFT轨道。

作者
Author: hnile 时间: 2022-3-22 10:42
本帖最后由 hnile 于 2022-3-22 10:45 编辑

sobereva 发表于 2022-3-17 03:12
帖子里别没事写那么多叹号，太浮躁

没有叫mltiwfn的程序

谢谢，以后不加！了。确实浮躁，受教
我试试

作者
Author: hnile 时间: 2022-3-22 10:44

zjxitcc 发表于 2022-3-17 09:45
注意你得说出你用的激发态方法，并不是所有激发态方法的组态系数都能被当前的Multiwfn所分析。常见的CIS/TD ...

谢谢回复，用的格林函数方法，自己写组态系数

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