| 8 | 1/1 | 返回列表 |
| 查看: 1740 | 回復: 7 | ||
[求助]
求大神們幫忙看看這個程序哪兒需要改啊 已有1人參與
|
|
初學molpro,學著計算NO基態(tài)的勢能,想要單點能掃描,這下面是計算程序,但是出來的結果不對啊,能不能請大神們幫忙看看,萬分感謝! ***,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、你說的沒有一個阱的結構,是因為你的態(tài)算錯了,你算的這個態(tài)不是成鍵的態(tài),所以鍵長越長,能量越低。 解決辦法有兩個, 一是不使用對稱性,算出來的就是最低的態(tài)。在geom上一行加上nosym,然后occ后面的參數(shù)改成8就行了。需要注意的是,這種方法不一定能得到正確的結果,需要旋轉軌道嘗試很多次。 二是把態(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ā)表 | |
|---|---|---|---|---|
|
[考研] 268求調(diào)劑 +3 | 簡單點0 2026-03-17 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 考研求調(diào)劑 +3 | 橘頌. 2026-03-17 | 4/200 |
|
|
[考研] 301求調(diào)劑 +4 | A_JiXing 2026-03-16 | 4/200 |
|
|
[考研] 【0856】化學工程(085602)313 分,本科學科評估A類院;瘜W工程與工藝,誠求調(diào)劑 +7 | 小劉快快上岸 2026-03-11 | 8/400 |
|
|
[考研] 梁成偉老師課題組歡迎你的加入 +8 | 一鴨鴨喲 2026-03-14 | 10/500 |
|
|
[考研] 289求調(diào)劑 +6 | 步川酷紫123 2026-03-11 | 6/300 |
|
|
[考研] 藥學383 求調(diào)劑 +3 | 藥學chy 2026-03-15 | 4/200 |
|
|
[考研] 環(huán)境工程調(diào)劑 +6 | 大可digkids 2026-03-16 | 6/300 |
|
|
[考研] 326求調(diào)劑 +4 | 諾貝爾化學獎覬?/a> 2026-03-15 | 7/350 |
|
|
[考研] 304求調(diào)劑 +4 | ahbd 2026-03-14 | 4/200 |
|
|
[考研] 304求調(diào)劑 +3 | 曼殊2266 2026-03-14 | 3/150 |
|
|
[考研] 中科大材料與化工319求調(diào)劑 +3 | 孟鑫材料 2026-03-14 | 3/150 |
|
|
[考研] 材料與化工 323 英一+數(shù)二+物化,一志愿:哈工大 本人本科雙一流 +4 | 自由的_飛翔 2026-03-13 | 5/250 |
|
|
[考研] 學碩285求調(diào)劑 +13 | Wisjxn 2026-03-12 | 46/2300 |
|
|
[考研] 330求調(diào)劑 +3 | ?醬給調(diào)劑跪了 2026-03-13 | 3/150 |
|
|
[考研] 一志愿中科院,化學方向,295求調(diào)劑 +4 | 一氧二氮 2026-03-11 | 4/200 |
|
|
[考研] 290求調(diào)劑 +9 | ADT 2026-03-11 | 9/450 |
|
|
[考研] 290求調(diào)劑 +7 | ADT 2026-03-12 | 7/350 |
|
|
[考研] 328化工專碩求調(diào)劑 +4 | 。,。,。,。i 2026-03-12 | 4/200 |
|
|
[考研] 321求調(diào)劑(食品/專碩) +3 | xc321 2026-03-12 | 6/300 |
|