| 查看: 2781 | 回复: 19 | |||
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): 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 |
回复此楼» 猜你喜欢
伦敦大学学院Benjamin Moss博士招收人工智能、原位光谱、催化相关26年博士生
已经有11人回复
-
分子动力学模拟合作需求
已经有4人回复
-
物理化学论文润色/翻译怎么收费?
已经有242人回复
-
中国科学院成都山地灾害与环境研究所2026年特别研究助理（博士后）招聘启事
已经有0人回复
-
湖北文理学院材料与化工专硕还有调剂名额，欢迎大家推荐考生！
已经有1人回复
-
湖南大学刘巧玲课题组2026年博士研究生招生信息
已经有28人回复
-
南方科技大学表面催化反应动力学研究课题组招收2026届博士生
已经有1人回复
-
北京邮电大学物理学院老师招收博士生
已经有3人回复
-
求助：三乙酰丙酮铁（Fe(acac)₃）标准红外(IR)与XRD谱图/数据
已经有1人回复
-
Comsol 6.0 桶内肥料产热产湿过程模拟求助
已经有1人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- G09计算问题
已经有9人回复
- siesta计算弱相互作用VDW
已经有17人回复
- 【整理】今天整理硬盘上的资料，偶然发现的好玩的小东西就直接发上来吧
已经有106人回复
- 【求助】BSSE of O2运行出错
已经有17人回复
- 【求助】两个分子通过Cu2+配位的结合能问题
已经有7人回复
- 【求助】关于oniom计算的bsse问题
已经有5人回复
- 【求助】BSSE计算出错
已经有16人回复
- 【求助】用counter做bsse修正三个原子片段时出错
已经有7人回复
- 【求助】用BSSE校正计算分子片之间的相互作用能出现错误
已经有3人回复
wanfei
金虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 505
- 散金: 5
- 帖子: 105
- 在线: 11.3小时
- 虫号: 888862
- 注册: 2009-10-30
- 性别: GG
- 专业: 理论和计算化学
2楼2010-04-04 21:49:10
aaq2800
铁杆木虫 (知名作家)
- 应助: 51 (初中生)
- 金币: 28585.3
- 散金: 1522
- 红花: 29
- 沙发: 1
- 帖子: 9614
- 在线: 2190小时
- 虫号: 598401
- 注册: 2008-09-10
- 性别: GG
- 专业: 原子和分子物理
3楼2010-04-04 21:49:29
aaq2800
铁杆木虫 (知名作家)
- 应助: 51 (初中生)
- 金币: 28585.3
- 散金: 1522
- 红花: 29
- 沙发: 1
- 帖子: 9614
- 在线: 2190小时
- 虫号: 598401
- 注册: 2008-09-10
- 性别: GG
- 专业: 原子和分子物理
4楼2010-04-04 21:50:46
zhangmt
至尊木虫 (著名写手)
我叫MT
- QC强帖: 5
- 应助: 99 (初中生)
- 金币: 6961.8
- 散金: 10406
- 红花: 49
- 帖子: 1761
- 在线: 763.2小时
- 虫号: 880392
- 注册: 2009-10-22
- 性别: GG
- 专业: 理论和计算化学
★ ★
aaq2800(金币+2):谢谢，昨天晚上回去想了想，应该就是这里出错了， 2010-04-05 08:24
aylayl08(金币+2):谢谢专家指教 2010-04-05 20:47
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重态的，所以出错。 |
5楼2010-04-04 23:22:58
zhangmt
至尊木虫 (著名写手)
我叫MT
- QC强帖: 5
- 应助: 99 (初中生)
- 金币: 6961.8
- 散金: 10406
- 红花: 49
- 帖子: 1761
- 在线: 763.2小时
- 虫号: 880392
- 注册: 2009-10-22
- 性别: GG
- 专业: 理论和计算化学
6楼2010-04-04 23:31:31
aaq2800
铁杆木虫 (知名作家)
- 应助: 51 (初中生)
- 金币: 28585.3
- 散金: 1522
- 红花: 29
- 沙发: 1
- 帖子: 9614
- 在线: 2190小时
- 虫号: 598401
- 注册: 2008-09-10
- 性别: GG
- 专业: 原子和分子物理
7楼2010-04-05 08:25:50
8楼2010-04-05 10:39:25
aaq2800
铁杆木虫 (知名作家)
- 应助: 51 (初中生)
- 金币: 28585.3
- 散金: 1522
- 红花: 29
- 沙发: 1
- 帖子: 9614
- 在线: 2190小时
- 虫号: 598401
- 注册: 2008-09-10
- 性别: GG
- 专业: 原子和分子物理
9楼2010-04-05 14:51:06
yjcmwgk
禁虫 (文坛精英)
密度泛函·小卒
- QC强帖: 38
- 应助: 165 (高中生)
- 贵宾: 16.999
- 金币: 555.4
- 散金: 39880
- 红花: 395
- 沙发: 2
- 帖子: 10406
- 在线: 1528.1小时
- 虫号: 448368
- 注册: 2007-11-01
- 专业: 理论和计算化学
- 管辖: 科研工具资源
10楼2010-04-05 14:55:18
20