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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Molpro/Molcas » 关于Molpro自旋轨道耦合输出的问题
101/1返回列表
查看: 3207  |  回复: 9
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
分享，学习，提高
1楼 2011-05-04 14:31:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fatpig8832

铁杆木虫 (著名写手)

【答案】应助回帖

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

分号及mrci后面的逗号可以直接用空格代替
赞一下(1人) 回复此楼
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的回帖

马永涛的虫

新虫 (小有名气)

【答案】应助回帖

physics7778(金币+4): 谢谢啦，我再试试看. 2011-05-05 19:24:29
好像是do循环的问题。数组中数太多了，试一下分两次计算。
赞一下 回复此楼
5楼2011-05-05 18:50:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

physics7778

木虫 (正式写手)

初入文坛
引用回帖:
Originally posted by 马永涛的虫 at 2011-05-05 18:50:39:
好像是do循环的问题。数组中数太多了，试一下分两次计算。

减少数组中的数可以算下去，现在有个问题就是如何让自旋轨道耦合的计算结果打印成表格的形式，以便绘图。
赞一下 回复此楼
分享，学习，提高
6楼2011-05-06 08:33:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

马永涛的虫

新虫 (小有名气)

【答案】应助回帖

引用回帖:
Originally posted by physics7778 at 2011-05-06 08:33:17:
减少数组中的数可以算下去，现在有个问题就是如何让自旋轨道耦合的计算结果打印成表格的形式，以便绘图。

不好意思，最近有点忙，今天才看到。不知道你解决了没有。
感觉需要在你的LSINT命令后加上保存能量的命令。
比如你在下面进行了ci的计算，有没有想过保存ci计算的能量呢？
然后再最后的table命令里加上ci的能量，就可以保存你的自旋轨道计算结果。
不知这是不是你要的？
呵呵。
赞一下 回复此楼
7楼2011-05-10 20:27:36
已阅   回复此楼   关注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的回帖

mchen10

金虫 (正式写手)

【答案】应助回帖

用
r0=0.6
dr=0.2
do i=1,40
   r(i) = r0 + dr
   [your code here]
   e(i)=energy
enddo

then do the e(i)-r(i) table
赞一下 回复此楼
更痛苦的是梦没醒路已经不能走了
9楼2011-05-10 23:43:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

physics7778

木虫 (正式写手)

初入文坛
引用回帖:
Originally posted by mchen10 at 2011-05-10 23:43:45:
用
r0=0.6
dr=0.2
do i=1,40
   r(i) = r0 + dr
   [your code here]
   e(i)=energy
enddo

then do the e(i)-r(i) table

谢谢大家的应助和提醒！
赞一下 回复此楼
分享，学习，提高
10楼2011-05-12 08:39:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 physics7778 的主题更新
101/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +37 医学老男孩 2026-08-13 86/4300 2026-08-18 00:18 by longsantaizi
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[基金申请] 时间戳变了，能看出什么问题？ +4 基诺咪客 2026-08-17 4/200 2026-08-17 16:10 by Vivilian
[基金申请] 今天系统多次维护，明天很可能放榜！ +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[基金申请] 时间戳又变了8-15 +13 archvillain 2026-08-15 25/1250 2026-08-16 20:13 by zhaosm1982
[基金申请] 2027广东省杰青 +3 奶牛小黑 2026-08-15 6/300 2026-08-16 20:07 by 奶牛小黑
[基金申请] 快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。 +3 Tide man 2026-08-14 3/150 2026-08-16 17:47 by jurkat.1640
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 26/1300 2026-08-16 12:35 by 启萌科技
[基金申请] 有时候，自然基金真的不能太认真 （我的申报经验） +10 majunge000 2026-08-11 12/600 2026-08-16 08:18 by xli1984
[精细化工] 招聘 金属平磨液，抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 各位道友，我要去昆明玩几天，回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 是这周出结果还是下周出结果？ +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[硕博家园] 读博的好处 +4 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[基金申请] 重要来源：本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] Filecode 又变了，巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 2019年青年基金涵评意见，大家看看几个A，几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 什么时候出结果，有咨询渠道？？？ +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
信息提示
关闭
请填处理意见
关闭