physics7778

木蟲 (正式寫手)

初入文壇

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

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

做自旋軌道耦合計(jì)算時(shí)，總是出現(xiàn)下面的問(wè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

出錯(cuò)信息
**********************************************************************************************************************************

*** 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ù)此樓

» 猜你喜歡

0854電子信息272求調(diào)劑！已經(jīng)有0人回復(fù)
化學(xué)-化工-材料-食品-生物-環(huán)境2026年研究生調(diào)劑（可以不考數(shù)學(xué)）已經(jīng)有0人回復(fù)
物理化學(xué)論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有69人回復(fù)
化學(xué)-化工-材料-食品-生物-環(huán)境-農(nóng)業(yè) 2026年研究生調(diào)劑（可以不考數(shù)學(xué)）已經(jīng)有1人回復(fù)
070300一志愿211，312分求調(diào)劑院校已經(jīng)有7人回復(fù)
26年儲(chǔ)能、電池方向博士申請(qǐng) 已經(jīng)有0人回復(fù)
070300化學(xué)求調(diào)劑已經(jīng)有6人回復(fù)
青島大學(xué)化學(xué)化工學(xué)院2026年招收調(diào)劑碩士研究生已經(jīng)有0人回復(fù)
北方工業(yè)大學(xué)085808儲(chǔ)能技術(shù)-02方向（電池）接收調(diào)劑生已經(jīng)有0人回復(fù)
272求調(diào)劑，接受跨專業(yè)調(diào)劑！已經(jīng)有1人回復(fù)
有一區(qū)TOP論文，求調(diào)劑已經(jīng)有1人回復(fù)

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

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

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

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

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

分號(hào)及mrci后面的逗號(hào)可以直接用空格代替

贊一下(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小時(shí)
蟲號(hào): 763064
注冊(cè): 2009-05-04
性別: GG
專業(yè): 原子和分子物理

引用回帖:

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

分號(hào)及mrci后面的逗號(hào)可以直接用空格代替

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

贊一下

回復(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小時(shí)
蟲號(hào): 770403
注冊(cè): 2009-05-14

引用回帖:

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

那個(gè)錯(cuò)誤信息看不懂，我估計(jì)就是指格式錯(cuò)誤...

贊一下

回復(fù)此樓

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

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

馬永濤的蟲

新蟲 (小有名氣)

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

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

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

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

贊一下

回復(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小時(shí)
蟲號(hào): 763064
注冊(cè): 2009-05-04
性別: GG
專業(yè): 原子和分子物理

引用回帖:

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

減少數(shù)組中的數(shù)可以算下去，現(xiàn)在有個(gè)問(wèn)題就是如何讓自旋軌道耦合的計(jì)算結(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小時(shí)
蟲號(hào): 1097251
注冊(cè): 2010-09-13
性別: GG
專業(yè): 理論和計(jì)算化學(xué)

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

引用回帖:

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

不好意思，最近有點(diǎn)忙，今天才看到。不知道你解決了沒(méi)有。
感覺需要在你的LSINT命令后加上保存能量的命令。
比如你在下面進(jìn)行了ci的計(jì)算，有沒(méi)有想過(guò)保存ci計(jì)算的能量呢？
然后再最后的table命令里加上ci的能量，就可以保存你的自旋軌道計(jì)算結(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小時(shí)
蟲號(hào): 763064
注冊(cè): 2009-05-04
性別: GG
專業(yè): 原子和分子物理

引用回帖:

Originally posted by 馬永濤的蟲 at 2011-05-10 20:27:36:
不好意思，最近有點(diǎn)忙，今天才看到。不知道你解決了沒(méi)有。
感覺需要在你的LSINT命令后加上保存能量的命令。
比如你在下面進(jìn)行了ci的計(jì)算，有沒(méi)有想過(guò)保存ci計(jì)算的能量呢？
然后再最后的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小時(shí)
蟲號(hào): 695973
注冊(cè): 2009-02-04
性別: GG
專業(yè): 理論和計(jì)算化學(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ù)此樓

更痛苦的是夢(mèng)沒(méi)醒路已經(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小時(shí)
蟲號(hào): 763064
注冊(cè): 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ā)表

[考研] 285求調(diào)劑 +6	AZMK 2026-03-30	10/500	2026-03-31 11:41 by AZMK
[考研] 284求調(diào)劑 +5	小熊～～ 2026-03-31	5/250	2026-03-31 10:22 by sophie2180
[考研] 085600材料與化工調(diào)劑 +16	kikiki7 2026-03-30	16/800	2026-03-31 10:03 by 氯化亞硝酰
[考研] 0856調(diào)劑 +6	曲聽筠 2026-03-30	6/300	2026-03-31 09:44 by Jaylen.
[考研] 一志愿華東師范大學(xué)有機(jī)化學(xué)專業(yè)，初試351分，復(fù)試被刷求調(diào)劑! +7	真名有冰 2026-03-29	8/400	2026-03-30 21:08 by maaj02
[考研] 332求調(diào)劑 +14	032500 2026-03-25	14/700	2026-03-30 20:03 by 無(wú)際的草原
[考研] 367求調(diào)劑 +5	芋泥啵�！� 2026-03-28	5/250	2026-03-30 19:56 by 無(wú)際的草原
[考研] 0703本科鄭州大學(xué)求調(diào)劑 +7	nhj_ 2026-03-25	7/350	2026-03-30 12:44 by fangnagu
[考研] 296求調(diào)劑 +10	彼岸t 2026-03-29	10/500	2026-03-30 10:50 by 探123
[考研] 求調(diào)劑 +10	張zz111 2026-03-27	11/550	2026-03-30 09:17 by 無(wú)際的草原
[考研] 求調(diào)劑 +4	QiMing7 2026-03-25	5/250	2026-03-29 21:10 by 唐沐兒
[考研] 總分293求調(diào)劑 +8	加一一九 2026-03-25	11/550	2026-03-29 19:53 by 無(wú)際的草原
[考研] 數(shù)一英一271專碩（085401）求調(diào)劑，可跨 +7	前行必有光 2026-03-28	8/400	2026-03-28 23:22 by 小木蟲tim
[考研] 394求調(diào)劑 +3	好事多磨靜候佳?/a> 2026-03-26	5/250	2026-03-28 14:24 by 唐沐兒
[考研] 299求調(diào)劑 +7	嗯嗯嗯嗯2 2026-03-27	7/350	2026-03-28 13:09 by 唐沐兒
[考研] 材料277求調(diào)劑 +7	min3 2026-03-24	7/350	2026-03-28 11:39 by xuxiang
[考研] 339求調(diào)劑，想調(diào)回江蘇 +6	烤麥芽 2026-03-27	8/400	2026-03-28 10:40 by 烤麥芽
[考研] 266分求材料化工冶金礦業(yè)等專業(yè)的調(diào)劑 +4	哇呼哼呼哼 2026-03-26	4/200	2026-03-27 17:02 by zhyzzh
[考研] 求調(diào)劑一志愿本科北科大化學(xué) 343 +6	13831862839 2026-03-24	7/350	2026-03-26 22:57 by 不吃魚的貓
[考研] 化學(xué)調(diào)劑一志愿上海交通大學(xué)336分-本科上海211 +4	小魚愛有機(jī) 2026-03-25	4/200	2026-03-26 10:19 by aa331100