计算化学公社

标题: Warnings of single point calculation [打印本页]

作者
Author: fantasticqhl 时间: 2015-2-7 06:07
标题: Warnings of single point calculation
做一个单点能的计算, 输入文件如下:

! SCS-MP2 UHF ECP{SDD=Ce,def2-TZVP/J} def2-TZVP def2-TZVP/J TightSCF NoFrozenCore
%pal nprocs 16
   end
%mp2 MaxCore 5000
      RI  true
      end
* xyz 1 2
  V                2.75373700 0.00001900 0.50261100
  Ce             -0.21277100 -1.71582100 -0.44976500
  O                3.46222500 0.00014900 2.02117700
*
(xyz 文件有删减)

在out文件的最后有下面的一些warning, 任务提交在服务器上面,一直在计算,
只是文件夹下面的文件一直没有在更新了, 期间没有查看过文件,只是copy到
另外一个文件夹下再查看过out文件.

Checking for AutoStart:
The File: Ce2V.gbw exists
Trying to determine its content:
   ... Fine, the file contains a basis set
   ... The file does not contain calculation information - skipping AutoStart
   ... Fine, the file contains calculation information
Warning: Info after ECPs in GBWFile!
   ... Fine, the calculation information was read
   ... Fine, the file contains a basis set
   ... The file does not contain calculation information - skipping AutoStart
   ... Fine, the file contains calculation information

现在不知道这个单点的计算状态, 是已经死了呢还是在继续算呢? 另外这些warning对结果会有影响吗?
因为另外的一个单点又有如下的warning:
WARNING!!!!!!!
Your GBWFile is either corrupt or from a different ORCA version!
Please be VERY careful with your calculation results!!!

Data_size, sizeof( TOrcaInfo ): 443764, 49704
   ... Fine, the calculation information was read
   ... The file does not contain a basis set - skipping AutoStart
   ... Fine, the file contains calculation information

看到这个warning,就更担心结果了! 期间也在google里面搜过,只是没有搜到什么类似的, 想来这个warning不常见?
而在out文件中又有这么一句话:
INFORMATION: There is an auxiliary basis set but the global RI flag is NOT set

是不是在混合基组下,不能使用RI加速呢? 我在另外的一个帖子里面看到一个网友只是讲
他没有测试混合基组下面的RI近似.

问题有点多, 麻烦熟悉的人给提点一下, 十分感谢!

作者
Author: sobereva 时间: 2015-2-7 11:01
! RI-SCS-MP2 UHF ECP{def-TZVP=Ce,def-TZVP/JK} RIJK def2-TZVP def2-TZVP/JK def2-TZVP/C TightSCF NoFrozenCore noautostart

作者
Author: fantasticqhl 时间: 2015-2-7 15:40

sobereva 发表于 2015-2-7 11:01
! RI-SCS-MP2 UHF ECP{def-TZVP=Ce,def-TZVP/JK} RIJK def2-TZVP def2-TZVP/JK def2-TZVP/C TightSCF NoFro ...

谢谢SOB，我再查看一下你帮我改动的关键字！

作者
Author: fantasticqhl 时间: 2015-2-10 16:39

sobereva 发表于 2015-2-7 11:01
! RI-SCS-MP2 UHF ECP{def-TZVP=Ce,def-TZVP/JK} RIJK def2-TZVP def2-TZVP/JK def2-TZVP/C TightSCF NoFro ...

再次请教一下SOB：

我算得这个单点，按照你给我修正后的关键词运算了，但是特别慢，2天多了，输出文件一点没有更新。

体系有2个Ce原子，2个V原子，9个O原子，用上面的方法与基组计算。
计算是用一个节点的16核计算，内存是5G

请问一点没有更新，是由于体系太大算不动吗？如果增加核数，并行的效率怎么样呀？十分感谢！

作者
Author: sobereva 时间: 2015-2-10 17:34
如果你已经把你在QQ上问我的并行的问题正常解决了，随便算一个体系都能正常运行，但依然算不动，那么把基组减小为def2-SVP再看看，并且仔细观察输出文件看卡在哪步了。

作者
Author: fantasticqhl 时间: 2015-2-10 22:06

sobereva 发表于 2015-2-10 17:34
如果你已经把你在QQ上问我的并行的问题正常解决了，随便算一个体系都能正常运行，但依然算不动，那么把基组 ...

谢谢sob！并行应该没有问题，我都用上了绝对路径了！我再测试一些小体系吧！

PS：QQ上面问你的，是我朋友，不是我，嘿嘿！

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