 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 4 | 1/1 | 返回列表 |
| 查看: 1547 | 回復(fù): 3 | |||
muxiachuixue木蟲 (正式寫手)
|
[求助]
NBO 電荷總數(shù)不為零
|
|
小弟最近算一個(gè)63個(gè)原子的分子，有C/H/O/N/S幾種原子，分子未經(jīng)高斯優(yōu)化，用b3lyp/6-311++g(d,p)算nbo，電荷總數(shù)為-1.8，可分子明明是中性的啊，輸入文件也規(guī)定了charge=0；再將分子用b3lyp/6-31+g(d,p)優(yōu)化，b3lyp/6-311++g(d,p)算nbo，電荷總數(shù)差不多為0 故請(qǐng)教大俠，為何未優(yōu)化之前算出的nbo電荷總和不為0呢？ 以下貼出輸入文件和輸出文件 Summary of Natural Population Analysis: Natural Population Natural ----------------------------------------------- Atom No Charge Core Valence Rydberg Total ----------------------------------------------------------------------- C 1 0.05001 1.99919 3.92264 0.02816 5.94999 N 2 -0.45280 1.99943 5.43155 0.02182 7.45280 C 3 0.09288 1.99912 3.88022 0.02778 5.90712 C 4 -0.10952 1.99892 4.09574 0.01487 6.10952 C 5 -0.06187 1.99886 4.04467 0.01833 6.06187 C 6 -0.18951 1.99908 4.17558 0.01485 6.18951 C 7 -0.12920 1.99910 4.11328 0.01681 6.12920 C 8 -0.20789 1.99915 4.19040 0.01835 6.20789 C 9 -0.17112 1.99906 4.15306 0.01900 6.17112 C 10 -0.25147 1.99862 4.22640 0.02645 6.25147 S 11 2.16679 8.10518 3.66823 2.05979 13.83321 N 12 -0.93644 1.99938 5.90754 0.02952 7.93644 C 13 -0.02604 1.99921 4.00486 0.02196 6.02604 C 14 -0.58809 1.99930 4.57808 0.01071 6.58809 O 15 -0.91952 1.99980 6.90725 0.01247 8.91952 O 16 -1.38312 1.99980 5.57151 1.81181 9.38312 C 17 -0.17626 1.99919 4.15813 0.01894 6.17626 N 18 -0.20247 1.99919 5.17253 0.03075 7.20247 N 19 -0.11881 1.99955 5.08428 0.03498 7.11881 N 20 -0.25648 1.99955 5.23112 0.02580 7.25648 C 21 0.08265 1.99898 3.88401 0.03436 5.91735 C 22 -0.05744 1.99875 4.03747 0.02122 6.05744 C 23 -0.20684 1.99916 4.16820 0.03948 6.20684 N 24 -0.81971 1.99946 5.54323 0.27703 7.81971 C 25 -0.58863 1.99919 3.56164 1.02779 6.58863 C 26 -0.01296 1.99880 3.86661 0.14755 6.01296 C 27 -0.09641 1.99886 4.06605 0.03149 6.09641 C 28 -0.39495 1.99920 4.37714 0.01861 6.39495 C 29 -0.21182 1.99923 4.15193 0.06066 6.21182 C 30 -0.06698 1.99875 4.04751 0.02072 6.06698 C 31 0.12853 1.99892 3.84499 0.02757 5.87147 N 32 -0.56214 1.99905 5.50775 0.05534 7.56214 C 33 -0.19784 1.99911 4.18066 0.01806 6.19784 C 34 -0.22406 1.99918 4.20635 0.01853 6.22406 C 35 -0.21721 1.99919 4.19905 0.01897 6.21721 C 36 -0.22166 1.99909 4.20590 0.01667 6.22166 H 37 0.19605 0.00000 0.80172 0.00224 0.80395 H 38 0.19762 0.00000 0.80102 0.00136 0.80238 H 39 0.23123 0.00000 0.76619 0.00258 0.76877 H 40 0.20967 0.00000 0.78812 0.00221 0.79033 H 41 0.21467 0.00000 0.78330 0.00203 0.78533 H 42 0.23104 0.00000 0.76526 0.00370 0.76896 H 43 0.43595 0.00000 0.56060 0.00345 0.56405 H 44 0.21792 0.00000 0.77990 0.00218 0.78208 H 45 0.20541 0.00000 0.79297 0.00163 0.79459 H 46 0.19461 0.00000 0.80352 0.00187 0.80539 H 47 0.22500 0.00000 0.77299 0.00201 0.77500 H 48 0.22843 0.00000 0.76997 0.00160 0.77157 H 49 0.23057 0.00000 0.76731 0.00212 0.76943 H 50 0.22866 0.00000 0.76901 0.00233 0.77134 H 51 0.19969 0.00000 0.79020 0.01011 0.80031 H 52 0.17823 0.00000 0.81560 0.00617 0.82177 H 53 0.01128 0.00000 0.92878 0.05994 0.98872 H 54 0.11182 0.00000 0.82727 0.06091 0.88818 H 55 0.20988 0.00000 0.78604 0.00407 0.79012 H 56 0.20393 0.00000 0.79256 0.00351 0.79607 H 57 0.22122 0.00000 0.77247 0.00631 0.77878 H 58 0.14636 0.00000 0.84176 0.01189 0.85364 H 59 0.39280 0.00000 0.60310 0.00410 0.60720 H 60 0.20923 0.00000 0.78918 0.00159 0.79077 H 61 0.20118 0.00000 0.79693 0.00189 0.79882 H 62 0.20196 0.00000 0.79625 0.00179 0.79804 H 63 0.19995 0.00000 0.79774 0.00231 0.80005 ======================================================================= * Total * -1.80404 78.07562 181.42533 6.30308 265.80404 |
回復(fù)此樓» 本帖附件資源列表
-
歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺(tái)，不對(duì)該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識(shí)產(chǎn)權(quán)問題，其責(zé)任在于用戶本人，如對(duì)版權(quán)有異議，請(qǐng)聯(lián)系郵箱：xiaomuchong@tal.com - 附件 1 : A1B7_nbo .gjf
- 附件 2 : A1B7_nbo.out
2012-02-17 10:19:02, 4.56 K
2012-02-17 10:19:04, 230.39 K
» 收錄本帖的淘帖專輯推薦
 模擬，理論 |
