24小时热门版块排行榜    

查看: 348  |  回复: 3
当前主题已经存档。

su571

[交流] [求助]Rh-.多重度1计算能量错误

计算Rh-.多重度1计算能量
输入文件是:
%chk=Rh.chk
% rwf=1, 200mw, 2, 200mw, 3, 200mw, 4, 100mw, 5, 100mw, 6, 100mw, 7, 100mw
%nosave
# b3lyp/lanl2dz  scf(maxcyc=300) int=ultrafine
Rh
-1 1
Rh

--Link1--
%chk=Rh.chk
% rwf=1, 200mw, 2, 200mw, 3, 200mw, 4, 100mw, 5, 100mw, 6, 100mw, 7, 100mw
%nosave
# bpw91/genecp  scf(maxcyc=500,tight)  Guess=Read geom=allcheck
int=ultrafine

Rh 0
S 9 1.00
0.29293500E+02 0.66660000E-02
0.16682700E+02 -0.49734000E-01
0.10423500E+02 0.13155900E+00
0.30952600E+01 -0.52438200E+00
0.70797100E+00 0.78637900E+00
0.32256700E+00 0.46395800E+00
0.69085000E-01 0.20486000E-01
0.27415000E-01 -0.33670000E-02
0.29840000E+00 0.00000000E+00
S 9 1.00
0.29293500E+02 -0.13580000E-02
0.16682700E+02 0.12392000E-01
0.10423500E+02 -0.37572000E-01
0.30952600E+01 0.17410900E+00
0.70797100E+00 -0.33850900E+00
0.32256700E+00 -0.33541700E+00
0.69085000E-01 0.65814700E+00
0.27415000E-01 0.54051000E+00
0.29840000E+00 0.00000000E+00
S 9 1.00
0.29293500E+02 0.11242000E-01
0.16682700E+02 -0.57401000E-01
0.10423500E+02 0.11064800E+00
0.30952600E+01 -0.37694300E+00
0.70797100E+00 0.12026360E+01
0.32256700E+00 -0.48758200E+00
0.69085000E-01 -0.18636390E+01
0.27415000E-01 0.17964270E+01
0.29840000E+00 0.00000000E+00
S 1 1.00
0.27415000E-01 0.10000000E+01
S 1 1.00
0.29840000E+00 0.10000000E+01
P 8 1.00
0.90271500E+01 0.20396000E-01
0.40242700E+01 -0.15366200E+00
0.99067400E+00 0.43208200E+00
0.47279700E+00 0.49834800E+00
0.21973400E+00 0.19574200E+00
0.80088000E-01 0.14083000E-01
0.30640000E-01 0.60000000E-04
0.33270000E+00 0.00000000E+00
P 8 1.00
0.90271500E+01 -0.56190000E-02
0.40242700E+01 0.46511000E-01
0.99067400E+00 -0.14952700E+00
0.47279700E+00 -0.22433700E+00
0.21973400E+00 -0.11470000E-01
0.80088000E-01 0.63649700E+00
0.30640000E-01 0.47462500E+00
0.33270000E+00 0.00000000E+00
P 8 1.00
0.90271500E+01 0.11131000E-01
0.40242700E+01 -0.96301000E-01
0.99067400E+00 0.29868500E+00
0.47279700E+00 0.88921900E+00
0.21973400E+00 -0.14854780E+01
0.80088000E-01 -0.11311200E+00
0.30640000E-01 0.89298500E+00
0.33270000E+00 0.00000000E+00
P 1 1.00
0.30640000E-01 0.10000000E+01
P 1 1.00
0.33270000E+00 0.10000000E+01
D 6 1.00
0.43045000E+01 -0.43910000E-02
0.13201900E+01 0.11968200E+00
0.56140500E+00 0.29627400E+00
0.22670200E+00 0.39087300E+00
0.87949000E-01 0.35396100E+00
0.32305000E-01 0.15016200E+00
D 6 1.00
0.43045000E+01 0.44330000E-02
0.13201900E+01 -0.14991000E+00
0.56140500E+00 -0.37329300E+00
0.22670200E+00 -0.31298700E+00
0.87949000E-01 0.42189700E+00
0.32305000E-01 0.59366000E+00
D 1 1.00
0.32305000E-01 0.10000000E+01
D 1 1.00
0.98300000E+00 0.10000000E+01
F 1 1.00
0.84990000E+00 0.10000000E+01
F 1 1.00
0.20280000E+00 0.10000000E+01
****

Rh 0
ECP28MDF,4,28
G-Komponente
1
2 1.000000 0.000000
S-G
2
2 7.858275 135.134974
2 3.382128 15.411632
P-G
4
2 6.849791 29.251437
2 6.710092 58.508363
2 3.042159 3.780243
2 2.937330 7.676547
D-G
4
2 5.416315 11.849911
2 5.333416 17.778103
2 1.976212 2.062383
2 1.961111 3.075654
F-G
2
2 5.028590 -6.928078
2 5.005582 -9.155099
输出文件是
dumping /fiocom/, unit = 2 NFiles =     1 SizExt =         0 WInBlk =       512
                  defal = F LstWrd =       66048 FType=2 FMxFil=10000

Number           0
Base         20480
End          66048
End1         66048
Wr Pntr      20480
Rd Pntr      20480
Length       45568

dumping /fiocom/, unit = 3 NFiles =     1 SizExt =    524288 WInBlk =       512
                  defal = T LstWrd =       66048 FType=2 FMxFil=10000

Number           0
Base         20480
End          66048
End1         66048
Wr Pntr      20480
Rd Pntr      20480
Length       45568
Error termination in NtrErr:
NtrErr Called from FileIO.
搜过以前的帖子,没看出解决门道,哪位能指点下,实在是非常感谢
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cuihang

铁杆木虫 (正式写手)

VASPer


lei0736(金币+1,VIP+0):谢谢
电荷数不是0,那么多重读就不该是1了
聪明人知道自己是傻子,傻子以为自己是聪明人
2楼2008-03-13 11:02:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Roseshu

金虫 (正式写手)

★ ★
lei0736(金币+1,VIP+0):多谢
zzgyb(金币+1,VIP+0):再奖励一个
"%chk=Rh.chk
% rwf=1, 200mw, 2, 200mw, 3, 200mw, 4, 100mw, 5, 100mw, 6, 100mw, 7, 100mw
%nosave
# b3lyp/lanl2dz  scf(maxcyc=300) int=ultrafine"

你第一个计算的chk文件都nosave了,后面的link当然读不到check文件了。
3楼2008-03-13 12:46:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

su571

谢谢Roseshu,第一步是原来编的,没注意,改了以后有结果了,十分感谢,最近一直在这些东西上困扰啊.
同时也谢谢cuihang,Rh带负一价电荷,并不是未成对电子为一
4楼2008-03-13 15:34:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 su571 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见