[求助] gaussian09計(jì)算1-丁基-3-甲基咪唑陽(yáng)離子電荷布居的問(wèn)題

大家好，本人是量化新手，最近學(xué)習(xí)用gaussian09進(jìn)行計(jì)算，在計(jì)算1-丁基-3-甲基咪唑陽(yáng)離子的電荷布居時(shí)遇到了一些問(wèn)題，想向大家請(qǐng)教一下。
我用了兩種方法，第一種是NBO，輸入文件和結(jié)果如下：
%nprocshared=8
%chk=bmim.chk
%mem=1000MW
# opt=tight b3lyp/6-311++g(d,p) pop=nbo sp test

bmim

1 1
C                2.36134200 -0.30291700 0.56377200
C                1.10194200 -0.62271700 0.14137200
N                0.33804200 0.51268300 0.30847200
C                1.07744200 1.50548300 0.81267200
N                2.31024200 1.01358300 0.96987200
C                3.47578200 1.78626400 1.51253000
C                -1.09057000 0.61694700 -0.02182900
C                -1.58383800 -0.71504400 -0.57138800
C                -3.06395700 -0.60702100 -0.91359800
C                -3.55722600 -1.93901200 -1.46315700
H                3.24051300 -0.96388600 0.57715300
H                0.76480600 -1.59256200 -0.25326800
H                0.73875000 2.52458100 1.05084700
H                4.37683100 1.13325600 1.53399800
H                3.67201000 2.66702900 0.86098100
H                3.24345300 2.13166700 2.54475600
H                -1.66513700 0.87466900 0.89589700
H                -1.23711500 1.41118600 -0.78764700
H                -1.00927100 -0.97276700 -1.48911400
H                -1.43729300 -1.50928400 0.19443000
H                -3.63852400 -0.34929800 0.00412800
H                -3.21050300 0.18721900 -1.67941500
H                -4.63888900 -1.86006900 -1.71324200
H                -2.98265900 -2.19673500 -2.38088300
H                -3.41068000 -2.73325200 -0.69733900

1 2 2.0 5 1.5 11 1.0
2 3 1.5 12 1.0
3 4 1.5 7 1.0
4 5 1.5 13 1.0
5 6 1.0
6 14 1.0 15 1.0 16 1.0
7 8 1.0 17 1.0 18 1.0
8 9 1.0 19 1.0 20 1.0
9 10 1.0 21 1.0 22 1.0
10 23 1.0 24 1.0 25 1.0
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
——————————————————下面是結(jié)果
                                    Natural Population
            Natural  -----------------------------------------------
Atom  No Charge       Core    Valence Rydberg    Total
-----------------------------------------------------------------------
   C 1 -0.15060    1.99911    4.13301 0.01848    6.15060
   C 2 -0.15109    1.99913    4.13376 0.01820    6.15109
   N 3 -0.25333    1.99926    5.23878 0.01529    7.25333
   C 4 0.33509    1.99915    3.64526 0.02050    5.66491
   N 5 -0.24413    1.99926    5.23166 0.01322    7.24413
   C 6 -0.60923    1.99924    4.59923 0.01076    6.60923
   C 7 -0.30384    1.99919    4.28929 0.01536    6.30384
   C 8 -0.64776    1.99923    4.63549 0.01304    6.64776
   C 9 -0.60424    1.99925    4.59313 0.01185    6.60424
   C 10 -0.92808    1.99929    4.91998 0.00881    6.92808
   H 11 0.31324    0.00000    0.68532 0.00143    0.68676
   H 12 0.31641    0.00000    0.68199 0.00160    0.68359
   H 13 0.29205    0.00000    0.70654 0.00141    0.70795
   H 14 0.29814    0.00000    0.70066 0.00120    0.70186
   H 15 0.29119    0.00000    0.70754 0.00127    0.70881
   H 16 0.29121    0.00000    0.70752 0.00127    0.70879
   H 17 0.29077    0.00000    0.70733 0.00190    0.70923
   H 18 0.29082    0.00000    0.70729 0.00190    0.70918
   H 19 0.31321    0.00000    0.68473 0.00206    0.68679
   H 20 0.31326    0.00000    0.68468 0.00206    0.68674
   H 21 0.30551    0.00000    0.69256 0.00194    0.69449
   H 22 0.30556    0.00000    0.69251 0.00193    0.69444
   H 23 0.31523    0.00000    0.68345 0.00132    0.68477
   H 24 0.31029    0.00000    0.68817 0.00155    0.68971
   H 25 0.31032    0.00000    0.68813 0.00155    0.68968
