计算化学公社

标题: ORCA做TDDFT计算出现orca_cis_mpi.exe出现错误 [打印本页]

作者
Author: calvados80 时间: 2021-3-2 17:43
标题: ORCA做TDDFT计算出现orca_cis_mpi.exe出现错误
如题，本人用ORCA 4.2.0 （win版）做tddft 算一个73原子的体系，当nroots设到200时，总是失败。! BP86 def2-SV(P) def2/J RIJCOSX grid4 gridx4 tightSCF noautostart miniprint nopop
%maxcore  1000
%pal nprocs  4 end
%tddft
nroots 200
TDA false
end
* xyz 2 1
C    -2.35087818 -3.48080345 -2.18119008
C    -2.39266249 -2.21840701 -2.75206614
C    -1.70430181 -1.15991224 -2.13699641
N    -1.00405130 -1.32568883 -1.01729307
C    -1.51887088 -4.97327750 -0.29981614
C    -1.62882695 -3.68771973 -0.98562363

...

以下是错误信息
*** TD-DFT CALCULATION INITIALIZED ***

-------------------
RPA-DIAGONALIZATION
-------------------

Dimension of the (A+B)(A-B) eigenvalue problem ...  62445
Number of roots to be determined             ... 200
Maximum size of the expansion space          ... 6000
Convergence tolerance for the residual       ... 2.500e-07
Convergence tolerance for the energies       ... 2.500e-07
Orthogonality tolerance                      ... 1.000e-14
Constructing the preconditioner             ... o.k.
Building the initial guess                   ...
job aborted:
[ranks] message

[0-3] process exited without calling finalize

---- error analysis -----

[0-3] on NEO-PC
E:\ORCA\ORCApack\orca\orca_cis_mpi.exe ended prematurely and may have crashed. exit code 64

---- error analysis -----

ORCA finished by error termination in CIS
Calling Command: mpiexec -np 4  E:\ORCA\ORCApack\orca\orca_cis_mpi.exe G:\phen\phen-uv-vis\Fe2\BP86-0%\tddft\PBE0.cisinp.tmp G:\phen\phen-uv-vis\Fe2\BP86-0%\tddft\PBE0
[file orca_tools/qcmsg.cpp, line 458]:
  .... aborting the run

做的时候控制1000MB pal core，无论是4核 8核还是12核，nroots 设在150时都能跑出来。测试的泛函包括BP86, TPSSh和 PBE0。但是一上到200个态就不行了，
即便把核控制在4个，1000MB也不行，如上所示。而且我在算BP86这个的时候，全程监视，没有出现爆内存的情况。主要是感兴趣的态波长比较高，感觉没出来，特来讨教是什么原因导致aborting the run。谢谢！

作者
Author: sobereva 时间: 2021-3-2 17:46
看清楚板块再发贴
学术求助帖不要发到公社大厅板块，这次给你移动了，下次再发错将删帖+记违规+扣分处理

作者
Author: calvados80 时间: 2021-3-2 17:48
不好意思，我下次注意

作者
Author: sobereva 时间: 2021-3-2 17:48
找个内存大的机子，把maxcore设大再试。还不行改用Linux版再试

强烈不建议用纯泛函算TDDFT，大概率明显低估激发能。通常首先应当考虑的是PBE0

作者
Author: calvados80 时间: 2021-3-2 18:00
谢谢sob老师，我试试。
我现在在重复文献中的数据，文献用的BP86做的，UV-VIS上550nm左右出峰，对应是2.4eV左右的HOMO-LUMO gap。我用TPSSh结果在2.7eV左右，怕上到PBE0偏的更多。算LUMO-HOMO gap和做TDDFT计算UV-VIS要求的泛函可否不同，我看到您发的帖子，介绍算有机体系LUMO-HOMO gap 用B3LYP，TDDFT用PBE0，只是和我的实验数据差距过大了。

作者
Author: sobereva 时间: 2021-3-2 20:09

calvados80 发表于 2021-3-2 18:00
谢谢sob老师，我试试。
我现在在重复文献中的数据，文献用的BP86做的，UV-VIS上550nm左右出峰，对应是2.4e ...

纯泛函算HOMO-LUMO gap也有明显低估。一般适合用20~25% HF成份的泛函算这个。

PBE0算HOMO-LUMO gap亦可，和B3LYP很接近，因为HF成份只差5%而已。

重复文献数据和方法选择的合理性是两码事，单纯为了重复文献数据BP86用也就用了。

作者
Author: calvados80 时间: 2021-3-3 09:16
好的，谢谢sob老师，我不打算用纯泛函做了

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