» 猜你喜歡
PbS量子點(diǎn)紫外吸收
已經(jīng)有0人回復(fù)
-
比利時(shí)魯汶大學(xué)與國(guó)家留學(xué)基金委共同資助博士研究生CSC-KU Leuven PhD Scholarship
已經(jīng)有0人回復(fù)
-
物理化學(xué)論文潤(rùn)色/翻譯怎么收費(fèi)?
已經(jīng)有142人回復(fù)
-
[電子教材]天津大學(xué)物理化學(xué)（第七版）上、下冊(cè)
已經(jīng)有3人回復(fù)
-
陰離子交換膜電解二氧化碳還原
已經(jīng)有0人回復(fù)
-
倫敦大學(xué)學(xué)院Benjamin Moss博士招收人工智能、原位光譜、催化相關(guān)26年博士生
已經(jīng)有11人回復(fù)
-
中國(guó)科學(xué)院成都山地災(zāi)害與環(huán)境研究所2026年特別研究助理（博士后）招聘啟事
已經(jīng)有0人回復(fù)
-
九江學(xué)院2026年最新高層次人才招聘公告
已經(jīng)有0人回復(fù)
-
深圳信息職業(yè)技術(shù)學(xué)院2026年教師招聘公告（最新）
已經(jīng)有0人回復(fù)
-
研究生招生
已經(jīng)有0人回復(fù)
-
山西大學(xué)資源與環(huán)境工程研究所化工分離/環(huán)境功能材料小組招收2026年博士生
已經(jīng)有1人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- NBO正常結(jié)束但是沒有分析結(jié)果
已經(jīng)有11人回復(fù)
- 蛋白質(zhì)的電荷不為整數(shù)問題
已經(jīng)有3人回復(fù)
- 2-二硝基甲苯的4位和6位碳的NBO電荷
已經(jīng)有3人回復(fù)
- 用gaussian計(jì)算負(fù)一價(jià)離子的NBO時(shí)，輸出文件中電荷顯示為何為0
已經(jīng)有7人回復(fù)
- 高斯計(jì)算，用不同基組算出的NBO軌道相位不同
已經(jīng)有3人回復(fù)
- 關(guān)于NBO及Mulliken電荷分析的疑問
已經(jīng)有10人回復(fù)
- NBO電荷比較親電活性大小
已經(jīng)有8人回復(fù)
- MultiWFN想要生成用于NBO分析的cub時(shí)不能載入.37文件
已經(jīng)有3人回復(fù)
- 流式細(xì)胞儀測(cè)蛋白熒光能不能樣品處理好后放過夜第二天再測(cè)
已經(jīng)有3人回復(fù)
- 求教，用簡(jiǎn)便的化學(xué)方法鑒別下列幾種酸：草酸、丙酮酸、丙酸
已經(jīng)有4人回復(fù)
- 急求幫助，請(qǐng)各位幫忙！謝謝計(jì)算NBO電荷分析出現(xiàn)如下問題
已經(jīng)有6人回復(fù)
- NBO 電荷分析中電荷總和不為零 分析 處理
已經(jīng)有6人回復(fù)
- Mulliken 布局分析任何分析電荷分布、轉(zhuǎn)移和化學(xué)性質(zhì)
已經(jīng)有10人回復(fù)
- 【求助】NBO軌道指的是nbo電荷的分布情況嗎？
已經(jīng)有3人回復(fù)
- 電荷密度圖 等高線不對(duì)稱
已經(jīng)有13人回復(fù)
- 大家做NBO分析都用其中的什么內(nèi)容呢？總覺得自己分析不夠透徹...
已經(jīng)有21人回復(fù)
- VESTA為什么不能導(dǎo)出 2D 電荷密度的數(shù)據(jù)？
已經(jīng)有5人回復(fù)
- 【求助】密立根電荷和NBO電荷的差異
已經(jīng)有4人回復(fù)
- 【求助】請(qǐng)問NPA電荷和NBO電荷有什么區(qū)別？
已經(jīng)有7人回復(fù)
- 【求助】Mulliken電荷分布總和不為0是什么原因
已經(jīng)有7人回復(fù)
quantumor
金蟲 (著名寫手)
快樂兔子
- 應(yīng)助: 13 (小學(xué)生)
- 金幣: 8506.3
- 散金: 10
- 紅花: 14
- 帖子: 1898
- 在線: 215.9小時(shí)
- 蟲號(hào): 630331
- 注冊(cè): 2008-10-19
- 性別: GG
- 專業(yè): 理論和計(jì)算化學(xué)
- 管轄: 量子化學(xué)
2樓2012-02-17 13:55:17
muxiachuixue
木蟲 (正式寫手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 2068.6
- 散金: 2
- 紅花: 5
- 帖子: 315
- 在線: 178.7小時(shí)
- 蟲號(hào): 1285058
- 注冊(cè): 2011-05-04
- 性別: GG
- 專業(yè): 合成藥物化學(xué)
|
你好，我發(fā)郵件問了高斯官方，高斯官方給的建議是加上IOp(6/76=8)來提高線性相關(guān)的標(biāo)準(zhǔn)，經(jīng)測(cè)試，電荷總和已正常，以下是gaussian的回復(fù)： If you look at the beginning of the NBO analysis you will see the problem, Remove ICRem= 289 ICRemI= 290 L=0 PNAO number 289 atom 14 type 51 NPNAO= 520 due to symmetric orthogonalization 2, value= 0.00D+00. Computed number of electrons 148.18290763 does not match expected value of 148. where the program was having some problems orthogonalizing the basis set. We have seen this with some other calculations with diffuse basis sets and sometimes tightening the criteria for linear dependence will help. Using # rb3lyp/6-311++g(d,p) pop=npa IOp(6/76=8) which changes it from 10^-6 to 10^-8 will still be stable and yet correct the removal of important contributions. |
3樓2012-03-08 18:46:10
quantumor
金蟲 (著名寫手)
快樂兔子
- 應(yīng)助: 13 (小學(xué)生)
- 金幣: 8506.3
- 散金: 10
- 紅花: 14
- 帖子: 1898
- 在線: 215.9小時(shí)
- 蟲號(hào): 630331
- 注冊(cè): 2008-10-19
- 性別: GG
- 專業(yè): 理論和計(jì)算化學(xué)
- 管轄: 量子化學(xué)
4樓2012-03-09 09:37:09
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 muxiachuixue 的主題更新
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研]
 271材料工程求調(diào)劑
