应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 【求助】BSSE计算出错
201/2返回列表12下一页
查看: 2238  |  回复: 19
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

aaq2800

铁杆木虫 (知名作家)

[交流] 【求助】BSSE计算出错已有8人参与

对gauss不熟悉，一直用ms，审搞人希望用gauss算算，但计算老是出错，请大人们帮我看看问题在哪。吸附的H2分子，弱作用，所以用了bsse。但是总是计算了好久了才出错。好象是说自旋多重度出错了，但不知道怎么处理。
输入文件
% mem=196mb  
#MP2/6-31G(d,p) scf=(maxcyc=600,xqc,conver=6) Counterpoise=2 nosymm

Opt job

0 2                                          
Al            2.848013            1.332099            1.303473    1
Al            2.288651            1.258135           -1.450525    1
Al            3.056950           -1.251095            2.160476    1
Al           -0.126431           -0.019232           -1.845458    1
Al            2.680607           -2.999411           -0.004668    1
Al            0.121452            1.668301            0.332793    1
Al            0.074777           -2.655518           -0.982495    1
Al            0.602885           -2.549631            1.895398    1
Al           -1.119792           -0.733932            0.665542    1
Al            4.042832           -0.587847           -0.326286    1
Al            2.240691           -1.358932           -2.249606    1
Al            0.712283            0.159701            2.553704    1
Si            1.450164           -0.643121            0.166434    1
Sc            2.095926            3.813722           -0.012795    1
H             2.533942            3.512321           -1.950298    2
H             2.496478            4.375593           -1.906499    2
H             3.391091            3.739833            1.559429    2
H             3.246425            4.567051            1.488878    2
H             0.140955            4.302727            0.617929    2
H             0.536578            4.984349            0.793110    2
H             2.381207            5.791306            0.233175    2
H             2.362714            5.729223           -0.599320    2

输出出错信息：
**** Warning!!: The largest beta MO coefficient is  0.10518428D+02

Disk-based method using OVN memory for 22 occupieds at a time.
Estimated scratch disk usage=   618143568 words.
Actual    scratch disk usage=   600654428 words.
JobTyp=2 Pass  1:  I=   1 to  22 NPSUse=  1 ParTrn=F ParDer=F DoDerP=F.
Actual    scratch disk usage=   603252048 words.
JobTyp=3 Pass  1:  I=   1 to  21 NPSUse=  1 ParTrn=F ParDer=F DoDerP=F.
Spin components of T(2) and E(2):
     alpha-alpha T2 =       0.2047990010D+00 E2=     -0.1883034873D+00
     alpha-beta  T2 =       0.7270600968D+00 E2=     -0.7519783603D+00
     beta-beta   T2 =       0.1971814433D+00 E2=     -0.1817972823D+00
(S**2,0)=  0.30869D+01           (S**2,1)=  0.28188D+01
E(PUHF)=      -0.39515486005D+04        E(PMP2)=      -0.39526694195D+04
ANorm=    0.1459123210D+01
E2 =    -0.1122079130D+01 EUMP2 =    -0.39525324250944D+04

。。。。。。。。。。。。。。。。。。。。。。。。
。。。。。。。。中间省略
。。。。。。。。。。。。。。。。。。。。。。。

Mulliken atomic spin densities:
              1
     1  Al   0.030663
     2  Al  -0.074323
     3  Al   0.322903
     4  Al  -0.399333
     5  Al   0.404929
     6  Al  -0.045667
     7  Al  -0.506348
     8  Al   0.398207
     9  Al  -0.312802
    10  Al   0.330209
    11  Al  -0.300362
    12  Al   0.361376
    13  Si   0.020553
    14  Sc   0.725689
    15  H    0.027951
    16  H   -0.001400
    17  H   -0.048594
    18  H    0.014630
    19  H    0.010623
    20  H    0.005703
    21  H    0.019521
    22  H    0.015869
Sum of Mulliken spin densities=   1.00000
Electronic spatial extent (au):  =  7314.8254
Charge=     0.0000 electrons
Dipole moment (field-independent basis, Debye):
    X=     0.5966    Y=     2.4132    Z=     0.0430  Tot=     2.4862
Quadrupole moment (field-independent basis, Debye-Ang):
   XX=  -189.9423   YY=  -188.9083   ZZ=  -191.9076
   XY=     3.2793   XZ=    -1.0367   YZ=     0.4921
Traceless Quadrupole moment (field-independent basis, Debye-Ang):
   XX=     0.3104   YY=     1.3444   ZZ=    -1.6549
   XY=     3.2793   XZ=    -1.0367   YZ=     0.4921
Octapole moment (field-independent basis, Debye-Ang**2):
  XXX=  -869.2207  YYY=   108.7606  ZZZ=   -84.0144  XYY=  -276.4596
  XXY=    34.5194  XXZ=   -30.9336  XZZ=  -291.8271  YZZ=    28.5548
  YYZ=   -26.1225  XYZ=     2.2026
