计算化学公社

标题: LMO-EDA MP2能量分解内存分配问题 [打印本页]

作者
Author: chunlinxxx 时间: 2019-3-17 09:52
标题: LMO-EDA MP2能量分解内存分配问题
各位好！最近尝试用GAMESS-US对一个大体系(含约100左右轻元素)在MP2/6-31g(d)水平进行LMO-EDA计算，期望得到一些定性的结果。但由于是GAMESS不熟悉，对其内存设置很是头疼，期间也在本论坛搜索了很多相关的帖子，也阅读了Sob老师的帖子：在GAMESS-US中设定内存分配方式 http://bbs.keinsci.com/thread-681-1-1.html，但是由于缺乏计算机硬件常识，设置的内存分配参数总是提示错误。在此请教各位，如果本人可用资源是56核心，1TB物理内存的高性能计算节点，该如何较大的设置以保证计算顺利进行，附上自己的输入文件，谢谢各位了。

空格$BASIS GBASIS=N31 NGAUSS=6 NDFUNC=1 NPFUNC=1 $END
空格$CONTRL SCFTYP=RHF RUNTYP=EDA MPLEVL=2 MAXIT=100 ispher=1 MULT=1 $END
空格$System TIMLIM=99999999 memddi=20000 mwords=600 $END
空格$lmoeda matom(1)=48,48 mcharg(1)=0,0 mmult(1)=1,1 $end
空格$DATA
Title
C1
...(96原子坐标)
空格$END

作者
Author: sobereva 时间: 2019-3-17 10:11
具体是什么体系？如果是分解弱相互作用，用PSI4做SAPT比用GAMESS-US强百倍