+3
|
.6lL 2026-03-18 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 331求調(diào)劑（0703有機(jī)化學(xué) +7 | ZY-05 2026-03-13 | 8/400 |
|
|
|
[考博] 環(huán)境領(lǐng)域全國(guó)重點(diǎn)實(shí)驗(yàn)室招收博士1-2名 +3 | QGZDSYS 2026-03-13 | 5/250 |
|
|
|
[考研] 工科材料085601 279求調(diào)劑 +6 | 困于星晨 2026-03-17 | 6/300 |
|
|
|
[考研] 307求調(diào)劑 +3 | 冷笙123 2026-03-17 | 3/150 |
|
|
|
[考研] 環(huán)境工程調(diào)劑 +8 | 大可digkids 2026-03-16 | 8/400 |
|
|
|
[考研] 085601求調(diào)劑 +4 | Du.11 2026-03-16 | 4/200 |
|
|
|
[考研] 【0856】化學(xué)工程（085602）313 分，本科學(xué)科評(píng)估A類院�；瘜W(xué)工程與工藝，誠(chéng)求調(diào)劑 +7 | 小劉快快上岸 2026-03-11 | 8/400 |
|
|
|
[考研] 有沒有道鐵/土木的想調(diào)劑南林，給自己招師弟中～ +3 | TqlXswl 2026-03-16 | 7/350 |
|
|
|
[考研] 梁成偉老師課題組歡迎你的加入 +8 | 一鴨鴨喲 2026-03-14 | 10/500 |
|
|
|
[基金申請(qǐng)] 國(guó)自科面上基金字體 +6 | iwuli 2026-03-12 | 7/350 |
|
|
|
[考研] 333求調(diào)劑 +3 | 文思客 2026-03-16 | 7/350 |
|
|
|
[考研] 326求調(diào)劑 +4 | 諾貝爾化學(xué)獎(jiǎng)覬?/a> 2026-03-15 | 7/350 |
|
|
|
[考研] 304求調(diào)劑 +4 | ahbd 2026-03-14 | 4/200 |
|
|
|
[考研] 中科大材料與化工319求調(diào)劑 +3 | 孟鑫材料 2026-03-14 | 3/150 |
|
|
|
[基金申請(qǐng)] 現(xiàn)在如何回避去年的某一個(gè)專家，不知道名字 +3 | zk200107 2026-03-12 | 6/300 |
|
|
|
[考研] 255求調(diào)劑 +3 | 李嘉慧， 2026-03-12 | 4/200 |
|
|
|
[考研] 330求調(diào)劑 +3 | ?醬給調(diào)劑跪了 2026-03-13 | 3/150 |
|
|
|
[考研] 一志愿山大07化學(xué) 332分 四六級(jí)已過 本科山東雙非 求調(diào)劑！ +3 | 不想理你 2026-03-12 | 3/150 |
|
|
|
[考研] 070303一志愿西北大學(xué)學(xué)碩310找調(diào)劑 +3 | d如愿上岸 2026-03-13 | 3/150 |
|