计算化学公社

标题: ORCA中使用def2-TZVP计算石墨烯负载Fe3O4团簇分别吸附十二烷烃，菲的单点能时出现报错 [打印本页]

作者
Author: 狍鸮 时间: 2025-6-16 15:24
标题: ORCA中使用def2-TZVP计算石墨烯负载Fe3O4团簇分别吸附十二烷烃，菲的单点能时出现报错
本帖最后由狍鸮于 2025-6-16 15:43 编辑

各位老师好！首先本人是个初学者，使用的电脑配置为CPU为7950X（16核），内存64GB，基本为sob老师推荐的入门级别的配置
该计算我参考了sob老师的两篇帖子genmer：生成团簇初始构型结合molclus做团簇结构搜索的超便捷工具 - 量子化学（Quantum Chemistry） - 计算化学公社和使用量子化学程序基于簇模型计算金属表面吸附问题 - 思想家公社的门口：量子化学·分子模拟·二次元，首先对石墨烯负载Fe3O4团簇模吸附苯酚，菲及十二烷烃进行团簇结构搜索找到能量最低的结构，然后使用ORCA在PBE0-D3（BJ）/def2-SV（P）下进行优化，最终得到了以下结构以及ORCA的输入文件。然后当我使用PBE0-D3（BJ）/def2-TZVP对吸附菲和十二烷烃进行单点能计算时，出现了以下报错。
然后我尝试了很多方法，比如更换泛函，增加单核内存减少并行核数，提升基组等级到def2-QZVP但是还是无法解决。最终我发现使用def2-TZVP(-f)，阉割掉f函数时该报错就解决了。
对此我有几个问题想请教一下各位老师

1.在相同输入结构下，为什么计算吸附苯酚的单点能能成功，而吸附菲以及十二烷烃的单点能却报错
2.如果使用阉割f函数的def2-TZVP或def2-QZVP会减少对该体系Fe3O4团簇的d轨道的描述减少，这样是否会对单点能的计算有影响
3.如果有影响，那么我使用什么方法去解决

本人第一次发帖提问，有任何不妥之处麻烦请各位老师批评指正，再次感谢各位老师

以下是对吸附菲的报错内容
ob aborted:----------[ranks] messagePACKAGE----------------------[0-11] process exited without calling finalizeNumber of atoms ... 170---- error analysis ----- ... 4432Number of shells ... 1636[0-11] on DESKTOP-O222UOT ... 3C:\ORCA_6.0.1\orca_startup_mpi.exe ended prematurely and may have crashed. exit code 0xc0000094RI-J (if used) integral strategy ... SPLIT-RIJ (Revised 2003 algorithm where possible)---- error analysis ----- ... 1Contraction scheme used ... SEGMENTED contractionORCA finished by error termination in Startup.. SCHWARTZCalling Command: mpiexec -np 12 C:\ORCA_6.0.1\orca_startup_mpi.exe D:\ORCA_run\Gr_C104_Fe3O4_3A_phenanthrene_1\Single_point\Gr_C104_Fe3O4_3A_phenanthrene.int.tmp D:\ORCA_run\Gr_C104_Fe3O4_3A_phenanthrene_1\Single_point\Gr_C104_Fe3O4_3A_phenanthrene ... 1.000e-11[file orca_tools/qcmsg.cpp, line 394]: ... 1.000e-11 .... aborting the runn ... NOT USEDExchange Range Separation ... NOT USEDHerobrine@DESKTOP-O222UOT /d/ORCA_run/Gr_C104_Fe3O4_3A_phenanthrene_1/Single_point/Gr_C104_Fe3O4_3A_phenanthreneλ nite Nucleus Model ... NOT USEDCABS basis ... NOT availableAuxiliary Coulomb fitting basis ... AVAILABLE # of basis functions in Aux-J ... 7206 # of shells in Aux-J ... 2234 Maximum angular momentum in Aux-J ... 4Auxiliary J/K fitting basis ... NOT availableAuxiliary Correlation fitting basis ... NOT availableAuxiliary 'external' fitting basis ... NOT availableChecking pre-screening integrals ... done ( 0.0 sec) Dimension = 1636Check shell pair data ... done ( 0.0 sec)Shell pair informationShell pair cut-off parameter TPreSel ... 1.0e-11Total number of shell pairs ... 1339066Shell pairs after pre-screening ... 459079Total number of primitive shell pairs ... 4492820Primitive shell pairs kept ... 855503 la=0 lb=0: 96967 shell pairs la=1 lb=0: 131973 shell pairs la=1 lb=1: 44332 shell pairs la=2 lb=0: 68978 shell pairs la=2 lb=1: 47695 shell pairs la=2 lb=2: 12138 shell pairs la=3 lb=0: 27005 shell pairs la=3 lb=1: 19272 shell pairs la=3 lb=2: 9026 shell pairs la=3 lb=3: 1693 shell pairsCalculating one electron integrals ... done ( 1.2 sec)Calculating RI/J V-Matrix + Cholesky decomp.... done ( 2.2 sec)Calculating Nuclear repulsion ... done ( 0.0 sec) ENN= 42622.465752641627 EhDiagonalization of the overlap matrix:[file orca_tools/qcmsg.cpp, line 394]: .... aborting the run

作者
Author: wzkchem5 时间: 2025-6-16 15:30
检查openmpi设置，检查和openmpi相关的PATH变量、LD_LIBRARY_PATH设置，检查调用orca所用的脚本的正确性（是否正确设置了所有环境变量，sob老师的电脑和你的电脑的所有不同之处你是否都考虑到了）

作者
Author: 狍鸮 时间: 2025-6-16 15:46

wzkchem5 发表于 2025-6-16 15:30
检查openmpi设置，检查和openmpi相关的PATH变量、LD_LIBRARY_PATH设置，检查调用orca所用的脚本的正确性（ ...

谢谢老师帮助，不好意思老师，我第一次发帖操作失误了，刚刚没编写完成，麻烦老师了

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