版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

返回列表

hdh912

金蟲 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 1180.9
帖子: 842
在線: 72.9小時(shí)
蟲號(hào): 387166

[交流] 【求助】請(qǐng)教molpro計(jì)算中出現(xiàn)的一個(gè)問(wèn)題

如題，我在對(duì)Na2分子的x1sigma+g, a3pi, A1sgma+u,B1PIU, C1PIU, D1PIU, E1PIU進(jìn)行掃描的時(shí)候，我的輸入卡如下：
  memory,30,m
geometry={na1;na2,na1,r1(i)}
distance= [2.1,2.3,2.5,2.6,2.8,2.9,2.95,3.0,3.05,3.1,3.15,3.2,3.42,3.62,3.9,4.2,4.5,4.8,5.3,5.7,5.8,5.9,6.0,6.2,6.4,6.6,6.8,7.0,7.2,7.4,7.6,7.8,8.0,8.2,8.4,8.6,8.8,9.0] ang
basis=VTZ
r1(1)=distance(1)
hf,occ,4,1,1,0,3,1,1,0;core,3,1,1,0,3,1,1,0;
i=0
do ir1=1,#distance
i=i+1
r1(i)=distance(ir1)
multi;
occ,4,1,1,0,3,1,1,0;core,3,1,1,0,3,1,1,0;wf,22,1,0;state,2,refstate,1,2;wf,22,5,0;state,1;wf,22,2,0,state,4,refstate,1,2,3,4;lquant,1;
ecas(i)=energy(1)
ci;occ,4,1,1,0,3,1,1,0;core,3,1,1,0,3,1,1,0;wf,22,1,0;state,2;wf,22,5,0;state,1;wf,22,2,0,state,4
emrci(i)=energy(1)
emrciq(i)=energd(1)
enddo
但是計(jì)算的結(jié)果始終出現(xiàn)如下信息：
State symmetry 3

  State symmetry 3:    Projection for operator LZ    squared value = 1

Number of electrons:  2 Spin symmetry=Singlet    Space symmetry=2
Number of states:    1
Number of CSFs:    0 (    0 determinants,    1 intermediate states)

Valence orbitals read from record    2100.2  Type=RHF/CANONICAL (state 1.1)
Frozen core orbitals read from record 2100.2  Type=RHF/CANONICAL (state 1.1)

EXTRA SYMMETRY OF MOS IN SYMMETRY 1: 0 0 0 1 1 1 1 2 1 1 2 2 1 1 1
EXTRA SYMMETRY OF MOS IN SYMMETRY 2: 0 1 1 1 1 1 2 1
EXTRA SYMMETRY OF MOS IN SYMMETRY 3: 0 1 1 1 1 1 2 1
EXTRA SYMMETRY OF MOS IN SYMMETRY 5: 0 0 0 1 1 1 1 2 1 1 2 1 2 1 1
EXTRA SYMMETRY OF MOS IN SYMMETRY 6: 0 1 1 1 1 2 1 1
EXTRA SYMMETRY OF MOS IN SYMMETRY 7: 0 1 1 1 1 2 1 1

Wavefunction dump at record          2140.2

Convergence thresholds  0.10E-01 (gradient)  0.10E-05 (energy)  0.10E-02 (step length)

Weight factors for state symmetry  1: 0.25000 0.25000
Weight factors for state symmetry  2: 0.25000
Weight factors for state symmetry  3: 0.25000

Number of orbital rotations    11    ( 0 Core/Active 0 Core/Virtual 0 Active/Active 11 Active/Virtual)
Total number of variables =    13

ITER. MIC  NCI  NEG    ENERGY(VAR)    ENERGY(PROJ) ENERGY CHANGE    GRAD(0)  GRAD(ORB) GRAD(CI)    STEP    TIME

MORE STATES THAN CONFIGURATIONS SPECIFIED:  1  2
懇請(qǐng)高手幫我修改下輸入卡，看怎么樣才能消除這個(gè)錯(cuò)誤，萬(wàn)分感謝。

回復(fù)此樓

» 猜你喜歡

