| 查看: 2103 | 回复: 7 | |||
nnipp金虫 (小有名气)
|
[求助]
如何得到指定原子上的电荷密度差
|
| 如题，得到这个体系的电子密度cube以后，如何指定原子得到某些原子上的电子密度差分图？ |
回复此楼» 收录本帖的淘帖专辑推荐
 量化计算 |
» 猜你喜欢
伦敦大学学院Benjamin Moss博士招收人工智能、原位光谱、催化相关26年博士生
已经有11人回复
-
钴酸锂半电池小倍率容量上不去
已经有1人回复
-
物理化学论文润色/翻译怎么收费?
已经有192人回复
-
吉林大学材料物理本科生求问调剂信息
已经有23人回复
-
2026第二届光电子与半导体器件前沿技术研讨会——光电子赋能·半导体创芯！
已经有1人回复
-
光电子赋能·半导体创芯！
已经有1人回复
-
光电子赋能·半导体创芯！
已经有0人回复
-
QE计算电声耦合的时候报错Error in routine lambda (100)wrong or too many modes
已经有2人回复
-
中国科学院成都山地灾害与环境研究所2026年特别研究助理（博士后）招聘启事
已经有0人回复
-
九江学院2026年最新高层次人才招聘公告
已经有0人回复
-
深圳信息职业技术学院2026年教师招聘公告（最新）
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 电荷差分密度
已经有9人回复
- 再次求问，差分电荷密度图的定性分析
已经有17人回复
- 两个原子的局域态密度在费米处重叠
已经有4人回复
- 怎么使用bader电荷分析出原子轨道上的电子数，或者其它能够得出轨道电子的软件？
已经有29人回复
- Gaussian中原子净电荷的指定
已经有3人回复
- 【求教】电子密度，差分密度，原子电荷这三个之间有什么关系？
已经有7人回复
- 请教有机分子中O原子和H原子电荷数怎么计算，O-H键键级怎么计算！谢谢！
已经有4人回复
- 请问：如何l利用xcrysden画电荷密度
已经有5人回复
- 【求助】差分电荷密度
已经有19人回复
- 【求助】原子的电荷密度分布怎么求？
已经有3人回复
- 【zhou2009个人文集】关于HF、H2O的电荷计算
已经有15人回复
zhou2009
版主 (著名写手)
- QC强帖: 15
- 应助: 337 (大学生)
- 贵宾: 4.053
- 金币: 42654.6
- 红花: 153
- 帖子: 1729
- 在线: 1274.7小时
- 虫号: 676800
- 注册: 2008-12-18
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 量子化学
2楼2013-03-29 17:11:06
nnipp
金虫 (小有名气)
- 应助: 3 (幼儿园)
- 金币: 1003.5
- 散金: 16
- 红花: 3
- 帖子: 172
- 在线: 234.5小时
- 虫号: 376554
- 注册: 2007-05-19
- 性别: GG
- 专业: 理论和计算化学
3楼2013-03-29 17:50:44
zhou2009
版主 (著名写手)
- QC强帖: 15
- 应助: 337 (大学生)
- 贵宾: 4.053
- 金币: 42654.6
- 红花: 153
- 帖子: 1729
- 在线: 1274.7小时
- 虫号: 676800
- 注册: 2008-12-18
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 量子化学
4楼2013-03-29 18:19:28
sobereva
至尊木虫 (著名写手)
本人已永久离开小木虫
- QC强帖: 17
- 应助: 226 (大学生)
- 贵宾: 5.4
- 金币: 2180.5
- 红花: 217
- 帖子: 1241
- 在线: 366.8小时
- 虫号: 281271
- 注册: 2006-09-23
- 专业: 宇宙学
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
zhou2009: 金币+5 2013-03-30 12:10:37
nnipp: 金币+20, ★★★★★最佳答案, 谢谢耐心回复 2013-03-30 23:08:40
感谢参与，应助指数 +1
zhou2009: 金币+5 2013-03-30 12:10:37
nnipp: 金币+20, ★★★★★最佳答案, 谢谢耐心回复 2013-03-30 23:08:40
|
在Multiwfn中，载入波函数文件后，可以利用选项-3来只把感兴趣的波函数信息保留。 假设你要对A和B两个状态求密度差，可以先读入B的波函数文件，然后进入-3，输入活性中心部分的原子数和其中的各个原子编号。然后进入主功能6选0将波函数保存为new.wfn，其中就没有了周边原子的信息。 然后再读入A，也是进入选项-3并且输入活性中心部分的原子数和其中的各个原子编号。然后进入主功能5来照常计算密度差格点数据，其中被减的就是new.wfn。这样，所得格点文件就只有活性中心那部分的密度差了。 注意A和B的几何结构必须一致。并且，最好不要给活性中心的原子用弥散函数。否则活性中心的密度差可能会延伸出去 另外，也可以按照普通方法计算密度差格点数据，然后把cube文件载入Multiwfn，之后选主功能13里的选项13。这个功能可以自行设定距离指定的一批原子较远的空间内的格点数据的数值，如果数值设为0，而且原子范围选的是活性中心区域的，就相当于把活性中心区域以外的格点数据内容都删了，作图时就看不出来了。 无论如何，用上述两种时办法务必搞清其物理意义。 |
5楼2013-03-30 12:07:55
nnipp
金虫 (小有名气)
- 应助: 3 (幼儿园)
- 金币: 1003.5
- 散金: 16
- 红花: 3
- 帖子: 172
- 在线: 234.5小时
- 虫号: 376554
- 注册: 2007-05-19
- 性别: GG
- 专业: 理论和计算化学
6楼2013-03-30 23:06:54
zhou2009
版主 (著名写手)
- QC强帖: 15
- 应助: 337 (大学生)
- 贵宾: 4.053
- 金币: 42654.6
- 红花: 153
- 帖子: 1729
- 在线: 1274.7小时
- 虫号: 676800
- 注册: 2008-12-18
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 量子化学
|
活性物种A与骨架B有氢键作用时，如果要分析B对A的影响，可以做AB与A、B的密度差。这又分两种情况： 1、求AB的wfn文件AB. Wfn，象上面sobereva博士讲的那样，在multiwfn中将AB分成两个分子片，得到A. wfn和B. wfn，做AB. Wfn-A. wfn-B. wfn。因为AB与A、B都是取的AB的MO系数，原子上的密度都抵消了，这样得到的只是两个分子片之间的重叠布居电荷的情况，没有掩盖的问题。 2、求AB的wfn文件AB. Wfn，然后从AB的分子坐标中，分别摘取A、B的坐标。得到A、B的坐标后，对A、B作为对立的分子进行计算，分别得到A. wfn和B. wfn，再作AB. Wfn-A. wfn-B. wfn，得到电子密度差。这时是将A、B作为对立的分子进行计算的，活性物种A与骨架B的刚性部分也是抵消了，得到的只是B对A的影响所引起的电子密度净变化。 3、如果是自己设计的AB坐标，A、B应是分开的，分别删除部分坐标形成A、B的坐标即可。如果是AB坐标从晶体中来，A、B的坐标是混杂的，需要在GV中将图展开，标记上原子的序列号，识别并记录A、B原子的坐标编号，在GV中进入坐标的表格表示，分两次从后往前按编号逐一删去A得B，删去B得A。 4、进一步可以将活性物种A与骨架B和被催化分子C作为一个系综，求BAC. Wfn-BA. wfn-C. wfn，这时可以看见催化剂BA对C的影响。 要注意的是：BA与C仅是范德华作用时，要考虑你的计算方法是否能胜任计算范德华作用。 5、分析2、4的电子密度差，可能很复杂，要将图像摆好方位，以便识别，还要结合MO系数来考察、认定。 [ Last edited by zhou2009 on 2013-3-31 at 12:05 ] |
» 本帖已获得的红花（最新10朵）
nnipp7楼2013-03-31 10:02:25
nnipp
金虫 (小有名气)
- 应助: 3 (幼儿园)
- 金币: 1003.5
- 散金: 16
- 红花: 3
- 帖子: 172
- 在线: 234.5小时
- 虫号: 376554
- 注册: 2007-05-19
- 性别: GG
- 专业: 理论和计算化学
8楼2013-03-31 13:12:23
10