版塊導(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í)熱門版塊排行榜

返回列表

當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖

kdswx

新蟲(chóng) (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1799
帖子: 167
在線: 47小時(shí)
蟲(chóng)號(hào): 2736957
注冊(cè): 2013-10-19

[求助] 求大神們幫忙看看這個(gè)程序哪兒需要改啊已有1人參與

初學(xué)molpro，學(xué)著計(jì)算NO基態(tài)的勢(shì)能，想要單點(diǎn)能掃描，這下面是計(jì)算程序，但是出來(lái)的結(jié)果不對(duì)啊，能不能請(qǐng)大神們幫忙看看，萬(wàn)分感謝！
***,NO potential

geometry={h
      n;
      o,n,r(i);
}
basis=vtz
distances=[0.4,0.41,0.42,0.43,0.44,0.45,0.46,0.47,0.48,0.49,0.5,0.51,0.52,0.53,0.54,0.55,0.56,0.57,0.58,0.59,0.6,0.61,0.62,0.63,0.64,0.65,0.66,0.67,0.68,0.69,0.7,0.71,0.72,0.73,0.74,0.75,0.76,0.77,0.78,0.79,0.8,0.81,0.82,0.83,0.84,0.85,0.86,0.87,0.88,0.89,0.9,0.91,0.92,0.93,0.94,0.95,0.96,0.97,0.98,0.99,1.0,1.01,1.02,1.03,1.04,1.05,1.06,1.07,1.08,1.09,1.1,1.11,1.12,1.13,1.14,1.15,1.16,1.17,1.18,1.19,1.2,1.21,1.22,1.23,1.24,1.25,1.26,1.27,1.28,1.29,1.3,1.31,1.32,1.33,1.34,1.35,1.36,1.37,1.38,1.39,1.4,1.41,1.42,1.43,1.44,1.45,1.46,1.47,1.48,1.49,1.5,1.51,1.52,1.53,1.54,1.55,1.56,1.57,1.58,1.59,1.6,1.61,1.62,1.63,1.64,1.65,1.66,1.67,1.68,1.69,1.7,1.71,1.72,1.73,1.74,1.75,1.76,1.77,1.78,1.79,1.8,1.81,1.82,1.83,1.84,1.85,1.86,1.87,1.88,1.89,1.9,1.91,1.92,1.93,1.94,1.95,1.96,1.97,1.98,1.99,2.0,2.01,2.02,2.03,2.04,2.05,2.06,2.07,2.08,2.09,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]
i=0
do ir=1,#distances
i=i+1
r(i)=distances(ir)

{hf;
wf,15,1,1
occ,4,2,2
}
escf(i)=energy

enddo

{table,r,escf
title,results for NO
}
---

回復(fù)此樓

» 猜你喜歡

PDF在線壓縮軟件已經(jīng)有5人回復(fù)
黑龍江省原子能研究院化學(xué)課題組招收2026級(jí)化學(xué)專業(yè)或化學(xué)相關(guān)專業(yè)碩士生已經(jīng)有0人回復(fù)
物理化學(xué)論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有139人回復(fù)
材料/化學(xué)相關(guān)專業(yè)2026級(jí)學(xué)術(shù)/專業(yè)型碩士研究生已經(jīng)有19人回復(fù)
求催化作用基礎(chǔ)第三版（甄開(kāi)吉) 課件已經(jīng)有3人回復(fù)
26年材料/化學(xué)相關(guān)專業(yè)博士研究生招生溫州大學(xué)化學(xué)與材料工程學(xué)院(26年秋季入學(xué)) 已經(jīng)有0人回復(fù)
求助電催化劑制備，怎么在碳布表面生長(zhǎng)PdCo合金已經(jīng)有0人回復(fù)
求助電催化劑制備，怎么在碳布表面生長(zhǎng)PdCo合金已經(jīng)有0人回復(fù)
304求調(diào)劑已經(jīng)有3人回復(fù)
上海華東理工大學(xué)機(jī)械與動(dòng)力工程學(xué)院先進(jìn)材料與制造研究組招聘公告已經(jīng)有0人回復(fù)

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