高效微納米儲(chǔ)能材料與產(chǎn)業(yè)化推進(jìn)實(shí)驗(yàn)室蘇州招收碩士研究生（含調(diào)劑）已經(jīng)有5人回復(fù)
高效微納米儲(chǔ)能材料與產(chǎn)業(yè)化推進(jìn)實(shí)驗(yàn)室蘇州招收碩士研究生（含調(diào)劑）已經(jīng)有0人回復(fù)
物理化學(xué)論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有89人回復(fù)
高效微納米儲(chǔ)能材料與產(chǎn)業(yè)化推進(jìn)實(shí)驗(yàn)室蘇州招收碩士研究生（含調(diào)劑）已經(jīng)有2人回復(fù)
化學(xué)化工學(xué)院招收化工、化學(xué)、材料等相關(guān)方向研究生（學(xué)碩、專碩都有調(diào)劑名額）已經(jīng)有1人回復(fù)
高效微納米儲(chǔ)能材料與產(chǎn)業(yè)化推進(jìn)實(shí)驗(yàn)室蘇州招收碩士研究生（含調(diào)劑）已經(jīng)有0人回復(fù)
溫州大學(xué)新能源材料微觀溯源課題組招收2026級(jí)博士研究生（材料、化學(xué)方向）已經(jīng)有6人回復(fù)
材料與化工 323 英一+數(shù)二+物化，一志愿：哈工大本人本科雙一流已經(jīng)有5人回復(fù)
26申博已經(jīng)有2人回復(fù)
巴黎西岱大學(xué)（Université Paris Cité）Fayna MAMMERI 導(dǎo)師尋找CSC公派博士生已經(jīng)有9人回復(fù)

» 本主題相關(guān)商家推薦: (我也要在這里推廣)

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

Molpro和Gaussian算cas，我該信任那一個(gè)呢？已經(jīng)有7人回復(fù)
Molpro不能產(chǎn)生molden格式文件已經(jīng)有7人回復(fù)
購(gòu)買molpro的問(wèn)題已經(jīng)有4人回復(fù)
Molpro內(nèi)存不足已經(jīng)有22人回復(fù)
用MOLPRO構(gòu)造原子自然軌道（ANO）基函數(shù) 已經(jīng)有3人回復(fù)
求助一個(gè)算法問(wèn)題已經(jīng)有8人回復(fù)
gauss 計(jì)算電偶極矩已經(jīng)有22人回復(fù)
Molpro中MRCI級(jí)別參考系數(shù)的區(qū)別已經(jīng)有3人回復(fù)
molpro 運(yùn)行出錯(cuò) 已經(jīng)有3人回復(fù)
計(jì)算出一個(gè)催化反應(yīng)的活化能為10kj/mol，這個(gè)數(shù)值可靠嗎？已經(jīng)有12人回復(fù)
molpro說(shuō)明書的例子運(yùn)行不了？？已經(jīng)有4人回復(fù)
molpro計(jì)算自旋軌道耦合輸入求助已經(jīng)有3人回復(fù)
全國(guó)第二屆理論與計(jì)算化學(xué)軟件Molpro Molcas培訓(xùn)班已經(jīng)有19人回復(fù)
第一原理動(dòng)力學(xué)的系綜選擇，和各種系綜的輸入文件Example 已經(jīng)有13人回復(fù)
請(qǐng)教用molpro打印氫分子離子的重疊積矩陣、單電子積分矩陣、偶極矩陣已經(jīng)有3人回復(fù)
求助 molpro運(yùn)行出錯(cuò)！急已經(jīng)有5人回復(fù)
關(guān)于Molpro自旋軌道耦合輸出的問(wèn)題已經(jīng)有9人回復(fù)
【求助】molpro試算已經(jīng)有4人回復(fù)
【求助】molpro中使用彌散函數(shù)出錯(cuò) 已經(jīng)有4人回復(fù)
【求助】molpro 計(jì)算問(wèn)題，請(qǐng)高手指教！�。�！已經(jīng)有12人回復(fù)
【求助】molpro頻率輸出問(wèn)題已經(jīng)有4人回復(fù)

» 搶金幣啦！回帖就可以得到:

查看全部散金貼

1樓 2011-02-28 00:57:48

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

hdh912

金蟲 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 1180.9
帖子: 842
在線: 72.9小時(shí)
蟲號(hào): 387166

怎么沒(méi)有兄弟會(huì)嗎？

贊一下

回復(fù)此樓

2樓2011-03-01 09:46:07

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

fatpig8832

鐵桿木蟲 (著名寫手)

應(yīng)助: 16 (小學(xué)生)
金幣: 7420.2
帖子: 1438
在線: 610.9小時(shí)
蟲號(hào): 770403

hdh912(金幣+1): 謝謝，修改后還是一樣的錯(cuò)誤。 2011-03-02 00:26:55

具體不清楚，不過(guò)下面這一行

wf,22,2,0,state,4,refstate,1,2,3,4;lquant,1;

前面既然寫了"state,4"，后面的 lquant 是不是應(yīng)該相應(yīng)的跟著寫上四個(gè)數(shù)字，而不是只寫一個(gè)呢？

贊一下

回復(fù)此樓

3樓2011-03-01 11:12:51

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

beefly

專家顧問(wèn) (職業(yè)作家)

hdh912(金幣+3): 謝謝老大的回答，但是我將refstate刪掉后出錯(cuò)信息還是一樣。 2011-03-02 00:26:19

