physics7778

木蟲 (正式寫手)

初入文壇

應(yīng)助: 4 (幼兒園)
金幣: 2178.1
紅花: 2
帖子: 507
在線: 465.5小時
蟲號: 763064
注冊: 2009-05-04
性別: GG
專業(yè): 原子和分子物理

[求助] 關(guān)于Molpro自旋軌道耦合輸出的問題

做自旋軌道耦合計算時，總是出現(xiàn)下面的問題：
輸入文件：
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 ]

回復(fù)此樓

» 猜你喜歡

2026湖北師范大學(xué)化學(xué)化工學(xué)院招收研究生調(diào)劑已經(jīng)有0人回復(fù)
26年材料/化學(xué)相關(guān)專業(yè)博士研究生招生溫州大學(xué)化學(xué)與材料工程學(xué)院(26年秋季入學(xué)) 已經(jīng)有0人回復(fù)
物理化學(xué)論文潤色/翻譯怎么收費? 已經(jīng)有294人回復(fù)
0703調(diào)劑已經(jīng)有2人回復(fù)
2026年博士研究生招生已經(jīng)有1人回復(fù)
304求調(diào)劑已經(jīng)有6人回復(fù)
關(guān)于研究生調(diào)劑的些許建議已經(jīng)有1人回復(fù)
塑料油分析測試已經(jīng)有1人回復(fù)

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

molpro計算自旋軌道耦合輸入求助已經(jīng)有3人回復(fù)
石墨烯中的自旋軌道耦合已經(jīng)有7人回復(fù)
全國第二屆理論與計算化學(xué)軟件Molpro Molcas培訓(xùn)班已經(jīng)有19人回復(fù)
ADF計算自旋軌道耦合已經(jīng)有7人回復(fù)
關(guān)于vasp計算自旋軌道耦合的問題已經(jīng)有4人回復(fù)
求助 molpro運行出錯！急已經(jīng)有5人回復(fù)
【討論】molpro, molcas和gaussian的優(yōu)缺點已經(jīng)有3人回復(fù)
【整理】今天整理硬盤上的資料，偶然發(fā)現(xiàn)的好玩的小東西就直接發(fā)上來吧已經(jīng)有105人回復(fù)
【求助】請教molpro計算中出現(xiàn)的一個問題已經(jīng)有7人回復(fù)
【問題已解決】vasp計算自旋軌道耦合讀問題已經(jīng)有10人回復(fù)
【求助】VASP考慮自旋軌道耦合的話，如何編譯已經(jīng)有17人回復(fù)
【求助】請問，把molpro和molcas 鉆研透徹之后，能有什么用處�。� 已經(jīng)有9人回復(fù)
【求助】Molpro如何執(zhí)行旋軌耦合計算如何執(zhí)行態(tài)平均計算? 已經(jīng)有10人回復(fù)
【求助】td計算三重態(tài)問題已經(jīng)有8人回復(fù)

分享，學(xué)習(xí)，提高

1樓 2011-05-04 14:31:35

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

fatpig8832

鐵桿木蟲 (著名寫手)

應(yīng)助: 16 (小學(xué)生)
金幣: 7420.2
散金: 31
紅花: 28
帖子: 1438
在線: 610.9小時
蟲號: 770403
注冊: 2009-05-14

【答案】應(yīng)助回帖

★
gmy1990(金幣+1): 非常感謝你對論壇的貢獻，獎勵一下！ 2011-05-04 18:31:53
physics7778(金幣+1): 2011-05-04 18:45:56

text那一行后面不應(yīng)該跟分號吧...

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

贊一下(1人)

回復(fù)此樓

2樓2011-05-04 18:08:05

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

physics7778

木蟲 (正式寫手)

初入文壇

應(yīng)助: 4 (幼兒園)
金幣: 2178.1
紅花: 2
帖子: 507
在線: 465.5小時
蟲號: 763064
注冊: 2009-05-04
性別: GG
專業(yè): 原子和分子物理

引用回帖:

Originally posted by fatpig8832 at 2011-05-04 18:08:05:
text那一行后面不應(yīng)該跟分號吧...

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

這個應(yīng)該不是問題的關(guān)鍵吧。
我想知道的是它最后提示的錯誤信息指的是什么。

