 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 2483 | 回復(fù): 10 | ||
[求助]
Multiwfn讀取NBO的.40文件是對應(yīng)MO嗎？
|
|
如題，manual里說NBO的.40文件是NLMO，我聽說這種軌道的性質(zhì)和MO是相同的，請問這種說法是否正確？ 另外，我試著做了.40的電子密度格點文件，又對同一結(jié)構(gòu)(相同坐標(biāo))再次用高斯計算能量，得到MO軌道的fch文件，同樣做電子密度格點。兩個格點文件的質(zhì)量選擇相同，二者相減會得到差值，也就是說.40和MO電子密度格點不一樣。這是正常的嗎？ |
回復(fù)此樓» 猜你喜歡
一志愿南大理論與計算化學(xué)方向（070300）342分，求調(diào)劑
已經(jīng)有1人回復(fù)
-
歡迎化學(xué)、化工、物理、材料、或計算機背景的考研咨詢
已經(jīng)有0人回復(fù)
-
物理化學(xué)論文潤色/翻譯怎么收費?
已經(jīng)有120人回復(fù)
-
PDF在線壓縮軟件
已經(jīng)有5人回復(fù)
-
材料/化學(xué)相關(guān)專業(yè)2026級學(xué)術(shù)/專業(yè)型碩士研究生
已經(jīng)有19人回復(fù)
-
求催化作用基礎(chǔ)第三版 （甄開吉) 課件
已經(jīng)有3人回復(fù)
-
26年材料/化學(xué)相關(guān)專業(yè)博士研究生招生 溫州大學(xué)化學(xué)與材料工程學(xué)院(26年秋季入學(xué))
已經(jīng)有0人回復(fù)
-
求助電催化劑制備，怎么在碳布表面生長PdCo合金
已經(jīng)有0人回復(fù)
-
求助電催化劑制備，怎么在碳布表面生長PdCo合金
已經(jīng)有0人回復(fù)
-
304求調(diào)劑
已經(jīng)有3人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 建議multiwfn開發(fā)GUI界面 文件打開什么的效率太低了
已經(jīng)有10人回復(fù)
- 如何用Multiwfn軟件做出P DOS?
已經(jīng)有4人回復(fù)
- 使用Multiwfn繪制原子軌道圖形、研究原子殼層結(jié)構(gòu)及相對論效應(yīng)的影響
已經(jīng)有17人回復(fù)
- Multiwfn的初步學(xué)習(xí)，越詳細(xì)越好。
已經(jīng)有6人回復(fù)
- MultiWFN想要生成用于NBO分析的cub時不能載入.37文件
已經(jīng)有3人回復(fù)
- 使用Multiwfn繪制NBO及相關(guān)軌道
已經(jīng)有30人回復(fù)
- 用Multiwfn、Gaussview、Molekel、VMD觀看龍蝦、盆景、骨盆、大腦
已經(jīng)有16人回復(fù)
- 使用Multiwfn做拓?fù)浞治鲆约坝嬎愎聦﹄娮咏嵌?/a>
已經(jīng)有18人回復(fù)
- 回答一些關(guān)于Multiwfn的疑問以及未來Multiwfn的發(fā)展打算
已經(jīng)有37人回復(fù)
sobereva
至尊木蟲 (著名寫手)
本人已永久離開小木蟲
- QC強帖: 17
- 應(yīng)助: 226 (大學(xué)生)
- 貴賓: 5.4
- 金幣: 2180.5
- 紅花: 217
- 帖子: 1241
- 在線: 366.8小時
- 蟲號: 281271
- 注冊: 2006-09-23
- 專業(yè): 宇宙學(xué)
【答案】應(yīng)助回帖
★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
gmy1990: 金幣+3 2013-06-13 21:05:02
klaus1987: 金幣+5, ★★★很有幫助 2013-06-14 09:54:20
感謝參與，應(yīng)助指數(shù) +1
gmy1990: 金幣+3 2013-06-13 21:05:02
klaus1987: 金幣+5, ★★★很有幫助 2013-06-14 09:54:20
|
.40是MO，.39才是NLMO（手冊寫錯了） NLMO和MO截然不同，前者可認(rèn)為是對后者酉變換產(chǎn)生的定域化軌道 .40和fch里的MO應(yīng)當(dāng)是一樣的，存在差值可能是微小的數(shù)值精度問題所致，或者什么其它地方?jīng)]有完全對應(yīng)上。最好直接通過在主功能0觀看等值面來比較二者。 |
2樓2013-06-13 18:10:14
3樓2013-06-14 09:59:40
4樓2013-06-14 10:49:50
5樓2013-06-14 10:51:59
sobereva
至尊木蟲 (著名寫手)
本人已永久離開小木蟲
- QC強帖: 17
- 應(yīng)助: 226 (大學(xué)生)
- 貴賓: 5.4
- 金幣: 2180.5
- 紅花: 217
- 帖子: 1241
- 在線: 366.8小時
- 蟲號: 281271
- 注冊: 2006-09-23
- 專業(yè): 宇宙學(xué)
【答案】應(yīng)助回帖
★ ★ ★ ★ ★ ★ ★
gmy1990: 金幣+2 2013-06-14 16:53:14
klaus1987: 金幣+5, ★★★★★最佳答案 2013-06-16 15:28:25
gmy1990: 金幣+2 2013-06-14 16:53:14
klaus1987: 金幣+5, ★★★★★最佳答案 2013-06-16 15:28:25
|
你可以去比較.40和fch當(dāng)中特定的分子軌道的密度，但是不要去比較它們計算的總密度。 用.40作為輸入后，選主功能6再選3，可以看到所有軌道都是占據(jù)的，也就是說，以.40為輸入時無法正確判斷軌道占據(jù)數(shù)，因此計算的總密度毫無意義。這是因為.40里面本身就沒有提供軌道占據(jù)信息。計算分子性質(zhì)的時候，不要用NBO plot文件，而應(yīng)當(dāng)直接用fch/wfn/wfx/molden作為輸入。 只有以.31(NLMO)作為輸入時，計算出的分子性質(zhì)才和實際一致，這是手冊中這句話的含義"Notice that among all types of the orbitals generated by NBO program, only using NLMO to calculate real space functions is meaningful!” |
6樓2013-06-14 14:11:24
7樓2013-06-16 15:29:56
8樓2013-06-16 17:18:12
sobereva
至尊木蟲 (著名寫手)
本人已永久離開小木蟲
- QC強帖: 17
- 應(yīng)助: 226 (大學(xué)生)
- 貴賓: 5.4
- 金幣: 2180.5
- 紅花: 217
- 帖子: 1241
- 在線: 366.8小時
- 蟲號: 281271
- 注冊: 2006-09-23
- 專業(yè): 宇宙學(xué)
★ ★ ★
gkf高: 金幣+3 2013-06-17 11:17:08
gkf高: 金幣+3 2013-06-17 11:17:08
|
這些.xx文件之間也不能相互比較計算出來的總的密度或其它實空間函數(shù) 原理上說，PNAO/NAO/PNHO/NHO/PNBO/NBO/PNLMO/NLMO/MO這些不管用哪種類型軌道，實空間函數(shù)（如總密度）的計算結(jié)果都是完全一致的，所以實空間函數(shù)的計算只要用MO就行了（.fch/.molden/.wfx/.wfn作為輸入文件），基于其它的軌道來計算由于結(jié)果和基于MO的完全一致，所以沒必要用它們。 而之所以我提到不要用NBO程序給出的軌道來計算（NLMO除外），一方面是因為這些NBO plot文件里往往沒有正確記錄或者根本沒記錄軌道占據(jù)數(shù)；另一方面是由于以某些類型軌道為基時（如NBO），整體的密度矩陣不是對角矩陣，此時以Multiwfn目前基于MO計算密度的方式來計算NBO表達(dá)的密度結(jié)果會和實際有一定偏差。 |
9樓2013-06-16 18:57:29
sobereva
至尊木蟲 (著名寫手)
本人已永久離開小木蟲
- QC強帖: 17
- 應(yīng)助: 226 (大學(xué)生)
- 貴賓: 5.4
- 金幣: 2180.5
- 紅花: 217
- 帖子: 1241
- 在線: 366.8小時
- 蟲號: 281271
- 注冊: 2006-09-23
- 專業(yè): 宇宙學(xué)
★ ★ ★
gkf高: 金幣+3 2013-06-17 11:17:18
gkf高: 金幣+3 2013-06-17 11:17:18
|
以下是我開發(fā)Multiwfn時針對這個問題寫的筆記，涉及到一些細(xì)節(jié)，可以參考 “在計算密度等屬性時，按照現(xiàn)有方法，將NBO占據(jù)數(shù)乘上NBO軌道波函數(shù)平方，再加和，得到的并不是真正的總密度。因為以NBO為基時密度矩陣并不是對角的，因為獲得NBO時并不是直接對整體密度矩陣對角化，而是依次將各個塊密度矩陣對角化。為了得到正確的密度，必須輸出以NBO為基的密度矩陣，將每個矩陣元乘上對應(yīng)的兩個NBO在此處的波函數(shù)值，然后整體加和，也就是必須把非對角元考慮進去。對于NAO也是如此。不過，按Multiwfn現(xiàn)有方法計算，對于NBO、NAO、NHO，密度積分倒是和期望的一致。因為它們都是相互正交的。使用NBO算的實空間函數(shù)和真實值偏差不大。 NLMO和普通方式獲得的定域化軌道一樣，以其為基時整體密度矩陣也是對角化的，因此可以直接將密度乘上波函數(shù)平方加和得到真實密度。正如NBCP文章里式子所述。 至于前頭帶P的就更不能這么算了，連正交都不滿足，因此積分值都不對�！� |
» 本帖已獲得的紅花（最新10朵）
klaus198710樓2013-06-16 19:01:50
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 klaus1987 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[基金申請] 被我言中：新模板不強調(diào)格式了，假專家開始管格式了 +4 | beefly 2026-03-14 | 4/200 |
|
|---|---|---|---|---|
|
|
[考研] 材料工程專碩調(diào)劑 +5 | 204818@lcx 2026-03-17 | 5/250 |
|
|
|
[碩博家園] 湖北工業(yè)大學(xué) 生命科學(xué)與健康學(xué)院-課題組招收2026級食品/生物方向碩士 +3 | 1喜春8 2026-03-17 | 5/250 |
|
|
[考研] 308求調(diào)劑 +4 | 是Lupa啊 2026-03-16 | 4/200 |
|
|
|
[考研] 085601求調(diào)劑 +4 | Du.11 2026-03-16 | 4/200 |
|
|
|
[考研] 求調(diào)劑，總分315，考的生物醫(yī)藥，一志愿湖南師范大學(xué)。調(diào)劑到任何專業(yè)都可以 +4 | 小丁想進步 2026-03-11 | 5/250 |
|
|
|
[考研] 289求調(diào)劑 +6 | 步川酷紫123 2026-03-11 | 6/300 |
|
|
|
[考研] 267一志愿南京工業(yè)大學(xué)0817化工求調(diào)劑 +6 | SUICHILD 2026-03-12 | 6/300 |
|
|
|
[考研] 藥學(xué)383 求調(diào)劑 +3 | 藥學(xué)chy 2026-03-15 | 4/200 |
|
|
|
[考研] 085600調(diào)劑 +5 | 漾漾123sun 2026-03-12 | 6/300 |
|
|
|
[考研] 085600材料與化工 求調(diào)劑 +13 | enenenhui 2026-03-13 | 14/700 |
|
|
|
[教師之家] 焦慮 +7 | 水冰月月野兔 2026-03-13 | 9/450 |
|
|
|
[考研] 326求調(diào)劑 +3 | mlpqaz03 2026-03-15 | 3/150 |
|
|
|
[考博] 東華理工大學(xué)化材專業(yè)26屆碩士博士申請 +6 | zlingli 2026-03-13 | 6/300 |
|
|
|
[考研] 341求調(diào)劑 +3 | 番茄頭--- 2026-03-10 | 3/150 |
|
|
|
[考研] 279求調(diào)劑 +3 | Dizzy123@ 2026-03-10 | 3/150 |
|
|
|
[考研] 304求調(diào)劑 +6 | Mochaaaa 2026-03-12 | 7/350 |
|
|
|
[考研] 材料與化工085600調(diào)劑求老師收留 +9 | jiaanl 2026-03-11 | 9/450 |
|
|
|
[考研] 085600材料與化工 309分請求調(diào)劑 +7 | dtdxzxx 2026-03-12 | 8/400 |
|
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +3 | 崔wj 2026-03-12 | 4/200 |
|