| 5 | 1/1 | 返回列表 |
| 查看: 2118 | 回復(fù): 4 | ||||
[交流]
【求助】如何讀取輸出文件中的鍵級
|
||||
高斯日常學(xué)習(xí)總結(jié) | Gaussian |
» 搶金幣啦!回帖就可以得到:
+2/100
+1/91
+1/86
+1/85
+1/85
+1/53
+1/31
+1/30
+1/19
+1/16
+1/14
+1/5
+1/5
+1/4
+1/2
+1/2
+1/1
+1/1
+1/1
+1/1
|
引用Sobereva姐的文章:談?wù)凬BO中用bndidx關(guān)鍵詞輸出的三種鍵級 希望對你有所幫助!! 1. 前言 用了bndidx關(guān)鍵詞后,NBO程序會進(jìn)行鍵級分析,輸出三種鍵級矩陣,包括: 1 Wiberg bond index matrix in the NAO basis 2 Atom-atom overlap-weighted NAO bond order 3 MO bond order 矩陣中(i,j)元素對應(yīng)i、j原子間的鍵級。對這三種不同鍵級,有些人對應(yīng)當(dāng)選取哪個感到疑惑,此文將談?wù)勥@三種鍵級的含義及其合理性。值得一提的是,經(jīng)常有人說“NBO鍵級”,這種說法是很不嚴(yán)格的,在NBO程序中,上面已經(jīng)看到有三種鍵級了,NBO還能產(chǎn)生LMO鍵級、NRT鍵級,另外,NBO程序搜索出的NBO們的占據(jù)數(shù)本身也可以用于衡量鍵級。所以,不應(yīng)當(dāng)使用含糊不清的“NBO鍵級”這個詞。 為了簡明且便于驗證,本文例子使用HF/STO-3G的基態(tài)水分子波函數(shù)。下面的是用DMNAO關(guān)鍵詞輸出的密度矩陣,將在后文用到。后文直接將此矩陣稱為DMNAO。 NAO 1 2 3 4 5 6 7 ---------- ------- ------- ------- ------- ------- ------- ------- 1. O 1 (S) 2.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 2. O 1 (S) 0.0000 1.7325 0.0000 0.0000 0.3538 0.4113 0.4113 3. O 1 (px) 0.0000 0.0000 2.0000 0.0000 0.0000 0.0000 0.0000 4. O 1 (py) 0.0000 0.0000 0.0000 1.1399 0.0000 0.7001 -0.7001 5. O 1 (pz) 0.0000 0.3538 0.0000 0.0000 1.5321 -0.5439 -0.5439 6. H 2 (S) 0.0000 0.4113 0.0000 0.7001 -0.5439 0.7977 -0.0623 7. H 3 (S) 0.0000 0.4113 0.0000 -0.7001 -0.5439 -0.0623 0.7977 2. 第一種鍵級(Wiberg bond index matrix in the NAO basis): 說起第一種鍵級,先要回顧一下Wiberg鍵級和Mayer鍵級的原始定義。Wiberg鍵級最早來自Tetrahedron 24,1083,鍵級W_AB=∑[i∈A]∑[j∈B]P(i,j)^2,其中P為密度矩陣,這個定義只適合正交基。Mayer在CPL 97,270里將其擴(kuò)展到非正交基,A、B間的鍵級定義為∑[i∈A]∑[j∈B]PS(i,j)*PS(j,i),這里PS(i,j)是指P和重疊矩陣S做矩陣乘所得矩陣的(i,j)元素。當(dāng)使用正交基后,變?yōu)椤芠i∈A]∑[j∈B]P(i,j)*P(j,i),由于P是對稱的,所以又還原到了Wiberg鍵級。Mayer鍵級數(shù)值比較符合一般化學(xué)觀念,比如單重鍵約為1,三重鍵約為3,苯的C-C鍵約為1.5。 NBO的各種分析主要基于由原始非正交基函數(shù)轉(zhuǎn)化成的自然原子軌道(NAO),NAO是正交基,所以可以直接使用Wiberg鍵級進(jìn)行分析,也就是NBO輸出的Wiberg bond index matrix in the NAO basis。這是三種鍵級中唯一推薦的鍵級。NBO給出的結(jié)果: Atom 1 2 3 ---- ------ ------ ------ 1. O 0.0000 0.9552 0.9552 2. H 0.9552 0.0000 0.0039 3. H 0.9552 0.0039 0.0000 可以很容易地通過NAO的密度矩陣DMNAO來計算此矩陣。此矩陣的(2,3)代表H-H的鍵級,由于6、7號NAO分別是兩個氫上的唯一基函數(shù),鍵級也就是DMNAO(6,7)^2=(-0.0623)^2=0.00388。再比如(1,3)代表的O-H鍵級,由于1~5號基函數(shù)對應(yīng)O,也就是∑[i=1,2,3,4,5]DMNAO(i,7)^2=0.95513。 在NBO輸出這個矩陣之后還會輸出Wiberg bond index, Totals by atom,指的就是每個原子形成鍵級的總和,也稱為原子價。對于第二種、第三種鍵級也是如此。 3. 第二種鍵級(Atom-atom overlap-weighted NAO bond order) 這種鍵級號稱是考慮了原子軌道間重疊,引入距離因素,并使得結(jié)果能看出成鍵反鍵,但實際結(jié)果不好,而且其定義我認(rèn)為沒有任何根據(jù)。A、B間鍵級定義為∑[i∈A]∑[j∈B]P(i,j)*S'(i,j),其中S'是PNAO重疊矩陣(注:PNAO是原始基函數(shù)向NAO轉(zhuǎn)換的中間基函數(shù),PNAO在原子內(nèi)正交而在原子間未必正交)。這種定義很莫名其妙,密度矩陣用NAO的,但重疊矩陣卻用PNAO的。從其形式上,可以想起Mulliken鍵級。Mulliken鍵級就是Mulliken布居分析得到的兩原子間交叉項的總和,∑[imo∈占據(jù)]∑[i∈A]∑[j∈B]2*占據(jù)數(shù)(imo)*C(i,imo)*C(j,imo)*S(i,j),寫成密度矩陣的形式就是∑[i∈A]∑[j∈B]2*P(i,j)*S(i,j)。所以,第二種鍵級實際上是一種既不完全是NAO、也不完全是PNAO下的不倫不類的Mulliken鍵級,而且還少乘了2。這種古怪的鍵級缺乏依據(jù),實際結(jié)果又差,所以不要使用。NBO輸出的結(jié)果為: Atom 1 2 3 ---- ------ ------ ------ 1. O 0.0000 0.5399 0.5399 2. H 0.5399 0.0000 -0.0145 3. H 0.5399 -0.0145 0.0000 這也容易驗證,使用SPNAO關(guān)鍵詞可以輸出PNAO重疊矩陣S' NAO 1 2 3 4 5 6 7 ---------- ------- ------- ------- ------- ------- ------- ------- 1. O 1 (S) 1.0000 0.0000 0.0000 0.0000 0.0000 0.0935 0.0935 2. O 1 (S) 0.0000 1.0000 0.0000 0.0000 0.0000 0.4675 0.4675 3. O 1 (px) 0.0000 0.0000 1.0000 0.0000 0.0000 0.0000 0.0000 4. O 1 (py) 0.0000 0.0000 0.0000 1.0000 0.0000 0.3205 -0.3205 5. O 1 (pz) 0.0000 0.0000 0.0000 0.0000 1.0000 -0.2265 -0.2265 6. H 2 (S) 0.0935 0.4675 0.0000 0.3205 -0.2265 1.0000 0.2328 7. H 3 (S) 0.0935 0.4675 0.0000 -0.3205 -0.2265 0.2328 1.0000 設(shè)T(i,j)=DMNAO(i,j)*S'(i,j),則T為: 1 2 3 4 5 6 7 1 2.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 2 0.00000 1.73250 0.00000 0.00000 0.00000 0.19228 0.19228 3 0.00000 0.00000 2.00000 0.00000 0.00000 0.00000 0.00000 4 0.00000 0.00000 0.00000 1.13990 0.00000 0.22438 0.22438 5 0.00000 0.00000 0.00000 0.00000 1.53210 0.12319 0.12319 6 0.00000 0.19228 0.00000 0.22438 0.12319 0.79770 -0.01450 7 0.00000 0.19228 0.00000 0.22438 0.12319 -0.01450 0.79770 求O-H鍵的這種鍵級也就是∑[i=1,2,3,4,5]T(i,7)=0+0.19228+0+0.22438+0.12319=0.53985,和NBO輸出的一致。 另外值得一提的是,所謂NPA電荷(無數(shù)人謬稱為NBO電荷),就是以NAO為基獲得的Mulliken電荷。A原子布居數(shù)=∑[i∈A]∑[j]P(i,j)S(i,j)=∑[i∈A]∑[j]P(i,j)δ(i,j)=∑[i∈A]P(i,i),所以氧的布居數(shù)∑[i=1,2,3,4,5]DMNAO(i,i)=2.0000+1.7325+2.0000+1.1399+1.5321=8.4045,凈電荷為8-8.4045=-0.4045,和NBO程序輸出的一致。 4. 第三種鍵級(MO bond order) 這個鍵級的名字十分令人疑惑,而且手冊上也沒說明。實際上它的定義極其簡單,A、B間的鍵級為∑[i∈A]∑[j∈B]P(i,j)。NBO輸出為 Atom 1 2 3 ---- ------ ------ ------ 1. O 0.0000 0.5675 -0.8328 2. H 0.5675 0.0000 -0.0623 3. H -0.8328 -0.0623 0.0000 求O-H2鍵級就是∑[i=1,2,3,4,5]DMNAO(i,7)=0.4113-0.7001-0.5439=-0.8327。結(jié)果十分荒誕,明明是典型sigma鍵,鍵級卻為負(fù)。而另一個等價的O-H3鍵鍵級卻是正值0.4113+0.7001-0.5439=0.5675!可見,這種鍵級從實際應(yīng)用上看毫無意義。從理論上同樣沒有根據(jù),從形式上看和Mulliken鍵級相仿佛,卻少乘了重疊矩陣元,另外NAO又是正交基,不存在交叉項,也不可能求Mulliken鍵級。而與Wiberg鍵級相比,又少求了平方。我認(rèn)為NBO程序中應(yīng)該去掉這種毫無意義的鍵級。 5. 用PNAO計算Mayer鍵級的嘗試 原理上,以非正交基PNAO做Mayer鍵級分析也是可以的,可以代替正交基NAO下的Wiberg鍵級分析,另外,做Mulliken分析通過交叉項還可以研究成鍵反鍵和各分子軌道貢獻(xiàn),但是NBO程序里沒有提供。筆者嘗試?yán)肞NAO做Mayer鍵級分析,這首先需要獲得PNAO密度矩陣,NBO3.1里不支持輸出它的DMPNAO關(guān)鍵詞,故通過.47文件在GENNBO5.0下輸出,所得矩陣如下: NAO 1 2 3 4 5 6 7 ---------- ------- ------- ------- ------- ------- ------- ------- 1. O 1( 1s ) 2.0000 0.0000 0.0000 0.0000 0.0000 0.1870 0.1870 2. O 1( 2s ) 0.0000 1.9767 0.0000 0.0000 0.0759 1.0236 1.0236 3. O 1( 2px) 0.0000 0.0000 2.0000 0.0000 0.0000 0.0000 0.0000 4. O 1( 2py) 0.0000 0.0000 0.0000 1.5602 0.0000 0.9391 -0.9391 5. O 1( 2pz) 0.0000 0.0759 0.0000 0.0000 1.7530 -0.7415 -0.7415 6. H 2( 1s ) 0.1870 1.0236 0.0000 0.9391 -0.7415 1.4612 0.3307 7. H 3( 1s ) 0.1870 1.0236 0.0000 -0.9391 -0.7415 0.3307 1.4612 然而,不幸的是,所得Mayer鍵級矩陣很糟 1 2 3 1 0.00000000 7.23697119 7.23697119 2 7.23697119 0.00000000 1.06885215 3 7.23697119 1.06885215 0.00000000 鍵級大得出奇。NBO5.0輸出的PNAO密度矩陣估計錯誤的,理論上Tr(PS)應(yīng)當(dāng)?shù)扔诜肿涌傠娮訑?shù),然而用PNAO的密度矩陣和重疊矩陣竟得到16.22606912。雖然測試這個方法由于NBO程序的bug失敗了,但原理上是合理的。 |
| 5 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 265求調(diào)劑 +4 | 梁梁校校 2026-03-17 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 329求調(diào)劑 +5 | 想上學(xué)吖吖 2026-03-19 | 5/250 |
|
|
[考研] 307求調(diào)劑 +9 | 冷笙123 2026-03-17 | 9/450 |
|
|
[考研] 288求調(diào)劑 +15 | 于海海海海 2026-03-19 | 15/750 |
|
|
[考研] 梁成偉老師課題組歡迎你的加入 +9 | 一鴨鴨喲 2026-03-14 | 11/550 |
|
|
[考研] 317求調(diào)劑 +3 | 申子申申 2026-03-19 | 6/300 |
|
|
[考研] 化學(xué)求調(diào)劑 +3 | 臨澤境llllll 2026-03-17 | 4/200 |
|
|
[考研] 一志愿 西北大學(xué) ,070300化學(xué)學(xué)碩,總分287,雙非一本,求調(diào)劑。 +3 | 晨昏線與星海 2026-03-19 | 3/150 |
|
|
[考研] 332求調(diào)劑 +3 | ydfyh 2026-03-17 | 3/150 |
|
|
[教師之家] 焦慮 +9 | 水冰月月野兔 2026-03-13 | 13/650 |
|
|
[考研] 0703化學(xué)調(diào)劑 +3 | 妮妮ninicgb 2026-03-17 | 3/150 |
|
|
[考研] 工科材料085601 279求調(diào)劑 +6 | 困于星晨 2026-03-17 | 6/300 |
|
|
[考研] 生物學(xué)071000 329分求調(diào)劑 +3 | 我愛生物生物愛?/a> 2026-03-17 | 3/150 |
|
|
[考研] 277調(diào)劑 +5 | 自由煎餅果子 2026-03-16 | 6/300 |
|
|
[考研] 一志愿蘇州大學(xué)材料工程(085601)專碩有科研經(jīng)歷三項國獎兩個實用型專利一項省級立項 +6 | 大火山小火山 2026-03-16 | 8/400 |
|
|
[論文投稿] 有沒有大佬發(fā)小論文能帶我個二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|
|
[考研] 304求調(diào)劑 +3 | 曼殊2266 2026-03-14 | 3/150 |
|
|
[考研] 本科南京大學(xué)一志愿川大藥學(xué)327 +3 | 麥田耕者 2026-03-14 | 3/150 |
|
|
[考研] 289求調(diào)劑 +4 | 這么名字咋樣 2026-03-14 | 6/300 |
|
|
[考研] 304求調(diào)劑 +7 | 7712b 2026-03-13 | 7/350 |
|