贊一下

回復(fù)此樓

分享，學(xué)習(xí)，提高

3樓2011-05-05 11:00:23

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

fatpig8832

鐵桿木蟲 (著名寫手)

應(yīng)助: 16 (小學(xué)生)
金幣: 7420.2
散金: 31
紅花: 28
帖子: 1438
在線: 610.9小時
蟲號: 770403
注冊: 2009-05-14

引用回帖:

Originally posted by physics7778 at 2011-05-05 11:00:23:
這個應(yīng)該不是問題的關(guān)鍵吧。
我想知道的是它最后提示的錯誤信息指的是什么。

那個錯誤信息看不懂，我估計就是指格式錯誤...

贊一下

回復(fù)此樓

4樓2011-05-05 13:47:42

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

馬永濤的蟲

新蟲 (小有名氣)

應(yīng)助: 3 (幼兒園)
金幣: 33.2
帖子: 70
在線: 20.8小時
蟲號: 1097251
注冊: 2010-09-13
性別: GG
專業(yè): 理論和計算化學(xué)

【答案】應(yīng)助回帖

physics7778(金幣+4): 謝謝啦，我再試試看. 2011-05-05 19:24:29

好像是do循環(huán)的問題。數(shù)組中數(shù)太多了，試一下分兩次計算。

贊一下

回復(fù)此樓

5樓2011-05-05 18:50:39

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

physics7778

木蟲 (正式寫手)

初入文壇

應(yīng)助: 4 (幼兒園)
金幣: 2178.1
紅花: 2
帖子: 507
在線: 465.5小時
蟲號: 763064
注冊: 2009-05-04
性別: GG
專業(yè): 原子和分子物理

引用回帖:

Originally posted by 馬永濤的蟲 at 2011-05-05 18:50:39:
好像是do循環(huán)的問題。數(shù)組中數(shù)太多了，試一下分兩次計算。

減少數(shù)組中的數(shù)可以算下去，現(xiàn)在有個問題就是如何讓自旋軌道耦合的計算結(jié)果打印成表格的形式，以便繪圖。

贊一下

回復(fù)此樓

分享，學(xué)習(xí)，提高

6樓2011-05-06 08:33:17

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

馬永濤的蟲

新蟲 (小有名氣)

應(yīng)助: 3 (幼兒園)
金幣: 33.2
帖子: 70
在線: 20.8小時
蟲號: 1097251
注冊: 2010-09-13
性別: GG
專業(yè): 理論和計算化學(xué)

【答案】應(yīng)助回帖

引用回帖:

Originally posted by physics7778 at 2011-05-06 08:33:17:
減少數(shù)組中的數(shù)可以算下去，現(xiàn)在有個問題就是如何讓自旋軌道耦合的計算結(jié)果打印成表格的形式，以便繪圖。

不好意思，最近有點忙，今天才看到。不知道你解決了沒有。
感覺需要在你的LSINT命令后加上保存能量的命令。
比如你在下面進行了ci的計算，有沒有想過保存ci計算的能量呢？
然后再最后的table命令里加上ci的能量，就可以保存你的自旋軌道計算結(jié)果。
不知這是不是你要的？
呵呵。

贊一下

回復(fù)此樓

7樓2011-05-10 20:27:36

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

physics7778

木蟲 (正式寫手)

初入文壇

應(yīng)助: 4 (幼兒園)
金幣: 2178.1
紅花: 2
帖子: 507
在線: 465.5小時
蟲號: 763064
注冊: 2009-05-04
性別: GG
專業(yè): 原子和分子物理

引用回帖:

Originally posted by 馬永濤的蟲 at 2011-05-10 20:27:36:
不好意思，最近有點忙，今天才看到。不知道你解決了沒有。
感覺需要在你的LSINT命令后加上保存能量的命令。
比如你在下面進行了ci的計算，有沒有想過保存ci計算的能量呢？
然后再最后的table命令里加上ci的 ...

謝謝老兄這么耐心的解答，我再按照你的辦法試試。

贊一下

回復(fù)此樓

分享，學(xué)習(xí)，提高

8樓2011-05-10 22:46:13

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

mchen10

金蟲 (正式寫手)

