版塊導航: 正在加載中...

登錄注冊

應《網(wǎng)絡安全法》要求，自2017年10月1日起，未進行實名認證將不得使用互聯(lián)網(wǎng)跟帖服務。為保障您的帳號能夠正常使用，請盡快對帳號進行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

返回列表

foclnbrsic

銅蟲 (初入文壇)

應助: 0 (幼兒園)
金幣: 61.5
帖子: 4
在線: 11.9小時
蟲號: 1758794
注冊: 2012-04-16
專業(yè): 無機合成和制備化學

[求助] gaussian09計算1-丁基-3-甲基咪唑陽離子電荷布居的問題

大家好，本人是量化新手，最近學習用gaussian09進行計算，在計算1-丁基-3-甲基咪唑陽離子的電荷布居時遇到了一些問題，想向大家請教一下。
我用了兩種方法，第一種是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，輸入文件沒有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
我的問題是
1）輸入文件中的sp test是什么意思？我在gaussian的文檔中沒有找到test關(guān)鍵詞。（輸入文件是我用chem3d保存構(gòu)型為gjf后，用gaussview載入的。setup選的是opt，tight，b3lyp和6-311++g(d,p)。gaussview自動在下邊添加了Additional Keywords：sp test。我試過如果沒有這兩個關(guān)鍵詞，程序無法運行）；
2）兩種方法算出來的電荷布居為何與直觀感覺相差很大？應該是咪唑環(huán)上正點強，氫上接近中性才比較符合直覺啊。而計算出來咪唑環(huán)是帶負電的。相比之下，我用Material Studio計算的Mulliken布居比較“正�！保ㄍ瑸閎3lyp）。是不是我基組和方法選擇的問題呢？
謝謝大家！

回復此樓

» 猜你喜歡

汕頭大學化學招調(diào)劑，相關(guān)內(nèi)容如下。已經(jīng)有4人回復
求助火焰封管的時候管子炸了已經(jīng)有1人回復
物理化學論文潤色/翻譯怎么收費? 已經(jīng)有84人回復
請問四氫呋喃溶解的聚合物用甲醇沉淀時，如何使沉淀過程加速？已經(jīng)有2人回復
七嗪類物質(zhì)合成求助已經(jīng)有0人回復
半透性容器進行高濕或者低濕研究已經(jīng)有0人回復
Chemical Bonding at Surfaces and Interfaces,最經(jīng)典的一本表面上化學相互作用教材已經(jīng)有0人回復
誠邀有志之士加入江蘇大學環(huán)境學院合成生物學課題組！(長期接收學碩專碩、碩博連讀) 已經(jīng)有16人回復
一志愿南大理論與計算化學方向（070300）342分，求調(diào)劑已經(jīng)有1人回復
歡迎化學、化工、物理、材料、或計算機背景的考研咨詢已經(jīng)有0人回復

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

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

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

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

凡士林護手霜

銅蟲 (正式寫手)

應助: 1 (幼兒園)
金幣: 1077.7
散金: 40
紅花: 1
帖子: 421
在線: 250.2小時
蟲號: 2911292
注冊: 2014-01-03
性別: GG
專業(yè): 污染控制化學

你好，請問你是做量子化學計算的研究嗎，你們單位購買了嗎

贊一下

回復此樓

夢想在心中，路在腳下。

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

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

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 foclnbrsic 的主題更新