————————————————————
第二種是Mulliken，輸入文件沒(méi)有pop=nbo，其他相同，結(jié)果是：
Mulliken atomic charges:
            1
   1  C -0.373840
   2  C 0.022657
   3  N 0.139338
   4  C 0.094221
   5  N 0.136618
   6  C -0.507527
   7  C -0.530899
   8  C -0.208323
   9  C -0.777532
10  C -1.008057
11  H 0.263268
12  H 0.261924
13  H 0.232876
14  H 0.249825
15  H 0.266894
16  H 0.266870
17  H 0.289971
18  H 0.289918
19  H 0.283122
20  H 0.283093
21  H 0.262814
22  H 0.262852
23  H 0.279577
24  H 0.260174
25  H 0.260165
我的問(wèn)題是
1）輸入文件中的sp test是什么意思？我在gaussian的文檔中沒(méi)有找到test關(guān)鍵詞。（輸入文件是我用chem3d保存構(gòu)型為gjf后，用gaussview載入的。setup選的是opt，tight，b3lyp和6-311++g(d,p)。gaussview自動(dòng)在下邊添加了Additional Keywords：sp test。我試過(guò)如果沒(méi)有這兩個(gè)關(guān)鍵詞，程序無(wú)法運(yùn)行）；
2）兩種方法算出來(lái)的電荷布居為何與直觀感覺(jué)相差很大？應(yīng)該是咪唑環(huán)上正點(diǎn)強(qiáng)，氫上接近中性才比較符合直覺(jué)啊。而計(jì)算出來(lái)咪唑環(huán)是帶負(fù)電的。相比之下，我用Material Studio計(jì)算的Mulliken布居比較“正�！保ㄍ瑸閎3lyp）。是不是我基組和方法選擇的問(wèn)題呢？
謝謝大家！

回復(fù)此樓

» 猜你喜歡

細(xì)胞培養(yǎng)，這22個(gè)細(xì)節(jié)一定要注意�。ㄒ唬� 已經(jīng)有0人回復(fù)
CSC訪學(xué)博后項(xiàng)目獲批，外方學(xué)校暫停合作該怎么辦？已經(jīng)有51人回復(fù)
物理化學(xué)論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有134人回復(fù)
大灣區(qū)大學(xué)劉天輝課題組招收2026級(jí)博士生（與中山大學(xué)物理學(xué)院聯(lián)合培養(yǎng)）已經(jīng)有22人回復(fù)
PbS量子點(diǎn)紫外吸收已經(jīng)有0人回復(fù)
比利時(shí)魯汶大學(xué)與國(guó)家留學(xué)基金委共同資助博士研究生CSC-KU Leuven PhD Scholarship 已經(jīng)有0人回復(fù)
半透性容器進(jìn)行高濕或者低濕研究已經(jīng)有0人回復(fù)
倫敦大學(xué)學(xué)院Benjamin Moss博士招收人工智能、原位光譜、催化相關(guān)26年博士生已經(jīng)有11人回復(fù)
分子動(dòng)力學(xué)模擬合作需求已經(jīng)有2人回復(fù)
中國(guó)科學(xué)院成都山地災(zāi)害與環(huán)境研究所2026年特別研究助理（博士后）招聘啟事已經(jīng)有0人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

怎么計(jì)算 atomic spin density 已經(jīng)有4人回復(fù)
請(qǐng)問(wèn)Gaussian計(jì)算frequency出現(xiàn)下面的問(wèn)題，如何解決？已經(jīng)有11人回復(fù)
電子軌道能量用什么方法計(jì)算比較好？已經(jīng)有6人回復(fù)
【小卒個(gè)人文集】寫(xiě)給初學(xué)者:小議形式電荷對(duì)量化初學(xué)者的誤導(dǎo) 已經(jīng)有27人回復(fù)
【求助】離子鍵和共價(jià)鍵的區(qū)分已經(jīng)有17人回復(fù)