輸入輸出明顯不符。refstate不僅在casscf計(jì)算中沒(méi)有意義，而且把你的casscf輸入結(jié)構(gòu)破壞了。

贊一下

回復(fù)此樓

4樓2011-03-01 13:01:36

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

hdh912

金蟲 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 1180.9
帖子: 842
在線: 72.9小時(shí)
蟲號(hào): 387166

引用回帖:

Originally posted by beefly at 2011-03-01 13:01:36:
輸入輸出明顯不符。refstate不僅在casscf計(jì)算中沒(méi)有意義，而且把你的casscf輸入結(jié)構(gòu)破壞了。

老哥，你好我的全部輸入卡修改為
**,Na2
memory,30,m
geometry={na1;na2,na1,r1(i)}
distance= [2.1,2.3,2.5,2.6,2.8,2.9,2.95,3.0,3.05,3.1,3.15,3.2,3.42,3.62,3.9,4.2,4.5,4.8,5.3,5.7,5.8,5.9,6.0,6.2,6.4,6.6,6.8,7.0,7.2,7.4,7.6,7.8,8.0,8.2,8.4,8.6,8.8,9.0] ang
basis=VTZ
r1(1)=distance(1)
hf,occ,4,1,1,0,3,1,1,0;core,3,1,1,0,3,1,1,0;
i=0
do ir1=1,#distance
i=i+1
r1(i)=distance(ir1)
multi;
occ,4,1,1,0,3,1,1,0;core,3,1,1,0,3,1,1,0;wf,22,1,0;state,2;wf,22,5,0;state,1;wf,22,2,0,state,4;lquant,1,1,1,1;
ecas(i)=energy(1)
ci;occ,4,1,1,0,3,1,1,0;core,3,1,1,0,3,1,1,0;wf,22,1,0;state,2;wf,22,5,0;state,1;wf,22,2,0,state,4
emrci(i)=energy(1)
emrciq(i)=energd(1)
enddo
table,r1,ecas,emrci,emrciq
head,r1,ecas,emrci,emrciq
save,na2.tab
輸出信息還是
State symmetry 3

  State symmetry 3:    Projection for operator LZ    squared value = 1

Number of electrons:  2 Spin symmetry=Singlet    Space symmetry=2
Number of states:    1
Number of CSFs:    0 (    0 determinants,    1 intermediate states)

Valence orbitals read from record    2100.2  Type=RHF/CANONICAL (state 1.1)
Frozen core orbitals read from record 2100.2  Type=RHF/CANONICAL (state 1.1)

EXTRA SYMMETRY OF MOS IN SYMMETRY 1: 0 0 0 1 1 1 1 2 1 1 2 2 1 1 1
EXTRA SYMMETRY OF MOS IN SYMMETRY 2: 0 1 1 1 1 1 2 1
EXTRA SYMMETRY OF MOS IN SYMMETRY 3: 0 1 1 1 1 1 2 1
EXTRA SYMMETRY OF MOS IN SYMMETRY 5: 0 0 0 1 1 1 1 2 1 1 2 1 2 1 1
EXTRA SYMMETRY OF MOS IN SYMMETRY 6: 0 1 1 1 1 2 1 1
EXTRA SYMMETRY OF MOS IN SYMMETRY 7: 0 1 1 1 1 2 1 1

Wavefunction dump at record          2140.2

Convergence thresholds  0.10E-01 (gradient)  0.10E-05 (energy)  0.10E-02 (step length)

Weight factors for state symmetry  1: 0.25000 0.25000
Weight factors for state symmetry  2: 0.25000
Weight factors for state symmetry  3: 0.25000

Number of orbital rotations    11    ( 0 Core/Active 0 Core/Virtual 0 Active/Active 11 Active/Virtual)
Total number of variables =    13

ITER. MIC  NCI  NEG    ENERGY(VAR)    ENERGY(PROJ) ENERGY CHANGE    GRAD(0)  GRAD(ORB) GRAD(CI)    STEP    TIME

MORE STATES THAN CONFIGURATIONS SPECIFIED:  1  2
懇請(qǐng)老哥幫忙。謝謝

贊一下

回復(fù)此樓

5樓2011-03-02 00:29:23

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

fatpig8832

鐵桿木蟲 (著名寫手)

應(yīng)助: 16 (小學(xué)生)
金幣: 7420.2
帖子: 1438
在線: 610.9小時(shí)
蟲號(hào): 770403

hdh912(金幣+1): 謝謝兄弟不過(guò)刪掉后還是不行 2011-03-02 23:08:12

從這一段來(lái)看
Weight factors for state symmetry  1: 0.25000 0.25000
Weight factors for state symmetry  2: 0.25000
Weight factors for state symmetry  3: 0.25000

symmetry 3 只算了一個(gè)態(tài)，與你之前設(shè)定的
wf,22,2,0,state,4；
不符