作者
Author: chunlinxxx 时间: 2019-3-17 10:15
结构是氢键二聚体，PSI4计算中心无，自己电脑跑不动，同时我比较想得到delta POL项。
分子坐标如下：
P 15.0 -1.91380    0.50660 -0.29630
O    8.0 -2.15560    0.36230    1.31160
O    8.0 -1.31450 -0.79210 -0.74630
N    7.0 -0.95440    1.84680 -0.54290
C    6.0 -3.46740    1.01000 -1.00250
C    6.0 -2.45310    1.49540    2.06100
C    6.0 -1.37330    3.13500 -0.80810
H    1.0    0.01530    1.68630 -0.25670
C    6.0 -3.63370    2.30790 -1.34740
C    6.0 -4.56350    0.02510 -1.08150
C    6.0 -3.77760    1.78610    2.34560
C    6.0 -1.41680    2.29880    2.51260
C    6.0 -2.67510    3.37220 -1.23770
C    6.0 -0.45030    4.22200 -0.65880
H    1.0 -4.59810    2.59460 -1.76150
C    6.0 -5.87930    0.40670 -0.79790
C    6.0 -4.31270 -1.30730 -1.41890
C    6.0 -4.07040    2.91880    3.09480
H    1.0 -4.55360    1.12070    1.98720
C    6.0 -1.72400    3.43200    3.25270
H    1.0 -0.39480    2.02590    2.27550
C    6.0 -3.06830    4.69650 -1.58400
C    6.0 -0.88230    5.52730 -1.00590
C    6.0    0.85670    4.04780 -0.14210
C    6.0 -6.92120 -0.49940 -0.87670
H    1.0 -6.08380    1.42410 -0.48820
C    6.0 -5.35430 -2.22440 -1.48640
H    1.0 -3.30220 -1.62980 -1.63500
C    6.0 -3.04730    3.74470    3.54220
H    1.0 -5.10080    3.15440    3.32950
H    1.0 -0.92430    4.07190    3.60320
C    6.0 -2.20730    5.73760 -1.48270
H    1.0 -4.08230    4.85270 -1.93290
C    6.0    0.01550    6.60740 -0.85730
C    6.0    1.69850    5.11940    0.00330
H    1.0    1.20920    3.06870    0.15880
C    6.0 -6.67130 -1.82480 -1.22120
H    1.0 -7.93490 -0.19380 -0.65490
C    6.0 -5.05960 -3.58860 -1.82110
H    1.0 -3.28060    4.62910    4.12070
H    1.0 -2.51330    6.74130 -1.74920
C    6.0    1.27700    6.41130 -0.36400
H    1.0 -0.32070    7.59980 -1.13350
H    1.0    2.68950    4.97920    0.41510
C    6.0 -7.75450 -2.76050 -1.28950
N    7.0 -4.79680 -4.67630 -2.07710
H    1.0    1.95420    7.24740 -0.24600
N    7.0 -8.63690 -3.49290 -1.34060
H    1.0 -0.02240 -1.68990    0.24880
N    7.0    0.94770 -1.85020    0.53430
P 15.0    1.90600 -0.50930    0.28820
C    6.0    1.36740 -3.13850    0.79830
O    8.0    2.14540 -0.36270 -1.31990
O    8.0    1.30670    0.78850    0.74060
C    6.0    3.46090 -1.01290    0.99140
C    6.0    2.66970 -3.37540    1.22620
C    6.0    0.44460 -4.22570    0.64930
C    6.0    2.44240 -1.49470 -2.07110
C    6.0    3.62820 -2.31090    1.33510
C    6.0    4.55690 -0.02760    1.06910
C    6.0    3.06380 -4.69980    1.57140
C    6.0    0.87740 -5.53100    0.99530
C    6.0 -0.86300 -4.05170    0.13410
C    6.0    3.76680 -1.78470 -2.35750
C    6.0    1.40600 -2.29790 -2.52260
H    1.0    4.59340 -2.59750    1.74740
C    6.0    5.87220 -0.40840    0.78190
C    6.0    4.30640    1.30430    1.40850
C    6.0    2.20300 -5.74110    1.47060
H    1.0    4.07820 -4.85590    1.91910
C    6.0 -0.02020 -6.61140    0.84720
C    6.0 -1.70460 -5.12350 -0.01090
H    1.0 -1.21610 -3.07250 -0.16600
C    6.0    4.05910 -2.91640 -3.10840
H    1.0    4.54290 -1.11950 -1.99900
C    6.0    1.71270 -3.43010 -3.26450
H    1.0    0.38410 -2.02560 -2.28400
C    6.0    6.91390    0.49790    0.85930
H    1.0    6.07610 -1.42540    0.47070
C    6.0    5.34770    2.22170    1.47460
H    1.0    3.29620    1.62610    1.62740
H    1.0    2.50960 -6.74480    1.73640
C    6.0 -1.28230 -6.41540    0.35530
H    1.0    0.31660 -7.60380    1.12250
H    1.0 -2.69610 -4.98350 -0.42170
C    6.0    3.03570 -3.74200 -3.55570
H    1.0    5.08930 -3.15140 -3.34450
H    1.0    0.91280 -4.06980 -3.61480
C    6.0    6.66430    1.82290    1.20590
H    1.0    7.92720    0.19300    0.63470
C    6.0    5.05340    3.58540    1.81150
H    1.0 -1.95940 -7.25170    0.23760
H    1.0    3.26870 -4.62560 -4.13550
C    6.0    7.74730    2.75880    1.27270
N    7.0    4.79080    4.67280    2.06930
N    7.0    8.62960    3.49150    1.32250

作者
Author: sobereva 时间: 2019-3-17 10:22
计算中心没有就在上面装一个。最终你会发现装PSI4、用PSI4算，比你最后折腾一通GAMESS-US省事得多，耗时还更低，结果还更有意义。鼓捣GAMESS-US完全是绕弯路
弱相互作用的能量分解，GAMESS-US已经完全没有意义了

获得极化部分，用Multiwfn结合Gaussian也能直接实现，看Multiwfn手册4.100.8的例子。

作者
Author: chunlinxxx 时间: 2019-3-17 10:26
哦，谢谢Sob老师的指点，我下去学习下PSI4。此体系我若用SAPT0基组怎么选择呢，此分析是不是基组水平要求高

作者
Author: sobereva 时间: 2019-3-17 10:36
仔细看了这篇就明白了
THE JOURNAL OF CHEMICAL PHYSICS 140, 094106 (2014)

作者
Author: chunlinxxx 时间: 2019-3-17 10:36
谢谢老师

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