應(yīng)助: 13 (小學(xué)生)
金幣: 1508.2
散金: 41
紅花: 6
帖子: 425
在線: 159.5小時
蟲號: 695973
注冊: 2009-02-04
性別: GG
專業(yè): 理論和計算化學(xué)

【答案】應(yīng)助回帖

用
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

贊一下

回復(fù)此樓

更痛苦的是夢沒醒路已經(jīng)不能走了

9樓2011-05-10 23:43:45

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

physics7778

木蟲 (正式寫手)

初入文壇

應(yīng)助: 4 (幼兒園)
金幣: 2178.1
紅花: 2
帖子: 507
在線: 465.5小時
蟲號: 763064
注冊: 2009-05-04
性別: GG
專業(yè): 原子和分子物理

引用回帖:

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

謝謝大家的應(yīng)助和提醒！

贊一下

回復(fù)此樓

分享，學(xué)習(xí)，提高

10樓2011-05-12 08:39:27

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 一志愿a區(qū)211，085601-307分求調(diào)劑 +7	黨嘉豪 2026-03-31	15/750	2026-03-31 14:19 by JourneyLucky
[考研] 08工科，295，接受跨專業(yè)調(diào)劑 +3	lmnlzy 2026-03-31	3/150	2026-03-31 12:36 by 探123
[考研] 求調(diào)劑 +4	圖鑒212 2026-03-30	4/200	2026-03-31 10:20 by cal0306
[考研] 南京大學(xué)化學(xué)調(diào)劑 +11	景隨風(fēng) 2026-03-29	16/800	2026-03-31 10:14 by herarysara
[考研] 085600，專業(yè)課化工原理，320分求調(diào)劑 +6	大饞小子 2026-03-29	6/300	2026-03-31 10:03 by 氯化亞硝酰
[有機交流] 甲基亞磺磺酸鈉和甲基磺酸酯反應(yīng)機理 10+3	kaobao456 2026-03-29	4/200	2026-03-30 23:16 by nBu鋰
[考研] 359求調(diào)劑 +5	王了個楠 2026-03-25	5/250	2026-03-30 19:36 by 源_2020
[考研] 334分一志愿武理材料求調(diào)劑 +16	李李不服輸 2026-03-26	16/800	2026-03-30 17:47 by wangjy2002
[考研] 考研調(diào)劑 +5	Sanmu-124 2026-03-26	5/250	2026-03-30 13:36 by chemdavid
[考研] 348求調(diào)劑 +6	小懶蟲不懶了 2026-03-28	6/300	2026-03-30 10:29 by Evan_Liu
[考研] 356求調(diào)劑 +4	gysy?s?a 2026-03-28	4/200	2026-03-29 10:32 by 唐沐兒
[考研] 305求調(diào)劑 +8	RuiFairyrui 2026-03-28	8/400	2026-03-29 08:22 by fmesaito
[考研] 求調(diào)劑 +7	爭取九點睡 2026-03-28	8/400	2026-03-28 21:07 by 爭取九點睡
[考研] 材料277求調(diào)劑 +7	min3 2026-03-24	7/350	2026-03-28 11:39 by xuxiang
[考研] 0856，材料與化工321分求調(diào)劑 +12	大饞小子 2026-03-27	13/650	2026-03-28 10:56 by self2008
[考研] 308求調(diào)劑 +7	墨墨漠 2026-03-27	7/350	2026-03-28 07:43 by 熱情沙漠
[考研] 086502化學(xué)工程342求調(diào)劑 +6	阿姨復(fù)古不過 2026-03-27	6/300	2026-03-28 07:06 by wangy0907
[考研] 266分求材料化工冶金礦業(yè)等專業(yè)的調(diào)劑 +4	哇呼哼呼哼 2026-03-26	4/200	2026-03-27 17:02 by zhyzzh
[考研] 復(fù)試調(diào)劑，一志愿南農(nóng)083200食品科學(xué)與工程 +5	XQTJZ 2026-03-26	5/250	2026-03-27 14:49 by 狂炫麥當(dāng)當(dāng)
[考研] 一志愿鄭大085600，310分求調(diào)劑 +5	李瀟可 2026-03-26	5/250	2026-03-27 11:14 by 不吃魚的貓