| 8 | 1/1 | 返回列表 |
| 查看: 1741 | 回復: 7 | ||
[求助]
求大神們幫忙看看這個程序哪兒需要改啊 已有1人參與
|
|
初學molpro,學著計算NO基態(tài)的勢能,想要單點能掃描,這下面是計算程序,但是出來的結(jié)果不對啊,能不能請大神們幫忙看看,萬分感謝! ***,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 } --- |
木蟲 (著名寫手)

木蟲 (著名寫手)

|
你這里有兩個問題。 1、Molpro中,內(nèi)坐標的鍵長單位默認是Bohr,你小于2Bohr的算了那么多都沒用的,只要看2到5Bohr的就差不多了。 2、你說的沒有一個阱的結(jié)構(gòu),是因為你的態(tài)算錯了,你算的這個態(tài)不是成鍵的態(tài),所以鍵長越長,能量越低。 解決辦法有兩個, 一是不使用對稱性,算出來的就是最低的態(tài)。在geom上一行加上nosym,然后occ后面的參數(shù)改成8就行了。需要注意的是,這種方法不一定能得到正確的結(jié)果,需要旋轉(zhuǎn)軌道嘗試很多次。 二是把態(tài)給選對。N-O的最大的abel群是c2v群,有4個不可約表示,需要搞清楚每個不可約表示的原子軌道占據(jù)數(shù),以及波函數(shù)屬于哪個不可約表示。 我用wf,15,1,1; occ,5,2,1,0或者wf,15,2,1; occ,5,2,1,0算出來都比你原來的要低,你的occ,4,2,2肯定是不對的。 [ 發(fā)自手機版 http://www.gaoyang168.com/3g ] |
|
1.5 -128.6871215 1.6 -128.9296432 1.7 -129.0904337 1.8 -129.1926184 1.9 -129.2528221 2.0 -129.2830743 2.1 -129.2921377 2.2 -129.2864358 2.3 -129.2707005 2.4 -129.2484223 2.5 -129.2221507 2.6 -129.1937002 2.7 -129.1643240 2.8 -129.1348963 2.9 -129.1061222 3.0 -129.0788913 3.1 -129.0547177 3.2 -129.0348454 3.3 -129.0191149 3.4 -129.0067277 3.5 -128.9968953 3.6 -128.9889930 3.8 -128.9772466 4.0 -128.9690866 4.2 -128.9632158 4.4 -128.9588936 4.6 -128.9556628 4.8 -128.9532197 5.0 -128.9513529 [ 發(fā)自手機版 http://www.gaoyang168.com/3g ] |
| 8 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 能源材料化學課題組招收碩士研究生8-10名 +3 | 脫穎而出 2026-03-16 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 材料專碩326求調(diào)劑 +6 | 墨煜姒莘 2026-03-15 | 7/350 |
|
|
[考研] 312求調(diào)劑 +4 | 陌宸希 2026-03-16 | 5/250 |
|
|
[考研] 08工科 320總分 求調(diào)劑 +4 | 梨花珞晚風 2026-03-17 | 4/200 |
|
|
[考博] 26申博 +4 | 八旬速覽 2026-03-16 | 4/200 |
|
|
[考研] 304求調(diào)劑 +4 | ahbd 2026-03-14 | 4/200 |
|
|
[考研] 0703化學調(diào)劑 +6 | 妮妮ninicgb 2026-03-15 | 9/450 |
|
|
[考研] 085600材料與化工 求調(diào)劑 +13 | enenenhui 2026-03-13 | 14/700 |
|
|
[考研] 中科大材料專碩319求調(diào)劑 +3 | 孟鑫材料 2026-03-13 | 3/150 |
|
|
[考研] 308 085701 四六級已過求調(diào)劑 +7 | 溫喬喬喬喬 2026-03-12 | 14/700 |
|
|
[考研] 材料與化工(0856)304求B區(qū)調(diào)劑 +6 | 邱gl 2026-03-12 | 7/350 |
|
|
[考研] 336求調(diào)劑 +6 | Iuruoh 2026-03-11 | 6/300 |
|
|
[考研] [0860]321分求調(diào)劑,ab區(qū)皆可 +4 | 寶貴熱 2026-03-13 | 4/200 |
|
|
[考研] 工科278分求調(diào)劑 +5 | 周慢熱啊 2026-03-12 | 7/350 |
|
|
[考研] 工科材料085601 279求調(diào)劑 +8 | 困于星晨 2026-03-12 | 10/500 |
|
|
[考研] 290求調(diào)劑 +7 | ADT 2026-03-12 | 7/350 |
|
|
[考研] 工科調(diào)劑 +4 | Jiang191123! 2026-03-11 | 4/200 |
|
|
[考研] 289求調(diào)劑 +3 | 李政瑩 2026-03-12 | 3/150 |
|
|
[考研] 070303一志愿西北大學學碩310找調(diào)劑 +3 | d如愿上岸 2026-03-13 | 3/150 |
|
|
[考研] 283求調(diào)劑,材料、化工皆可 +8 | 蘇打水7777 2026-03-11 | 10/500 |
|