不如刪去 lquant,1,1,1,1; 試試...

贊一下

回復(fù)此樓

6樓2011-03-02 09:49:32

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

beefly

專家顧問(wèn) (職業(yè)作家)

hdh912(金幣+5): 謝謝老哥，老哥能否說(shuō)仔細(xì)點(diǎn) 輸入卡應(yīng)該怎么改，真的麻煩了 2011-03-02 23:09:22

活性空間錯(cuò)誤：2個(gè)活性電子，一個(gè)活性軌道，實(shí)際上是一個(gè)RHF計(jì)算，只能有一個(gè)1A1態(tài)。

贊一下(1人)

回復(fù)此樓

7樓2011-03-02 10:56:44

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

知222

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
帖子: 12
在線: 9.9小時(shí)
蟲號(hào): 7783046

★
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖

引用回帖:

2樓: Originally posted by hdh912 at 2011-03-01 09:46:07
怎么沒(méi)有兄弟會(huì)嗎？

把hf改為{hf;wf,22,1,0;}再試試

贊一下

回復(fù)此樓

8樓2018-12-14 11:44:25

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

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 hdh912 的主題更新

返回列表

普通表情龍兔虎貓高級(jí)回復(fù) (可上傳附件)

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

[考研] 一志愿北京化工大學(xué) 070300 學(xué)碩 336分求調(diào)劑 +4	vv迷 2026-03-22	4/200	2026-03-22 23:29 by king123！
[考研] 317求調(diào)劑 +12	申子申申 2026-03-19	18/900	2026-03-22 22:23 by luoyongfeng
[考研] 0854電子信息求調(diào)劑 +3	α____ 2026-03-22	3/150	2026-03-22 21:28 by zhq0425
[考研] 一志愿中南化學(xué)（0703）總分337求調(diào)劑 +9	niko- 2026-03-19	10/500	2026-03-22 16:08 by ColorlessPI
[考研] 生物學(xué)調(diào)劑 +5	Surekei 2026-03-21	5/250	2026-03-22 14:39 by tcx007
[考研] 285求調(diào)劑 +6	ytter 2026-03-22	6/300	2026-03-22 12:09 by 星空星月
[考研] 初試 317 +7	半拉月丙 2026-03-20	7/350	2026-03-21 22:26 by peike
[考研] 工科0856求調(diào)劑 +3	沐析汀汀 2026-03-21	3/150	2026-03-21 18:30 by 學(xué)員8dgXkO
[考研] 265求調(diào)劑 +12	梁梁校校 2026-03-19	14/700	2026-03-21 13:38 by lature00
[考研] 085601調(diào)劑 358分 +3	zzzzggh 2026-03-20	4/200	2026-03-21 10:21 by luoyongfeng
[考研] 二本跨考鄭大材料306英一數(shù)二 +3	z1z2z3879 2026-03-17	3/150	2026-03-21 02:29 by JourneyLucky
[考研] 求調(diào)劑，一志愿:南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩，總分289分 +4	@taotao 2026-03-19	4/200	2026-03-20 22:14 by JourneyLucky
[考研] 工科材料085601 279求調(diào)劑 +7	困于星晨 2026-03-17	9/450	2026-03-20 17:38 by 無(wú)懈可擊111
[考研] 生物學(xué)調(diào)劑招人�。。� +3	山海天嵐 2026-03-17	4/200	2026-03-19 21:34 by 怎么釋懷
[考研] 085601材料工程專碩求調(diào)劑 +10	慕寒mio 2026-03-16	10/500	2026-03-19 15:26 by 丁丁*
[考研] 312求調(diào)劑 +8	陌宸希 2026-03-16	9/450	2026-03-18 12:39 by Linda Hu
[考博] 26博士申請(qǐng) +3	1042136743 2026-03-17	3/150	2026-03-17 23:30 by 輕松不少隨
[考研] 085601求調(diào)劑 +4	Du.11 2026-03-16	4/200	2026-03-17 17:08 by ruiyingmiao
[考研] 283求調(diào)劑 +3	聽(tīng)風(fēng)就是雨； 2026-03-16	3/150	2026-03-17 07:41 by 熱情沙漠
[考研] 333求調(diào)劑 +3	文思客 2026-03-16	7/350	2026-03-16 18:21 by 文思客

<td id="kfhah"></td>

24小時(shí)熱門版塊排行榜

[交流] 【求助】請(qǐng)教molpro計(jì)算中出現(xiàn)的一個(gè)問(wèn)題

» 猜你喜歡

» 本主題相關(guān)商家推薦: (我也要在這里推廣)

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

» 搶金幣啦！回帖就可以得到: 查看全部散金貼

» 搶金幣啦！回帖就可以得到:

查看全部散金貼