圖形求助，大神幫忙~~~ 已經(jīng)有19人回復(fù)
求大神幫忙翻譯已經(jīng)有2人回復(fù)
哪位大神能幫我把這個(gè)圖畫(huà)到word里哈？不是用截圖的形式哦~然后用word回復(fù)給我！已經(jīng)有5人回復(fù)
急急急急急急，TEM測(cè)試的圖太模糊了，求大神幫忙分析原因！已經(jīng)有19人回復(fù)
請(qǐng)大神幫忙看看這是什么植物？？？？？已經(jīng)有6人回復(fù)
新手求助~~求大神幫忙��！已經(jīng)有3人回復(fù)
積分求助求大神幫忙！已經(jīng)有7人回復(fù)
哪位大神幫幫忙已經(jīng)有4人回復(fù)
很簡(jiǎn)單的Fortran賦值出錯(cuò)啦，求大神們幫忙找bug... 已經(jīng)有11人回復(fù)
請(qǐng)各位大神幫忙看看�。�！急急已經(jīng)有7人回復(fù)
求大神幫忙解答！已經(jīng)有8人回復(fù)
請(qǐng)問(wèn)有沒(méi)有哪位環(huán)境類的大神做分子模擬之類的已經(jīng)有5人回復(fù)
熱重-差熱分析的圖，哪個(gè)大神幫忙解釋一下唄…… 已經(jīng)有6人回復(fù)
修改圖標(biāo)問(wèn)題，請(qǐng)各位大神幫忙！已經(jīng)有5人回復(fù)
問(wèn)題解決，謝謝關(guān)注！已經(jīng)有20人回復(fù)
第一次投英文的，作者指導(dǎo)有點(diǎn)沒(méi)看明白，求大神幫忙，謝謝已經(jīng)有17人回復(fù)
求兩個(gè)值中最大的值但編寫(xiě)的代碼不知道哪里錯(cuò)了，求大神幫忙。已經(jīng)有10人回復(fù)
求大神幫忙，怎么用tecplot出這種圖已經(jīng)有6人回復(fù)
求大神幫我改一下下面的程序！matlab迭代法的不知道錯(cuò)在哪里！已經(jīng)有7人回復(fù)

1樓 2014-05-06 18:36:56

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

清新俊逸

木蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 23 (小學(xué)生)
金幣: 415.6
散金: 5412
紅花: 15
帖子: 326
在線: 331.9小時(shí)
蟲(chóng)號(hào): 539264
注冊(cè): 2008-04-04
專業(yè): 化學(xué)動(dòng)力學(xué)

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

你這里有兩個(gè)問(wèn)題。
1、Molpro中，內(nèi)坐標(biāo)的鍵長(zhǎng)單位默認(rèn)是Bohr，你小于2Bohr的算了那么多都沒(méi)用的，只要看2到5Bohr的就差不多了。
2、你說(shuō)的沒(méi)有一個(gè)阱的結(jié)構(gòu)，是因?yàn)槟愕膽B(tài)算錯(cuò)了，你算的這個(gè)態(tài)不是成鍵的態(tài)，所以鍵長(zhǎng)越長(zhǎng)，能量越低。
解決辦法有兩個(gè)，
一是不使用對(duì)稱性，算出來(lái)的就是最低的態(tài)。在geom上一行加上nosym，然后occ后面的參數(shù)改成8就行了。需要注意的是，這種方法不一定能得到正確的結(jié)果，需要旋轉(zhuǎn)軌道嘗試很多次。
二是把態(tài)給選對(duì)。N-O的最大的abel群是c2v群，有4個(gè)不可約表示，需要搞清楚每個(gè)不可約表示的原子軌道占據(jù)數(shù)，以及波函數(shù)屬于哪個(gè)不可約表示。
我用wf,15,1,1; occ,5,2,1,0或者wf,15,2,1; occ,5,2,1,0算出來(lái)都比你原來(lái)的要低，你的occ,4,2,2肯定是不對(duì)的。

