24小时热门版块排行榜

南方科技大学公共卫生及应急管理学院2026级博士研究生招生报考通知（长期有效）

返回列表

foclnbrsic

铜虫 (初入文坛)

应助: 0 (幼儿园)
金币: 61.5
帖子: 4
在线: 11.9小时
虫号: 1758794
注册: 2012-04-16
专业: 无机合成和制备化学

[求助] gaussian09计算1-丁基-3-甲基咪唑阳离子电荷布居的问题

大家好，本人是量化新手，最近学习用gaussian09进行计算，在计算1-丁基-3-甲基咪唑阳离子的电荷布居时遇到了一些问题，想向大家请教一下。
我用了两种方法，第一种是NBO，输入文件和结果如下：
%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
——————————————————下面是结果
                                    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，其他相同，结果是：
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关键词。（输入文件是我用chem3d保存构型为gjf后，用gaussview载入的。setup选的是opt，tight，b3lyp和6-311++g(d,p)。gaussview自动在下边添加了Additional Keywords：sp test。我试过如果没有这两个关键词，程序无法运行）；
2）两种方法算出来的电荷布居为何与直观感觉相差很大？应该是咪唑环上正点强，氢上接近中性才比较符合直觉啊。而计算出来咪唑环是带负电的。相比之下，我用Material Studio计算的Mulliken布居比较“正常”（同为b3lyp）。是不是我基组和方法选择的问题呢？
谢谢大家！

回复此楼