计算化学公社
标题: ORCA做TDDFT计算时出现I/O Operation Failed问题 [打印本页]
作者Author: 11hanche 时间: 2021-3-3 09:34
标题: ORCA做TDDFT计算时出现I/O Operation Failed问题
各位老师好，本人计算小白，刚开始用ORCA做TDDFT计算，遇到“FATAL ERROR ENCOUTERED” "I/O OPERATION FAILED"问题，输入文件如下：
! B3LYP D3 def2-TZVP def2/J def2-TZVP/C RIJCOSX grid4 gridx4 tightSCF opt freq noautostart miniprint nopop
%maxcore 1000
%pal nprocs 8 end
%tddft
nroots 10
TDA false
end
* xyz 0 1
C 0.84484901 5.92903392 -9.34042271
C 0.78991981 7.13917346 -8.71639172
……
*
<span]ATI-Ex.proc6.scfgrad后就没有进行下去了。
<span]请问各位老师这是什么问题？该如何解决？谢谢各位老师！
作者Author: abin 时间: 2021-3-3 14:20
再来一次。
如果还这样，
可以考虑添加或者更换一个硬盘。
如果是多人共用的存储，
可以考虑修改读写目录。
或者还有其他可能，
先重跑一次吧。
另，RIJCOSX，gridx4，你是搞笑的吗？
或者建议，优化和能量计算，分开做。
作者Author: 11hanche 时间: 2021-3-3 15:53
请问RIJCOSX，grid4是有什么问题？开始还不太会写命令，用Multiwfn生成了一个输入文件。
作者Author: abin 时间: 2021-3-3 17:53
看手册，
我就不粘贴复制了。
搜索RIJCOSX就看到了。
作者Author: wzkchem5 时间: 2021-3-3 19:53
GridX4太小了。TDDFT单点可能还可以，但是opt freq肯定不行，容易结构优化不收敛，频率计算也容易有假虚频。
另外你确定你是要做激发态结构优化吗（而不是先做基态结构优化再在基态结构基础上算TDDFT）？如果只优化第一激发态的结构，不需要写nroots 10，一般可以把nroots设成1或者2（取决于第一激发态和第二激发态离得多远），占用硬盘会小很多。此外也可以考虑把基组减小到def2-SVP，对于一般的体系，尤其是def2-TZVP跑不动的那种大体系，def2-SVP对于结构优化来说还算可以接受，最后用triple-zeta基组算个单点就行了。
作者Author: sobereva 时间: 2021-3-3 23:05
如果是量化初学者，建议先从Gaussian上手，遇到问题的几率会小得多。优化第一激发态算最低三个态足矣（这也是Gaussian默认的）。ORCA不像G16那样有TDDFT的解析Hessian，算freq这一步耗时极高
激发态优化用低档次3-zeta基组诸如6-311G*或def-TZVP完全足够，没要带f极化的昂贵的3-zeta基组诸如def2-TZVP。一般def2-SVP也可以接受。
另外，TDDFT计算用PBE0多数情况比用B3LYP更合适
相关信息看
Gaussian中用TDDFT计算激发态和吸收、荧光、磷光光谱的方法
http://sobereva.com/314（http://bbs.keinsci.com/thread-2413-1-1.html）
谈谈量子化学中基组的选择
http://sobereva.com/336（http://bbs.keinsci.com/thread-3545-1-1.html）
简谈量子化学计算中DFT泛函的选择
http://sobereva.com/272（http://bbs.keinsci.com/thread-536-1-1.html）
作者Author: 11hanche 时间: 2021-3-4 08:51
谢谢解答，我重新算下试试。
基态结构优化后算的激发态，是想算吸收和发射光谱，所以先算了10个态。
作者Author: 11hanche 时间: 2021-3-4 09:00
谢谢sob老师，我再试试。
鉴于高斯版权问题，目前只能从开源软件入手，配合购买的MOMAP。是想计算分子的吸收和发射光谱以及其它激发态的性质，请问老师ORCA相比其它开源量化软件是否更适合做激发态的计算？
sob老师提到的三篇文章都看过了，没学到位，我会再详细看看。
作者Author: wzkchem5 时间: 2021-3-4 23:44
本帖最后由 wzkchem5 于 2021-3-4 23:46 编辑
那你输入文件写错了，opt freq和%tddft写在同一个文件里，意思是优化激发态结构然后算激发态频率，而不是做基态结构优化和频率分析之后再做个TDDFT。要优化基态结构然后在优化的结构上算激发态，opt freq和%tddft必须是两个独立的输入文件。又或者也可以用$new_job连接成一个输入文件，但是我个人不太喜欢那样做，容易引起歧义。如果你的意思是你已经把基态优化做完了，现在这个输入文件就是做激发态结构优化的（算发射光谱用），那么输入文件写得也有问题，因为一次只能优化一个态，而且默认优化第一个态，你现在的写法等于优化第一个激发态的结构，但是算10个态，这样会浪费很多计算量。
作者Author: sobereva 时间: 2021-3-5 00:03
ORCA根本不是开源程序
在免费程序里，TDDFT方面ORCA是首选
| 欢迎光临 计算化学公社 (http://bbs.keinsci.com/) |
Powered by Discuz! X3.3 |