[ 發(fā)自手機(jī)版 http://www.gaoyang168.com/3g ]

贊一下

回復(fù)此樓

6樓2014-05-09 21:49:48

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

查看全部 8 個(gè)回答

luqing6879

木蟲(chóng) (著名寫(xiě)手)

應(yīng)助: 8 (幼兒園)
金幣: 3567.5
散金: 249
紅花: 3
帖子: 1907
在線: 808.1小時(shí)
蟲(chóng)號(hào): 614833
注冊(cè): 2008-09-28
專業(yè): 理論和計(jì)算化學(xué)

★ ★
gmy1990: 金幣+2 2014-05-07 05:34:17

第一，你的循環(huán)是有問(wèn)題的；不需要i=i+1
第二，你開(kāi)始給的鍵長(zhǎng)太小了。molpro默認(rèn)的單位是bohr。

贊一下

回復(fù)此樓

cishan.goodweb.cn/慈善點(diǎn)擊，點(diǎn)滴幫助。

2樓2014-05-07 01:01:49

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

kdswx

新蟲(chóng) (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1799
帖子: 167
在線: 47小時(shí)
蟲(chóng)號(hào): 2736957
注冊(cè): 2013-10-19

引用回帖:

2樓: Originally posted by luqing6879 at 2014-05-07 01:01:49
第一，你的循環(huán)是有問(wèn)題的；不需要i=i+1
第二，你開(kāi)始給的鍵長(zhǎng)太小了。molpro默認(rèn)的單位是bohr。

我這樣的程序能出來(lái)數(shù)據(jù)，可問(wèn)題是數(shù)據(jù)畫(huà)出圖來(lái)是有能量隨著鍵長(zhǎng)減小然后趨于平緩，沒(méi)有再增加的趨勢(shì)，不符合勢(shì)能曲線啊，這是為什么呢？

贊一下

回復(fù)此樓

3樓2014-05-07 08:17:26

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

luqing6879

木蟲(chóng) (著名寫(xiě)手)

應(yīng)助: 8 (幼兒園)
金幣: 3567.5
散金: 249
紅花: 3
帖子: 1907
在線: 808.1小時(shí)
蟲(chóng)號(hào): 614833
注冊(cè): 2008-09-28
專業(yè): 理論和計(jì)算化學(xué)

引用回帖:

3樓: Originally posted by kdswx at 2014-05-07 08:17:26
我這樣的程序能出來(lái)數(shù)據(jù)，可問(wèn)題是數(shù)據(jù)畫(huà)出圖來(lái)是有能量隨著鍵長(zhǎng)減小然后趨于平緩，沒(méi)有再增加的趨勢(shì)，不符合勢(shì)能曲線啊，這是為什么呢？...

看看你的scf有沒(méi)有收斂吧。

贊一下

回復(fù)此樓

cishan.goodweb.cn/慈善點(diǎn)擊，點(diǎn)滴幫助。

4樓2014-05-07 08:44:42

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

查看全部 8 個(gè)回答

普通表情龍兔虎貓

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

[考研] 307求調(diào)劑 +3	冷笙123 2026-03-17	3/150	2026-03-18 09:55 by macy2011
[考研] 293求調(diào)劑 +11	zjl的號(hào) 2026-03-16	16/800	2026-03-18 08:10 by zhukairuo
[考研] 344求調(diào)劑 +4	knight344 2026-03-16	4/200	2026-03-17 17:27 by ruiyingmiao
[考研] 【0856】化學(xué)工程（085602）313 分，本科學(xué)科評(píng)估A類院�；瘜W(xué)工程與工藝，誠(chéng)求調(diào)劑 +7	小劉快快上岸 2026-03-11	8/400	2026-03-17 16:57 by ruiyingmiao
[考研] 085600材料與化工求調(diào)劑 +5	緒幸與子 2026-03-17	5/250	2026-03-17 16:40 by laoshidan
[考研] 有沒(méi)有道鐵/土木的想調(diào)劑南林，給自己招師弟中～ +3	TqlXswl 2026-03-16	7/350	2026-03-17 15:23 by TqlXswl
[考研] 285化工學(xué)碩求調(diào)劑（081700） +9	柴郡貓_ 2026-03-12	9/450	2026-03-17 10:18 by Sammy2
[考研] 化學(xué)調(diào)劑0703 +8	啊我我的 2026-03-11	8/400	2026-03-16 17:23 by 我的船我的海
[考研] 0856求調(diào)劑 +3	劉夢(mèng)微 2026-03-15	3/150	2026-03-16 10:00 by houyaoxu
[考研] 297求調(diào)劑 +4	學(xué)海漂泊 2026-03-13	4/200	2026-03-14 11:51 by 熱情沙漠
[考研] 材料工程調(diào)劑 +9	咪咪空空 2026-03-12	9/450	2026-03-13 22:05 by 星空星月
[考研] 0856材料與化工301求調(diào)劑 +5	奕束光 2026-03-13	5/250	2026-03-13 22:00 by 星空星月
[考研] 求b區(qū)學(xué)校調(diào)劑 +3	周56 2026-03-11	3/150	2026-03-13 16:20 by JourneyLucky
[考研] 310求調(diào)劑 +3	【上上簽】 2026-03-11	3/150	2026-03-13 16:16 by JourneyLucky
[考研] 工科278分求調(diào)劑 +5	周慢熱啊 2026-03-12	7/350	2026-03-13 15:49 by JourneyLucky
[考研] 工科材料085601 279求調(diào)劑 +8	困于星晨 2026-03-12	10/500	2026-03-13 15:42 by ms629
[考研] 314求調(diào)劑 +7	無(wú)懈可擊的巨人 2026-03-12	7/350	2026-03-13 15:40 by JourneyLucky
[考研] 一志愿山大07化學(xué) 332分四六級(jí)已過(guò) 本科山東雙非求調(diào)劑！ +3	不想理你 2026-03-12	3/150	2026-03-13 14:18 by JourneyLucky
[考研] 0856化學(xué)工程280分求調(diào)劑 +4	shenzxsn 2026-03-11	4/200	2026-03-13 11:55 by ymwdoctor
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +3	崔wj 2026-03-12	4/200	2026-03-12 19:33 by 求調(diào)劑zz