| 查看: 2716 | 回復(fù): 11 | ||||
科學(xué)島鐵蟲 (初入文壇)
|
[求助]
用NBO如何分析體系,特別是開殼層體系,
|
|||
|
最近瀏覽了很多關(guān)于用guassian計算NBO的帖子,但是對一些問題還是很迷茫, 1.開殼層體系的鍵級:有人說用總密度矩陣得到的wiberg 鍵級,也有人說用alpha和beta得到的鍵級之和的2倍,更有人說是alpha和beta的一半,現(xiàn)在我也不知道怎么處理了? 2.關(guān)于鍵級的離子鍵和共價鍵成分,據(jù)說NBO3.1版本不能計算NRT(自然共振理論),假如是這樣的話,那離子鍵成分能不能表示為i=abs(C(A)^2-C(B)^2)/(C(A)^2+C(B)^2),i=1為純離子鍵,i=0為純共價鍵? 3.那怎樣才能算NRT,我最近看了NBO的手冊,不是很明白,另外我也看了NRT的3篇文獻(xiàn), Natural Resonance Theory: I. Natural Bond Order and Valency Natural Resonance Theory: II. NaturalBond Order and Valency Natural Resonance Theory: III. NaturalBond Order and Valency 是否有單獨的NBO計算軟件,現(xiàn)在我想用高版本的軟件計算NRT,現(xiàn)在求軟件!。 有誰可以指點指點? 我就只有19個金幣,全部懸賞,還請大家多多指點啊,非常感謝。。。。。。。! |
納米金 | Gaussian |
鐵桿木蟲 (著名寫手)
至尊木蟲 (著名寫手)
本人已永久離開小木蟲
|
1 可以將Wiberg鍵級當(dāng)成Mayer鍵級在正交基下的特例(即重疊矩陣為單位矩陣的情況)。閉殼層、開殼層下的Mayer鍵級表達(dá)式在Multiwfn手冊3.11.1節(jié)明確地給出了。 靠總密度矩陣來得到開殼層體系的Wiberg/Mayer鍵級在原理上是不正確的,雖然結(jié)果從實際上來看也可以用。正確做法是靠alpha和beta密度矩陣分別計算alpha和beta電子的Wiberg/Mayer鍵級,然后再加和。 2 可以。但這只是一種表述方式,是否可靠得自己根據(jù)實際情況定奪,某些情況可能沒法這么討論。相對而言,靠NRT給出離子和共價成分更可靠。 3 用NBO5及以上版本(收費)。高斯自帶的3.1不支持。 |
至尊木蟲 (著名寫手)
本人已永久離開小木蟲
|
1 計算alpha和beta在加和,就是手冊里3.11.1節(jié)里第一個公式,這也是Mayer的原文所定義的。你引用的那句話只不過是另一種處理開殼層體系的方式。 2 計算Mayer/Wiberg鍵級時不要用彌散函數(shù)。你可以計算能量等性質(zhì)時用彌散函數(shù),但是計算鍵級時使用不加彌散函數(shù)的波函數(shù)。也可以改用模糊鍵級,原理上類似Mayer鍵級,但是是在實空間下計算的,所以受彌散函數(shù)影響很小,用彌散函數(shù)也沒關(guān)系,詳見手冊3.11.6節(jié),對于極性不強的鍵和Mayer鍵級在結(jié)果上接近。 3 NRT是很耗時的。比較大的體系往往很難算得動。 |
鐵蟲 (初入文壇)
鐵蟲 (初入文壇)
鐵蟲 (初入文壇)
|
1.我看了3.11.1章節(jié),是這么說的: For unrestricted or restricted open-shell system, there is another way to calculate total bond order rather than summing up alpha and beta bond orders, that is summing up alpha and beta density matrices to form total density matrix first and then calculate Mayer bond order by using restricted close-shell formula, this treatment is sometimes called “generalized Wiberg bond order“, these total bond orders are printed following the title “Mayer bond order from mixed alpha&beta density matrix”. 靠alpha和beta密度矩陣分別計算alpha和beta電子的Wiberg/Mayer鍵級,然后再加和 和章節(jié)里說的好像不一致,是我理解錯了嗎? 2.我算的是AuC的陰離子體系,按說最好用彌散基組吧,但手冊里又說:Similar to Mulliken population, Mayer bond order and the multi-center bond order described below are sensitive to basis-set, so don't use the basis-sets having diffusion functions, otherwise the bond order result will be unreliable.我現(xiàn)在又迷惑了。 3.NRT計算所需的時間長嗎?我只有NBO3.1版本的,所以先問問, |
鐵蟲 (初入文壇)
鐵蟲 (初入文壇)
至尊木蟲 (著名寫手)
本人已永久離開小木蟲
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿天津大學(xué)化學(xué)工藝專業(yè)(081702)315分求調(diào)劑 +10 | yangfz 2026-03-17 | 10/500 |
|
|---|---|---|---|---|
|
[考研] 一志愿985,本科211,0817化學(xué)工程與技術(shù)319求調(diào)劑 +7 | Liwangman 2026-03-15 | 7/350 |
|
|
[考研] 26調(diào)劑/材料/英一數(shù)二/總分289/已過A區(qū)線 +7 | 步川酷紫123 2026-03-13 | 7/350 |
|
|
[考研] 297求調(diào)劑 +8 | 戲精丹丹丹 2026-03-17 | 8/400 |
|
|
[考研] 302求調(diào)劑 +10 | 呼呼呼。。。。 2026-03-17 | 10/500 |
|
|
[考研] 312求調(diào)劑 +8 | 陌宸希 2026-03-16 | 9/450 |
|
|
[考研] 0703化學(xué)336分求調(diào)劑 +6 | zbzihdhd 2026-03-15 | 7/350 |
|
|
[基金申請] 被我言中:新模板不強調(diào)格式了,假專家開始管格式了 +4 | beefly 2026-03-14 | 4/200 |
|
|
[考研] 268求調(diào)劑 +7 | 好運連綿不絕 2026-03-12 | 8/400 |
|
|
[考研] 268求調(diào)劑 +8 | 一定有學(xué)上- 2026-03-14 | 9/450 |
|
|
[考研] 本人考085602 化學(xué)工程 專碩 +16 | 不知道叫什么! 2026-03-15 | 18/900 |
|
|
[考研] 332求調(diào)劑 +6 | Zz版 2026-03-13 | 6/300 |
|
|
[考研] 機械專碩325,尋找調(diào)劑院校 +3 | y9999 2026-03-15 | 5/250 |
|
|
[考研] 0703 物理化學(xué)調(diào)劑 +3 | 我可以上岸的對?/a> 2026-03-13 | 5/250 |
|
|
[考研] 中科大材料與化工319求調(diào)劑 +3 | 孟鑫材料 2026-03-14 | 3/150 |
|
|
[考研] 255求調(diào)劑 +3 | 李嘉慧, 2026-03-12 | 4/200 |
|
|
[考研] 材料080500調(diào)劑求收留 +3 | 一顆meteor 2026-03-13 | 3/150 |
|
|
[考研] 26調(diào)劑/材料科學(xué)與工程/總分295/求收留 +9 | 2026調(diào)劑俠 2026-03-12 | 9/450 |
|
|
[考研] 308求調(diào)劑 +3 | 是Lupa啊 2026-03-12 | 3/150 |
|
|
[論文投稿]
投稿問題
5+4
|
星光燦爛xt 2026-03-12 | 6/300 |
|