1樓 2015-06-10 16:27:58

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

凡士林護(hù)手霜

銅蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 1 (幼兒園)
金幣: 1077.7
散金: 40
紅花: 1
帖子: 421
在線: 250.2小時(shí)
蟲(chóng)號(hào): 2911292
注冊(cè): 2014-01-03
性別: GG
專(zhuān)業(yè): 污染控制化學(xué)

你好，請(qǐng)問(wèn)你是做量子化學(xué)計(jì)算的研究嗎，你們單位購(gòu)買(mǎi)了嗎

贊一下

回復(fù)此樓

夢(mèng)想在心中，路在腳下。

2樓2016-04-26 13:49:43

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 354求調(diào)劑 +4	Tyoumou 2026-03-18	7/350	2026-03-18 21:45 by Tyoumou
[考研] 321求調(diào)劑 +3	何潤(rùn)采123 2026-03-18	3/150	2026-03-18 21:27 by li123456789.
[考研] 一志愿天津大學(xué)化學(xué)工藝專(zhuān)業(yè)（081702）315分求調(diào)劑 +10	yangfz 2026-03-17	10/500	2026-03-18 20:14 by walc
[考研] 344求調(diào)劑 +6	knight344 2026-03-16	7/350	2026-03-18 20:13 by walc
[考研] 307求調(diào)劑 +3	冷笙123 2026-03-17	3/150	2026-03-18 09:55 by macy2011
[考研] 0703化學(xué)336分求調(diào)劑 +6	zbzihdhd 2026-03-15	7/350	2026-03-18 09:53 by zhukairuo
[考研] 326求調(diào)劑 +5	上岸的小葡 2026-03-15	6/300	2026-03-17 17:26 by ruiyingmiao
[考研] 材料專(zhuān)碩326求調(diào)劑 +6	墨煜姒莘 2026-03-15	7/350	2026-03-17 17:10 by ruiyingmiao
[考研] 302求調(diào)劑 +4	小賈同學(xué)123 2026-03-15	8/400	2026-03-17 10:33 by 小賈同學(xué)123
[考研] 東南大學(xué)364求調(diào)劑 +5	JasonYuiui 2026-03-15	5/250	2026-03-16 21:28 by 木瓜膏
[考研] 321求調(diào)劑 +5	大米飯！ 2026-03-15	5/250	2026-03-16 16:33 by houyaoxu
[考研] 070300化學(xué)學(xué)碩求調(diào)劑 +6	太想進(jìn)步了0608 2026-03-16	6/300	2026-03-16 16:13 by kykm678
[考研] 285求調(diào)劑 +6	ytter 2026-03-12	6/300	2026-03-16 15:05 by njzyff
[考研] 277材料科學(xué)與工程080500求調(diào)劑 +3	自由煎餅果子 2026-03-16	3/150	2026-03-16 14:10 by 運(yùn)氣yunqi
[考研] 0856專(zhuān)碩279求調(diào)劑 +5	加油加油！? 2026-03-15	5/250	2026-03-15 11:58 by 2020015
[碩博家園] 085600 260分求調(diào)劑 +3	天空還下雨么 2026-03-13	5/250	2026-03-13 18:46 by 天空還下雨么
[考研] 求調(diào)劑 +5	一定有學(xué)上- 2026-03-12	5/250	2026-03-13 18:31 by ms629
[考研] 290求調(diào)劑 +7	ADT 2026-03-12	7/350	2026-03-13 15:17 by JourneyLucky
[考研] 085600材料與化工 309分請(qǐng)求調(diào)劑 +7	dtdxzxx 2026-03-12	8/400	2026-03-13 14:43 by jxchenghu
[考博] 福州大學(xué)楊黃浩課題組招收2026年專(zhuān)業(yè)學(xué)位博士研究生，2026.03.20截止 +3	Xiangyu_ou 2026-03-12	3/150	2026-03-13 09:36 by duanwu655