请选择 进入手机版 | 继续访问电脑版

计算化学公社

 找回密码
 现在注册!
查看: 115|回复: 4

[综合交流] ACE-Tyrosine-NME + glucose的RESP电荷计算问题

[复制链接]

84

帖子

0

威望

1131

eV
积分
1215

Level 4 (黑子)

发表于 2020-11-18 22:27:38 | 显示全部楼层 |阅读模式
计算的分子是ACE-Tyrosine-NME + glucose,原子index分别是1-6 (ACE), 7-26 (Tyrosine), 27-32 (NME), 33-54 (glucose)
一共有五个构象,由ORCA优化得来,Multiwfn版本是3.7 linux版,做了电荷constraints,试过2次,对结果有一些疑问。
constraints1 : Tyrosine部分电荷为-0.194, glucose部分 0.194 (见附件resp_contraints1.out.gz),得到的结果是
Tyrosine的电荷和是-0.196986,glucose部分电荷是0.19356

glucose部分还好,Tyrosine部分的电荷的差异是不是有点大,还是正常的,还是可以设定更强的constraints?

constraints 2:在constraints 1 的基础上,添加了ACE,NME都为0,这次的计算都没有收敛,得到的结果肯定是不可靠的。
是不是不能这样添加constraints呢?(resp_constraints2.out.gz)

谢谢!

resp_constraints2.out.gz (7.51 KB, 下载次数: 3)

2万

帖子

25

威望

3万

eV
积分
60837

管理员

公社社长+计算化学玩家

发表于 2020-11-19 06:09:17 | 显示全部楼层
我不知道trysine和glucose是怎么相连的,如果是非共价方式形成的复合物,不应当把二者弄到一起算RESP电荷。

我没有你的具体输入文件有些情况不好判断。

对于情况1,也尝试不考虑多构象,或者用one-stage方式拟合看看是什么情况

把ACE和NME分别约束成0原理上没问题,但鉴于当前难收敛因而结果不能用,尝试把二者总和设为0,或者尝试其一设为0。

PS:当前上传的out文件本身很小,没必要压缩
北京科音自然科学研究中心http://www.keinsci.com  致力于计算化学的发展和传播,长期开办最高水准的各种量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训,是提升计算化学研究水平的最佳选择。欢迎加入“北京科音”公众号获取培训最新消息和计算化学资讯!培训相关信息见《北京科音办的培训班FAQ》(http://bbs.keinsci.com/thread-5098-1-1.html)。
欢迎加入人气最高、水准最高的综合性理论与计算化学交流QQ群:思想家公社QQ群1号:18616395,2号:466017436。合计6000人。两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一定会被拒绝加入。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

84

帖子

0

威望

1131

eV
积分
1215

Level 4 (黑子)

 楼主| 发表于 2020-11-19 08:13:13 | 显示全部楼层
本帖最后由 beyond 于 2020-11-19 08:15 编辑
sobereva 发表于 2020-11-19 06:09
我不知道trysine和glucose是怎么相连的,如果是非共价方式形成的复合物,不应当把二者弄到一起算RESP电荷。 ...

谢谢SOB! tyrosine的OH是与glucose的C1共价相连的!

我尝试了一下只用一个molden文件,另外ACE与NME的电荷限制为0,Tyrsine为-0.194,glucose为0.194,
得到的结果还是不合理,出现了多个原子电荷为0的情况,如ACE的C=O基团。

其实在对另外一个分子UDP-Glucose做电荷限制的RESP拟合,也出现类似情况.

附件是一个ACE-Tyrosine-NME-Glucose的molden文件, Tyr-Glc_5ps.molden.gz (9 MB, 下载次数: 4)

2万

帖子

25

威望

3万

eV
积分
60837

管理员

公社社长+计算化学玩家

发表于 2020-11-20 05:32:19 | 显示全部楼层
beyond 发表于 2020-11-19 08:13
谢谢SOB! tyrosine的OH是与glucose的C1共价相连的!

我尝试了一下只用一个molden文件,另外ACE与NME ...

注意进入选项6设置电荷约束的时候屏幕上有提示
Note: For standard two-stage RESP fitting, below options only take effect for the first stage
手册里对选项6的介绍也写了
Option 6: This option is use to set charge constraint in "one-stage fitting" or the first step of "two-stage RESP fitting".
因此对于两步式拟合,第二步的时候自动不用电荷约束,这是因为自动判断的等价性约束和自定义的电荷约束很容易发生冲突。我估计这个应该是造成两步式拟合结果稍微偏离0.194的原因。


这是我用一步拟合得到的输出,没有约束1-6,27-32 0,结果中7~26的加和精确为0.194
out.txt (9.63 KB, 下载次数: 1)
北京科音自然科学研究中心http://www.keinsci.com  致力于计算化学的发展和传播,长期开办最高水准的各种量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训,是提升计算化学研究水平的最佳选择。欢迎加入“北京科音”公众号获取培训最新消息和计算化学资讯!培训相关信息见《北京科音办的培训班FAQ》(http://bbs.keinsci.com/thread-5098-1-1.html)。
欢迎加入人气最高、水准最高的综合性理论与计算化学交流QQ群:思想家公社QQ群1号:18616395,2号:466017436。合计6000人。两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一定会被拒绝加入。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

84

帖子

0

威望

1131

eV
积分
1215

Level 4 (黑子)

 楼主| 发表于 7 天前 | 显示全部楼层
sobereva 发表于 2020-11-20 05:32
注意进入选项6设置电荷约束的时候屏幕上有提示
Note: For standard two-stage RESP fitting, below opt ...

谢谢SOB! 改用一步拟合之后,我也得到了预期的结果!

那我就改用你帖子里面的方法来做吧,有问题再请教!
您需要登录后才可以回帖 登录 | 现在注册!

本版积分规则

手机版|北京科音自然科学研究中心|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949-1号 )

GMT+8, 2020-11-28 16:35 , Processed in 0.154276 second(s), 26 queries .

快速回复 返回顶部 返回列表