| 查看: 865 | 回复: 2 | |||
[求助]
能否用multiwfn计算重叠布居电荷
|
|
我又厚脸皮来了，还是关于zhou2009老师的“关于HF、H2O的电荷计算”帖子。 原帖中先是用multiwfn做出了Δρ图，据图分析了分子轨道的性质，然后用GsGrid计算各个轨道的重叠布居电荷。最近multiwfn加了盆积分功能，我想试试用multiwfn代替GsGrid计算上述重叠布居电荷。(不想重新学个软件了) 我的做法是，首先按之前sobereva专家给出的方法做Δρ图，只是在“减掉波函数”时选择功能5而不是4，生成格点文件。再选择功能17，对内存里的格点文件计算盆积分。 我用hf/sto-3g算了水分子的mo3和mo4： mo3的盆积分得到了两个积分值为0.32的盆 mo4的盆积分得到了一个积分值为0.29，一个为-0.32的盆 zhou2009老师帖子里的计算结果为 “MO3中H上电荷为0.2649。从它的Δρ得到所谓重叠布居电荷为0.6298” “MO4的O、H1、H2三者之间的重叠布居电荷为0.3293” 这样来看，mo3我得到的两个盆积分加和大约就是zhou2009老师帖子里的计算结果。而mo4就不好说了，我这里两个盆积分值加和，正负抵消成-0.03了。不知道-0.32的盆是不是就对应了O的电荷集中转移到2H上的值，可不可以理解成zhou2009老师帖子里说的mo4上的OH重叠布局全划分给H了，那么-0.32和zhou老师给出的0.3293对应呢？ 我自己无法准确判断，主要是因为我用高斯计算的水分子的mulliken布居分析结果和zhou2009老师帖子里的不太一致，那么就没法和帖子里的数据直接对应了，特来发问，求指教。 另附我计算的水分子优化后的坐标和mulliken布居分析，高斯版本是09W A02 Standard orientation: --------------------------------------------------------------------- Center Atomic Atomic Coordinates (Angstroms) Number Number Type X Y Z --------------------------------------------------------------------- 1 8 0 0.000000 0.000000 0.127185 2 1 0 0.000000 0.758222 -0.508739 3 1 0 0.000000 -0.758222 -0.508739 --------------------------------------------------------------------- Full Mulliken population analysis: 1 2 3 4 5 1 1 O 1S 2.10788 2 2S -0.10780 2.00685 3 2PX 0.00000 0.00000 2.00000 4 2PY 0.00000 0.00000 0.00000 0.75079 5 2PZ 0.00000 0.00000 0.00000 0.00000 1.17276 6 2 H 1S -0.00112 -0.02504 0.00000 0.16111 0.11889 7 3 H 1S -0.00112 -0.02504 0.00000 0.16111 0.11889 6 7 6 2 H 1S 0.62634 7 3 H 1S -0.04537 0.62634 |
回复此楼» 收录本帖的淘帖专辑推荐
 SOB 集锦 |
