| 查看: 1159 | 回复: 1 | |||
[求助]
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）。是不是我基组和方法选择的问题呢？ 谢谢大家！ |
回复此楼» 猜你喜欢
中科院西安光机所——聚焦功能材料创新｜诚邀参与 2026 年度学术研讨会。
已经有1人回复
-
中科院西安光机所——促进学术交流，推动产学研融合，分享领域成果，加强校企合作。
已经有1人回复
-
物理化学论文润色/翻译怎么收费?
已经有128人回复
-
中科院西安光机所——促进学术交流，推动产学研融合，分享领域成果，加强校企合作。
已经有0人回复
-
Comsol 6.0 桶内肥料产热产湿过程模拟求助
已经有1人回复
-
分子动力学 Plumed教程求助
已经有0人回复
-
研究生学术新锐奖｜宁波微纳材料研讨会参会通知
已经有2人回复
-
申请2027年材料/化学类博士
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 怎么计算 atomic spin density
已经有4人回复
- 请问Gaussian计算frequency出现下面的问题，如何解决？
已经有11人回复
- 电子轨道能量用什么方法计算比较好？
已经有6人回复
- 【小卒个人文集】写给初学者:小议形式电荷对量化初学者的误导
已经有27人回复
- 【求助】离子键和共价键的区分
已经有17人回复
凡士林护手霜
铜虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 1077.7
- 散金: 40
- 红花: 1
- 帖子: 421
- 在线: 250.2小时
- 虫号: 2911292
- 注册: 2014-01-03
- 性别: GG
- 专业: 污染控制化学
2楼2016-04-26 13:49:43
100