Hexadecapole moment (field-independent basis, Debye-Ang**3):
XXXX= -5348.8369 YYYY= -5256.9854 ZZZZ= -2701.7338 XXXY=   -37.0818
XXXZ=  -145.0473 YYYX=   -26.9552 YYYZ=    77.7606 ZZZX=  -131.6933
ZZZY=    71.0272 XXYY= -1722.1014 XXZZ= -1350.0786 YYZZ= -1318.9702
XXYZ=    19.4781 YYXZ=   -42.4167 ZZXY=   -30.1793
N-N= 2.507531301254D+03 E-N=-1.444427151739D+04  KE= 3.950342541279D+03
Counterpoise: doing DCBS calculation for fragment   2 NewBq=T
Basis read from rwf:  (6D, 7F)
The combination of multiplicity 2 and     8 electrons is impossible.
Error termination via Lnk1e in C:\G03W\l301.exe at Sun Apr 04 18:45:47 2010.
Job cpu time:  0 days 21 hours 19 minutes 30.0 seconds.
File lengths (MBytes):  RWF=   5455 Int=      0 D2E=      0 Chk=      4 Scr=      1
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐，对您同样有帮助:
1楼2010-04-04 20:32:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wanfei

金虫 (小有名气)

#MP2

aaq2800(金币+2):谢谢吹捧，呵呵 2010-04-05 08:26
#MP2要空格# MP2
回复此楼
2楼2010-04-04 21:49:10
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

aaq2800

铁杆木虫 (知名作家)
输入的自旋多重度为2应该没错啊，其他地方有哪里不对呢？
赞一下 回复此楼
3楼2010-04-04 21:49:29
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

aaq2800

铁杆木虫 (知名作家)
引用回帖:
Originally posted by wanfei at 2010-04-04 21:49:10:
#MP2要空格# MP2

应该不会是这个问题把，看最后出错的信息，算着算着自旋多重度出问题，都已经算了20多个小时了

[ Last edited by aaq2800 on 2010-4-4 at 21:52 ]
赞一下 回复此楼
4楼2010-04-04 21:50:46
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangmt

至尊木虫 (著名写手)

我叫MT
★ ★
aaq2800(金币+2):谢谢，昨天晚上回去想了想，应该就是这里出错了， 2010-04-05 08:24
aylayl08(金币+2):谢谢专家指教 2010-04-05 20:47
去看高斯中文手册，66页关于counterpoise的说明和例子，你这个题目要给3套电荷和多重度，应该写成：

0 2 0 2 0 1
Al            2.848013            1.332099            1.303473    1
Al            2.288651            1.258135           -1.450525    1
…………
第一个0 2是指整个体系的电荷与多重度，
第二个0 2是指后面标注1的那部分的电荷与多重度，
第三套0 1指后面标注2的那部分的电荷与多重度，在你的体系中就是8个H
8个H在电荷为0的时候，是不可能有2重态的，所以出错。
赞一下(1人) 回复此楼
一群自以为正义凛然的年轻人将一切不能以科学解释的事情定性为封建迷信并大刀阔斧地进行消灭，其实这是修养不足学识浅薄的一种体现，也是可恶的偏执和愚蠢的自以
5楼2010-04-04 23:22:58
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangmt

至尊木虫 (著名写手)

我叫MT
aaq2800(金币+2):谢谢交流 2010-04-05 08:25
顺便多说一句，你的体系12个Al，外加1个Si，1个Sc，还用6-31g**基组，而且还是mp2计算，看来你得需要几台好机器了——好在是单点算，要是优化，估计你得算到毕业才能算完——不过你竟然还是在windows下面跑，还不到200M内存，真是累死机器不偿命啊，想不出错都难啊。
赞一下 回复此楼
一群自以为正义凛然的年轻人将一切不能以科学解释的事情定性为封建迷信并大刀阔斧地进行消灭，其实这是修养不足学识浅薄的一种体现，也是可恶的偏执和愚蠢的自以
6楼2010-04-04 23:31:31
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

aaq2800

铁杆木虫 (知名作家)
引用回帖:
Originally posted by zhangmt at 2010-04-04 23:31:31:
顺便多说一句，你的体系12个Al，外加1个Si，1个Sc，还用6-31g**基组，而且还是mp2计算，看来你得需要几台好机器了——好在是单点算，要是优化，估计你得算到毕业才能算完——不过你竟然还是在windows下面跑，还不 ...

没办法，人家审稿人要求算算，咱不敢违抗啊
赞一下 回复此楼
7楼2010-04-05 08:25:50
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fooo

金虫 (正式写手)
aaq2800(金币+1):谢谢吹捧 2010-04-05 14:51
换机器吧。
回复此楼
8楼2010-04-05 10:39:25
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

aaq2800

铁杆木虫 (知名作家)
引用回帖:
Originally posted by fooo at 2010-04-05 10:39:25:
换机器吧。

呵呵，还好，把H去掉，稍微降低一下精度，3个小时能跑完

[ Last edited by aaq2800 on 2010-4-5 at 14:58 ]
赞一下 回复此楼
9楼2010-04-05 14:51:06
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yjcmwgk

禁虫 (文坛精英)

密度泛函·小卒

优秀版主
aaq2800(金币+1):呵呵，你那么多bb啊 2010-04-05 14:59
引用回帖:
Originally posted by aaq2800 at 2010-04-05 08:25:50:

没办法，人家审稿人要求算算，咱不敢违抗啊

是审稿人要求你用波普尔型基组吗？如果是，你完全可以抗辩回去。这不是扯么
赞一下 回复此楼
10楼2010-04-05 14:55:18
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 aaq2800 的主题更新
201/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复(可上传附件)
信息提示
关闭
请填处理意见
关闭