» 猜你喜欢
湖南农业大学能源材料创新团队招收2025化学专业研究生（学硕）调剂
已经有0人回复
-
2026年西华大学材料学院-材料近净成形与表面技术研究团队-招收调剂研究生
已经有10人回复
-
物理化学论文润色/翻译怎么收费?
已经有243人回复
-
调剂信息
已经有0人回复
-
河北农大生物与医药专业招收调剂，名额充足
已经有0人回复
-
齐鲁工业大学 国家重点实验室-新能源、锂离子电池、材料、化学方向
已经有0人回复
-
杰青/长江团队招收调剂
已经有0人回复
-
中南林业科技大学罗勇锋教授团队2026年招收材料类调剂研究生
已经有1人回复
-
电子科技大学材料学院王翀教授拟招收工程博士一名（化学、材料方向）
已经有2人回复
-
湖南农业大学能源材料创新团队招收2026化学专业研究生（学硕）调剂
已经有6人回复
-
奥尔堡Xinxin Xiao课题组招收一名岗位制博士（酶基生物电化学）
已经有0人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 分子中电荷变化种种
已经有82人回复
- 化学键重叠布居数分析中的负值问题
已经有3人回复
- 建议multiwfn开发GUI界面 文件打开什么的效率太低了
已经有10人回复
- Gaussian计算 Mulliken overlap populations 和 bond overlap population
已经有8人回复
- 如何用Multiwfn软件做出P DOS?
已经有4人回复
- 使用Multiwfn绘制原子轨道图形、研究原子壳层结构及相对论效应的影响
已经有17人回复
- 使用Multiwfn绘制NBO及相关轨道
已经有30人回复
- 【求助】NBO轨道指的是nbo电荷的分布情况吗？
已经有3人回复
- 用Multiwfn、Gaussview、Molekel、VMD观看龙虾、盆景、骨盆、大脑
已经有16人回复
- 使用Multiwfn做拓扑分析以及计算孤对电子角度
已经有18人回复
- 回答一些关于Multiwfn的疑问以及未来Multiwfn的发展打算
已经有37人回复
- GAUSSIAN怎么计算overlap matrix？
已经有31人回复
- 谈谈Molekel做静电势填色等值面图的方法及误区
已经有50人回复
- 关于分子间范德华作用的随想
已经有68人回复
- 【求助】castep计算电荷布局数分析
已经有9人回复
- 【zhou2009个人文集】关于电荷计算的第二组分子
已经有4人回复
- 【zhou2009个人文集】关于HF、H2O的电荷计算
已经有15人回复
- 【sobereva个人文集】电子定域性的图形分析
已经有56人回复
sobereva
至尊木虫 (著名写手)
本人已永久离开小木虫
- QC强帖: 17
- 应助: 226 (大学生)
- 贵宾: 5.4
- 金币: 2180.5
- 红花: 217
- 帖子: 1241
- 在线: 366.8小时
- 虫号: 281271
- 注册: 2006-09-23
- 专业: 宇宙学
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
gmy1990: 金币+5 2013-04-01 20:13:10
klaus1987: 金币+10, ★★★★★最佳答案 2013-04-02 09:27:15
感谢参与，应助指数 +1
gmy1990: 金币+5 2013-04-01 20:13:10
klaus1987: 金币+10, ★★★★★最佳答案 2013-04-02 09:27:15
|
你不需要用主功能17来做盆积分。只要你用主功能5计算完了格点数据，屏幕上直接就会输出当前格点数据在整个空间的积分值。 如果你按照我之前给你回复的方法，但只是将主功能4改成了主功能5，那么这个积分值就是MO2当中O和两个氢这两个片段间的重叠布居值。 高斯输出的那个矩阵是考虑了所有轨道的，而且内容和你要分析的也不是一回事。 gsgrid早已并入Multiwfn作为主功能13了。主功能13比最后一个版本gsgrid的功能也要强大很多。所以没必要学也没必要用gsgrid。 如果你只想要得到某个MO当中某两个原子间的重叠布居值，那么最简单的方法就是载入fch文件，然后选7，再选5，之后选2。这时会输出gross atomic population matrix，这和高斯输出的Condensed to atoms矩阵是一致的，但这是考虑了所有轨道的结果。接下来选y，就会将每个轨道的贡献输出到当前目录的groatmdcp.txt当中。从中可以找到例如 ************ Orbital 2 Occ: 2.00000000 ************ 1 2 3 1 1.36473981 0.12815040 0.12815040 2 0.12815040 0.05514825 0.00618105 3 0.12815040 0.00618105 0.05514825 4 1.62104062 0.18947970 0.18947970 （最后一行是相应列的加和） 如果你想要得到MO2当中1和3号原子的重叠布居，那么就是0.12815040乘以2。 Multiwfn也可以获得特定MO当中基函数之间的重叠布居。对于Mulliken布居的分析，Multiwfn是现有任何程序中最强大的。 |
2楼2013-04-01 18:07:17
3楼2013-04-02 09:27:01
