计算化学公社

 找回密码 Forget password
 注册 Register

求助福井函数f+值和其等值面图相矛盾

查看数: 243 | 评论数: 14 | 收藏 Add to favorites 2
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2026-1-26 18:33

正文摘要:

本帖最后由 moiraaa 于 2026-1-28 19:41 编辑 请问各位老师,之前拜读社长的文章,我计算的福井函数f+值和其等值面图相矛盾: 等值面图显示8,9,17,18,30,7,38,16,26,25,28,34,32,37是红色的区域, ...

回复 Reply

moiraaa 发表于 Post on 2026-2-2 09:49:13
sobereva 发表于 2026-2-2 03:40
某种福井函数越正的地方,暗示(但不一定真的准确)相应类型反应在此处发生的可能性越大。
也要注意局限 ...

好的,谢谢社长!!
sobereva 发表于 Post on 2026-2-2 03:40:10
moiraaa 发表于 2026-2-1 16:41
谢谢社长的解答
1蓝色是正的,红色是负的
2绿色是正,蓝色是负

某种福井函数越正的地方,暗示(但不一定真的准确)相应类型反应在此处发生的可能性越大。
也要注意局限性(不限于福井函数),参考量子化学波函数分析与Multiwfn程序培训班(http://www.keinsci.com/WFN)这页幻灯片我的总结:


moiraaa 发表于 Post on 2026-2-1 16:41:14
sobereva 发表于 2026-2-1 00:39
1、2 每次提问都要描述完整、严谨,先说清楚绿色和蓝色是什么含义,哪个是正那个是负

3 你载入的是哪 ...

谢谢社长的解答
1蓝色是正的,红色是负的
2绿色是正,蓝色是负
sobereva 发表于 Post on 2026-2-1 00:39:01
moiraaa 发表于 2026-1-31 17:19
社长,想请教您几个问题:
1.是不是f+,f-,f0的填色等值面都要看蓝色越深的地方,反应活性越大?
2.在vmd ...

1、2 每次提问都要描述完整、严谨,先说清楚绿色和蓝色是什么含义,哪个是正那个是负

3 你载入的是哪个态的fch文件,Multiwfn产生的电子密度cub文件就对应哪个态的电子密度

4、5 都用N态的电子密度格点数据定义电子密度的等值面
moiraaa 发表于 Post on 2026-1-31 17:19:36
本帖最后由 moiraaa 于 2026-1-31 17:43 编辑
sobereva 发表于 2026-1-31 05:17
N.fch和N+1.fch里的原子序号根本不一样,进入Multiwfn主功能0一看便知
先确保所有波函数文件里相同序号 ...

社长,想请教您几个问题:
1.是不是f+,f-,f0的填色等值面都要看蓝色越深的地方,反应活性越大?
2.在vmd中看f+,f-,f0是不是就是看绿色的部分,反应活性越大?
3.画f+的等值面图时,按照博文402的操作得到的第一个density.cub是不是N+1.cub?
4.画f+的填色等值面图时,是把N+1.cub拖到gview里,载入f+吗?
5.画f0的填色等值面,是不是N.cub拖到gview里,载入f0吗
moiraaa 发表于 Post on 2026-1-31 10:11:33
sobereva 发表于 2026-1-31 05:17
N.fch和N+1.fch里的原子序号根本不一样,进入Multiwfn主功能0一看便知
先确保所有波函数文件里相同序号 ...

好的,社长,真的是一语道破梦中人,感谢社长
sobereva 发表于 Post on 2026-1-31 05:17:08
moiraaa 发表于 2026-1-30 21:01
社长,我去做了f+的等值面图,基本上吻合简缩福井函数但是22N那里的蓝色很少,这个点不是很吻合,请问社 ...

N.fch和N+1.fch里的原子序号根本不一样,进入Multiwfn主功能0一看便知
先确保所有波函数文件里相同序号的原子的坐标都一样再说,否则结果都是错的
moiraaa 发表于 Post on 2026-1-29 16:57:53
sobereva 发表于 2026-1-29 16:17
你先看Multiwfn直接显示的f+的等值面图,像下文那样,这最清楚直观。把f+等值面图的特征和简缩福井函数对上 ...

好的,谢谢社长,我先去做等值面图
sobereva 发表于 Post on 2026-1-29 16:17:29
你先看Multiwfn直接显示的f+的等值面图,像下文那样,这最清楚直观。把f+等值面图的特征和简缩福井函数对上了,再去试图对照f+等值面图和分子表面f+填色图。
使用Multiwfn超级方便地计算出概念密度泛函理论中定义的各种量
http://sobereva.com/484http://bbs.keinsci.com/thread-13328-1-1.html

当前你的格点数据看起来可能有问题。先确保能完整重复下文的福井函数填色图例子,避免细节弄错了
基于Multiwfn产生的cube文件在VMD和GaussView中绘制填色等值面图的方法
http://sobereva.com/402http://bbs.keinsci.com/thread-7756-1-1.html
我现在不方便检查你的文件

sobereva 发表于 Post on 2026-1-29 03:48:56
moiraaa 发表于 2026-1-28 19:44
谢谢社长的纠正,已经修改帖子。
色彩刻度是0.001和-0.001
因为同事要拼这张图,所以就算了那个简缩福 ...

色彩刻度要体现颜色变化方式
图片放到文章里必须给出color bar

“拼这张图”语义不明
都已经说了不建议用f+
简缩福井函数总能和分子表面福井函数填色图(恰当选取电子密度等值面数值和色彩刻度范围的情况下)对应上。另外,福井函数等值面图比填色图更常用
sobereva 发表于 Post on 2026-1-28 15:16:22
表格里那叫简缩福井函数不叫福井函数,认真看下文
计算化学中的一些常见不良写法和用词
http://sobereva.com/298http://bbs.keinsci.com/thread-1358-1-1.html
亲电取代反应中活性位点预测方法的比较(http://www.whxb.pku.edu.cn/CN/abstract/abstract28694.shtml
Comparative study on the methods for predicting the reactive site of nucleophilic reaction(https://link.springer.com/article/10.1007%2Fs11426-015-5494-7

gview里显示的是分子表面的f+的填色等值面图,不能叫等值面图,这和f+的等值面图完全是两码事。

若不给出色彩刻度,填色图则毫无意义。

另外,f+预测亲核反应位点本来就比较烂,准确度差,仔细看
Comparative study on the methods for predicting the reactive site of nucleophilic reaction(https://link.springer.com/article/10.1007%2Fs11426-015-5494-7
TCA上的一篇对比不同原子电荷预测反应位点、亲电/亲核性的文章(http://bbs.keinsci.com/thread-15512-1-1.html

还不如用Multiwfn算Hirshfeld电荷、LEAE(下文)、双描述符等讨论
使用Multiwfn通过局部电子附着能(LEAE)考察亲核反应位点、难易及弱相互作用
http://sobereva.com/676http://bbs.keinsci.com/thread-38416-1-1.html

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-17 01:11 , Processed in 0.204648 second(s), 25 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list