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

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 【求助】BSSE计算出错
52/1返回列表
查看: 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的回帖

yjcmwgk

禁虫 (文坛精英)

密度泛函·小卒

优秀版主
aaq2800(金币+1):bb快分完了 2010-04-05 19:29
引用回帖:
Originally posted by aaq2800 at 2010-04-05 14:58:22:

这种弱的范德瓦尔斯相互作用，DFT不能处理，MP2算的结果还比较可靠些，所以得用来验证些结果的可靠性。

晕！
我质疑的是你用的基组
又没质疑你用的方法……汗
赞一下 回复此楼
12楼2010-04-05 15:17:06
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 20 个回答

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的回帖
查看全部 20 个回答
普通表情 龙 兔 虎 猫 高级回复(可上传附件)
信息提示
关闭
请填处理意见
关闭