24小时热门版块排行榜    

查看: 3351  |  回复: 9
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

physics7778

木虫 (正式写手)

初入文坛

[求助] 关于Molpro自旋轨道耦合输出的问题

做自旋轨道耦合计算时,总是出现下面的问题:
输入文件:
rvec=[0.8,0.9,1.0,1.1,1.2,1.3,1.4,1.5,1.6,1.7,1.8,1.9,2.0,2.1,2.2,2.3,2.4,2.5,2.6,2.7,2.8,2.9,3.0,3.1,3.2,3.3,3.4,3.5,3.6,3.7,3.8,3.9,4.0,4.1,4.2,4.3,4.4,4.5,4.6,4.7,4.8,4.9,5.0,5.1,5.2,5.3,5.4,5.5,5.6,5.7,5.8,5.9,6.0,7.0,8.0] ang
do i=1,#rvec
rGaH=rvec(i)
basis=avqz
{hf;
occ,9,3,3,1;
wf,32,1,0;
}
ehf(i)=energy(1);
{multi;
occ,11,4,4,1;
FROZEN,5,2,2,0;
closed,7,3,3,1;   
wf,32,1,2;state,1;
wf,32,2,2;state,1;
wf,32,3,2;state,1;}
{ci;core,5,2,2,0;occ,11,4,4,1;closed,7,3,3,1;wf,32,1,2;state,1;save,3012.1;noexc}
{ci;core,5,2,2,0;occ,11,4,4,1;closed,7,3,3,1;wf,32,2,2;state,1;save,3022.1;noexc}
{ci;core,5,2,2,0;occ,11,4,4,1;closed,7,3,3,1;wf,32,3,2;state,1;save,3032.1;noexc}
{ci;core,5,2,2,0;occ,11,4,4,1;closed,7,3,3,1;wf,32,1,2;save,4012.1}
{ci;core,5,2,2,0;occ,11,4,4,1;closed,7,3,3,1;wf,32,2,2;save,4022.1}
{ci;core,5,2,2,0;occ,11,4,4,1;closed,7,3,3,1;wf,32,3,2;save,4032.1}
LSINT
text,mrci,occ,11,4,4,1;core,5,2,2,0;closed,7,3,3,1;
{ci;hlsmat,ls,4012.1,4022.1,4032.1;print,hls=2,vls=0}
enddo
table,rvec,ehf
save,GaH.tab

出错信息
**********************************************************************************************************************************


*** mrci,occ,11,4,4,1;core,5,2,2,0;closed,7,3,3,1;
  

?ERROR IN GLOBAL_CREATE_INT FOR ARRAY   nblockss   LEN=       1

[ Last edited by physics7778 on 2011-5-4 at 14:35 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

分享,学习,提高
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

physics7778

木虫 (正式写手)

初入文坛

引用回帖:
Originally posted by 马永涛的虫 at 2011-05-10 20:27:36:
不好意思,最近有点忙,今天才看到。不知道你解决了没有。
感觉需要在你的LSINT命令后加上保存能量的命令。
比如你在下面进行了ci的计算,有没有想过保存ci计算的能量呢?
然后再最后的table命令里加上ci的 ...

谢谢老兄这么耐心的解答,我再按照你的办法试试。
分享,学习,提高
8楼2011-05-10 22:46:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 10 个回答

fatpig8832

铁杆木虫 (著名写手)

【答案】应助回帖

★
gmy1990(金币+1): 非常感谢你对论坛的贡献,奖励一下! 2011-05-04 18:31:53
physics7778(金币+1): 2011-05-04 18:45:56
text那一行后面不应该跟分号吧...

分号及mrci后面的逗号可以直接用空格代替
2楼2011-05-04 18:08:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

physics7778

木虫 (正式写手)

初入文坛

引用回帖:
Originally posted by fatpig8832 at 2011-05-04 18:08:05:
text那一行后面不应该跟分号吧...

分号及mrci后面的逗号可以直接用空格代替

这个应该不是问题的关键吧。
我想知道的是它最后提示的错误信息指的是什么。
分享,学习,提高
3楼2011-05-05 11:00:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fatpig8832

铁杆木虫 (著名写手)

引用回帖:
Originally posted by physics7778 at 2011-05-05 11:00:23:
这个应该不是问题的关键吧。
我想知道的是它最后提示的错误信息指的是什么。

那个错误信息看不懂,我估计就是指格式错误...
4楼2011-05-05 13:47:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 5sJFDKp2kwuI 2026-11-06 11/550 2026-11-09 04:16 by l48hFd8Nwae2
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 5sJFDKp2kwuI 2026-11-06 13/650 2026-11-09 04:09 by l48hFd8Nwae2
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 5sJFDKp2kwuI 2026-11-06 12/600 2026-11-09 04:05 by l48hFd8Nwae2
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 5sJFDKp2kwuI 2026-11-06 12/600 2026-11-09 03:58 by l48hFd8Nwae2
[公派出国] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 5sJFDKp2kwuI 2026-11-06 8/400 2026-11-09 00:14 by l48hFd8Nwae2
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 5sJFDKp2kwuI 2026-11-06 8/400 2026-11-08 23:01 by l48hFd8Nwae2
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +10 ncTVBM6AOJHF 2026-11-04 22/1100 2026-11-06 17:18 by 5sJFDKp2kwuI
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +6 ELBZ6WwEaRj6 2026-11-04 14/700 2026-11-06 16:32 by 5sJFDKp2kwuI
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +6 ttH3f1oA6ace 2026-11-04 12/600 2026-11-06 16:10 by 5sJFDKp2kwuI
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +7 ncTVBM6AOJHF 2026-11-04 17/850 2026-11-06 13:47 by BkA9Bi5cMhkD
[硕博家园] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +6 ncTVBM6AOJHF 2026-11-04 11/550 2026-11-06 13:47 by BkA9Bi5cMhkD
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +8 ncTVBM6AOJHF 2026-11-04 18/900 2026-11-06 13:36 by BkA9Bi5cMhkD
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 ELBZ6WwEaRj6 2026-11-04 10/500 2026-11-06 05:32 by qL1cU46YIdK0
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +6 kvewwGtKktqB 2026-11-04 16/800 2026-11-06 05:00 by qL1cU46YIdK0
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 ELBZ6WwEaRj6 2026-11-04 8/400 2026-11-06 02:24 by qL1cU46YIdK0
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 ELBZ6WwEaRj6 2026-11-04 9/450 2026-11-06 02:15 by qL1cU46YIdK0
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 ttH3f1oA6ace 2026-11-04 7/350 2026-11-06 02:13 by qL1cU46YIdK0
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 hNpp4Y2CGvSC 2026-11-04 12/600 2026-11-06 01:53 by qL1cU46YIdK0
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +6 ncTVBM6AOJHF 2026-11-04 10/500 2026-11-06 01:51 by qL1cU46YIdK0
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 ncTVBM6AOJHF 2026-11-04 12/600 2026-11-06 01:42 by qL1cU46YIdK0
信息提示
请填处理意见