计算化学公社

 找回密码 Forget password
 注册 Register

阳离子表面静电势分析判断亲核部位求助

查看数: 11689 | 评论数: 20 | 收藏 Add to favorites 3
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2019-8-29 11:58

正文摘要:

如图是通过Gaussview给出的带一个正电荷的离子的表面静电势图,该静电势图基本都是蓝色的,是不是说明该离子没有亲核部位,但是这个离子上面有一个羟基,羟基是可以参加反应的(应该是亲核部位)。所以这个静电势图 ...

回复 Reply

sobereva 发表于 Post on 2023-8-12 03:40:18
guoxiaomeng 发表于 2023-8-11 17:10
我想问一下您这个刻度是怎么修改的呢,我现在遇到了同样的问题了

哪里改色彩刻度范围下文明确说了
使用Multiwfn+VMD快速地绘制静电势着色的分子范德华表面图和分子间穿透图(含视频演示)
http://sobereva.com/443http://bbs.keinsci.com/thread-11080-1-1.html
guoxiaomeng 发表于 Post on 2023-8-11 17:10:17
我想问一下您这个刻度是怎么修改的呢,我现在遇到了同样的问题了
sobereva 发表于 Post on 2019-9-5 01:32:17
air18238016558 发表于 2019-9-4 09:38
卢老师,借助于这个求助话题,再请教您一个关于AIM的小问题,上传的表是一篇参考文献中剪下来的,不知道 ...

他的图单位都是a.u.,Multiwfn的显示临界点信息的界面里也明确说了



显然是一样的
air18238016558 发表于 Post on 2019-9-3 11:58:59
明白了,再次感谢卢老师!
sobereva 发表于 Post on 2019-9-3 00:44:48
air18238016558 发表于 2019-9-2 09:11
卢老师好,还得请教您,之前我可能没把问题说清楚,因为我是通过iso命令显示静电势着色等值面的,然后还 ...

那不叫极大/极小,而是最大/最小
色彩刻度上下限并不需要非得和实际的静电势上下限一致。刻意把范围设得比实际小一点,颜色差异会更鲜明。而且,强行和实际静电势最大、最小值弄成一致,色彩刻度的标签不够整齐、上下限数值不对称,在文中表达的时候也不方便。
air18238016558 发表于 Post on 2019-9-2 09:11:23
sobereva 发表于 2019-9-2 04:00
用ESPext和ESPpt的情况是不一样的
对于ESPext显示成透明等值面,因为VMD此时是基于静电势的cube文件里的 ...

卢老师好,还得请教您,之前我可能没把问题说清楚,因为我是通过iso命令显示静电势着色等值面的,然后还用了ext命令显示极大值和极小值,默认的颜色刻度单位为-0.03-0.03a.u.(0.03a.u.换算成kcal/mol后不到19kcal/mol),但是我的极大值有48kcal/mol,极小值有-43kcal/mol(这是一个中性分子,不是之前上传的图片上的阳离子),这样感觉颜色刻度值和极值不一致,如果把颜色刻度在VMD里由-0.03-0.03a.u.改成-0.08-0.08a.u.(0.08a.u.差不多对应50kcal/mol),iso命令显示的静电势着色等值面图片颜色又变得非常非常浅,感觉颜色层次不分明?谢谢老师的耐心指点!
sobereva 发表于 Post on 2019-9-2 04:00:46
air18238016558 发表于 2019-9-1 13:36
谢谢卢老师,已按照更改,但是发现一个问题,您说“ESPpt.vmd脚本里默认色彩刻度是-50kcal/mol ~ 50 kcal ...

用ESPext和ESPpt的情况是不一样的
对于ESPext显示成透明等值面,因为VMD此时是基于静电势的cube文件里的数据来着色的,而此cube文件里的单位是原子单位,所以色彩刻度的单位总是a.u.
而对于ESPpt通过表面顶点形式显示表面的情况,VMD是基于pdb文件的B因子那一列的静电势数据来着色的,中性(一般情况下)是以kcal/mol作为单位。而对于带净电荷的体系,由于静电势数量级较大,B因子那几列记录不下,所以(一般情况下)是以数值大小更小的eV作为单位。如果不确定的话,直接用文本编辑器打开pdb文件,看开头提示的单位便知。

颜色都为灰色说明色彩刻度设得不合适。根据Multiwfn做定量分子表面分析时候输出的表面静电势最小值和最大值,取当前作图对应的色彩刻度单位,在VMD里作为色彩刻度下限和上限就完了。这点在443文章末尾那个截图里已经体现得超级清楚了。
air18238016558 发表于 Post on 2019-9-1 13:36:25
本帖最后由 air18238016558 于 2019-9-1 16:20 编辑
sobereva 发表于 2019-8-31 22:16
阳离子体系,表面各个位置静电势都为正是非常正常的。比较相对大小就完了
另外,静电势预测的结果未必一 ...

谢谢卢老师,已按照更改,但是发现一个问题,您说“ESPpt.vmd脚本里默认色彩刻度是-50kcal/mol ~ 50 kcal/mol,对应颜色按照“蓝-白-红”过渡。如果想在VMD的界面里直接修改显示设定,诸如点的尺寸、色彩刻度范围等”,而我在VMD里界面看到的色彩刻度范围默认的是-0.03~0.03,这个单位应该是a.u.;如何把这个a.u.改为像您在196一文中所示范的(VMD里界面看到的色彩刻度范围默认的是kcal/mol)?在ESPiso.vmd脚本里更改为kcal/mol,这样通过ESPiso.bat计算后,图形颜色全为灰色了,不知道该如何改?还请卢老师指点,谢谢!
sobereva 发表于 Post on 2019-8-31 22:16:40
air18238016558 发表于 2019-8-31 20:53
图片重新上传了,还望卢老师指点!

阳离子体系,表面各个位置静电势都为正是非常正常的。比较相对大小就完了
另外,静电势预测的结果未必一定靠谱,可以用Multiwfn考察Hirshfeld电荷、考察福井函数f+等

PS:你的图作得太难看了。直接用文中的ext命令绘制成透明等值面图效果比这个好太多
sobereva 发表于 Post on 2019-8-31 20:34:19
air18238016558 发表于 2019-8-31 16:21
卢老师好,我已经按照您的方法通过multiwfn与VMD做出了静电势图(如图),和中性分子不一样,阳离子的静 ...

你传图方式不对,看不到。怎么上传图片看论坛置顶的新社员必读
air18238016558 发表于 Post on 2019-8-31 16:21:18
sobereva 发表于 2019-8-29 20:53
Grapihcs - color里对Name修改颜色(假设你是按照默认的,通过Name作为coloring的依据)

卢老师好,我已经按照您的方法通过multiwfn与VMD做出了静电势图(如图),和中性分子不一样,阳离子的静电势虽然有极大值和极小值,但都是正值,不知道这个该如何分析它们的亲核部位或亲电部位,还望老师指点,谢谢!E:H-R.tif
air18238016558 发表于 Post on 2019-8-29 22:06:45
sobereva 发表于 2019-8-29 20:53
Grapihcs - color里对Name修改颜色(假设你是按照默认的,通过Name作为coloring的依据)

好的,谢谢卢老师!

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

GMT+8, 2024-11-27 16:53 , Processed in 0.173482 second(s), 26 queries , Gzip On.

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