请选择 进入手机版 | 继续访问电脑版
第八届北京科音分子动力学与GROMACS培训班将于2021年1月18~21日于北京举办,请点击此链接查看详情,这是系统性学习分子动力学模拟、掌握GROMACS程序使用的最好机会!(再下一届预计于3或4月举办)

计算化学公社

 找回密码
 现在注册!
查看: 223|回复: 5

[Multiwfn使用咨询] 关于使用Mutiwfn计算RESP的若干问题

[复制链接]

15

帖子

1

威望

89

eV
积分
124

Level 2 能力者

发表于 2020-11-29 21:58:12 | 显示全部楼层 |阅读模式
本帖最后由 YZH_jesse 于 2020-11-29 22:04 编辑

我阅读完Sobereva老师《R原子电荷的思想以及在Multiwfn中的计算》(http://bbs.keinsci.com/thread-16190-1-1.html) 后尝试了一下甲醛分子RESP2的计算,以下为我的计算步骤及计算过程中出现的几个主要问题:

第一步:在gview中随便画一个甲醛分子的初始结构,保存为 .gjf 格式的Gaussian输入文件,修改其内容为:
  1. %nprocshared=6
  2. %mem=2GB
  3. %chk=C:\Users\yxdhx\Desktop\Charge Distribution\Formaldehyde\opt.chk
  4. # B3LYP/TZVP em=GD3BJ opt

  5. niconiconi

  6. 0 1
  7. C                  0.52701500   -0.00001300    0.00042100
  8. H                  1.11571800   -0.93807900   -0.00084900
  9. H                  1.11357700    0.93925600   -0.00084800
  10. O                 -0.67392300   -0.00013700   -0.00010400

  11. --link1--
  12. %oldchk=C:\Users\yxdhx\Desktop\Charge Distribution\Formaldehyde\opt.chk
  13. %chk=C:\Users\yxdhx\Desktop\Charge Distribution\Formaldehyde\SP_gas.chk
  14. # B3LYP/def2TZVP em=GD3BJ geom=allcheck


  15. --link1--
  16. %oldchk=C:\Users\yxdhx\Desktop\Charge Distribution\Formaldehyde\opt.chk
  17. %chk=C:\Users\yxdhx\Desktop\Charge Distribution\Formaldehyde\SP_solv.chk
  18. # B3LYP/def2TZVP em=GD3BJ scrf=solvent=ethanol geom=allcheck
复制代码
第二步:将Gaussian计算得到的两个 .chk 文件通过 formchk.exe 程序转为 .fch 文件,输入到Multiwfn中,依次选择 7 / 18 / 1 来计算RESP,得到的结果与Sobereva老师博文中所得到结果的对比如下:
(PS:其中 Test1 为我画的第一个初始结构的计算结果,Test2为我画出的第二个初始结构的计算结果,Test sob为我用当前Multiwfn版本并使用http://bbs.keinsci.com/thread-16190-1-1.html中给出的甲醛初始结构算得的结果,最后一个为sob老师在网站上直接给出的计算结果,意在对比初始结构对于计算结果的影响)


问题1:初始结构对于计算结果还是会产生影响的(小数点后第三位),而且由于Lammps中给出的TIP3P水模型中 O 和 H 原子的电荷均精确到小数点后第三位(-0.830/0.415),所以我觉得这个误差还是有必要去深究一下;而且很有趣的是我随意画出的两个初始结构的计算结果 和sob老师给出的初始结构的计算结果有明显差异,我不理解这个差异产生的原因是什么(三者除了初始原子坐标不同 其他条件保持严格相同,而且从 .gjf 第一部分可以看出是先对甲醛分子进行过结构优化的)

问题2:依次选择 7 / 18 / 1 计算RESP后,均提示了 Stage 2 of standard RESP fitting is skipped since no atom needs to be fitted (包括sob老师在博文中给出的 .gjf 文件内容),但是更神奇的是,在同样不加额外的约束的前提下,我自己构建的模型中H原子的电荷并不是完全一致的,但是sob老师的模型中H原子RESP的计算结果就是完全一致的(可由插图中看出),又是一个不能理解的地方。

希望各位老师们能给予我这个巨菜无比的初学者一点“神谕”











计算结果对比图

计算结果对比图

2万

帖子

25

威望

3万

eV
积分
62470

管理员

公社社长+计算化学玩家

发表于 2020-11-30 00:52:07 | 显示全部楼层
1 从笛卡尔坐标上看,你的H2CO明显不满足C2v点群。有点群对称性的体系在计算前一定要先在GaussView里对称化成实际的点群,这样算得又快又精确又容易收敛

2 也是因为你的体系不满足C2v对称性所致
北京科音自然科学研究中心http://www.keinsci.com  致力于计算化学的发展和传播,长期开办最高水准的各种量子化学、分子动力学与GROMACS、波函数分析与Multiwfn程序等主题的培训,是快速提升计算化学研究水平的最佳途径。欢迎加入“北京科音”公众号获取培训最新消息和计算化学资讯!培训相关信息见《北京科音办的培训班FAQ》
欢迎加入人气最高、水准最高的综合性理论与计算化学交流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!

15

帖子

1

威望

89

eV
积分
124

Level 2 能力者

 楼主| 发表于 2020-11-30 10:24:08 | 显示全部楼层
sobereva 发表于 2020-11-30 00:52
1 从笛卡尔坐标上看,你的H2CO明显不满足C2v点群。有点群对称性的体系在计算前一定要先在GaussView里对称化 ...

感谢sob老师  在GaussView中设置对称化后很轻松地解决了上述问题。
虽然计算结果的精度满足我的需求,但是和您给出的甲醛分子初始笛卡尔坐标算得的结果并不是严格相同(小数点后第四位有差异)(由相同版本的Multiwfn算得),这个差异可以忽略不计吗?


屏幕截图 2020-11-30 102147.png

2万

帖子

25

威望

3万

eV
积分
62470

管理员

公社社长+计算化学玩家

发表于 2020-11-30 12:10:09 | 显示全部楼层
YZH_jesse 发表于 2020-11-30 10:24
感谢sob老师  在GaussView中设置对称化后很轻松地解决了上述问题。
虽然计算结果的精度满足我的需求,但 ...

你还是没有充分对称化,两个氢的电荷还是有差异。
北京科音自然科学研究中心http://www.keinsci.com  致力于计算化学的发展和传播,长期开办最高水准的各种量子化学、分子动力学与GROMACS、波函数分析与Multiwfn程序等主题的培训,是快速提升计算化学研究水平的最佳途径。欢迎加入“北京科音”公众号获取培训最新消息和计算化学资讯!培训相关信息见《北京科音办的培训班FAQ》
欢迎加入人气最高、水准最高的综合性理论与计算化学交流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!

15

帖子

1

威望

89

eV
积分
124

Level 2 能力者

 楼主| 发表于 2020-12-1 23:39:41 | 显示全部楼层
sobereva 发表于 2020-11-30 12:10
你还是没有充分对称化,两个氢的电荷还是有差异。

sob老师  Test 3 是我对称化之后计算的结果  H原子的电荷是一致的  (Test 1和Test 2是我之前的计算结果  我在回复中没有写清楚)
现在我在尝试计算 间苯二胺(MPD)的RESP,还是存在很多问题。

我的计算流程为:

1. 在GaussView中构建一个MPD分子,在点群对称中设置C2v对称(如果存在这个等级的话),保存为 .gjf 格式的文件,关键词为:(目的是对这一初始构象进行几何优化)
  1. %mem=2GB
  2. %nprocshared=6
  3. %chk=C:\MPD.chk
  4. # opt b3lyp/6-311g(d,p) em=gd3bj
复制代码

2. 将优化之后的结果分别进行四种不同的处理:

    2.1 直接将优化之后的结果放到Multiwfn中进行RESP计算 (7 / 18 / 1)

    2.2 将优化之后的结果根据关键词 #HF/6-31G* SCF=tight Test Pop=MK iop(6/33=2) iop(6/42=6) 在Gaussian中重新计算,将输出结果放到Multiwfn中进行RESP计算 (7 /18 / 8 / 1)

    2.3 将优化之后的结果根据您在 (http://bbs.keinsci.com/thread-16190-1-1.html)一文中提到的RESP2的模板进行计算,将输出结果(gas/solv)分别在Multiwfn中计算取平均(7 / 18 / 1)

    2.4 将优化之后的结果放到Antechamber (Amber 20 版本)中进行计算,命令行为:antechamber -i MPD.out -fi gout -o MPD.mol2 -fo mol2 -c resp

3. 结果分析:
屏幕截图 2020-12-01 231512.png


    B.R. 为参考文献通过AM1-BCC算得的结果 (文献:Computer simulations of water flux and salt permeability of the reverse osmosis FT-30 aromatic polyamide membrane)(之前社长说不去计较文献的结果,我想的是如果我算的结果和他相近的话肯定是我的结果更为准确,但是计算结果相差太大了)


    在大量的尝试过程中,不同的初始建模经过几何优化(在gaussian中采用相同的关键词,如1.中所述)之后得到的构象共有三种(还有一种没画图),三种构象的RESP结果均不相同,而且和文献的结果均相差甚远。



想到可能存在问题的地方就是:
1. 在几何优化之前需要得到MPD分子能量最低的几何构象
2. 需要在点群对称上继续处理一下(可我在GaussView中已经将对称性设置成C2v了,而且我尝试过在几何优化时增加symm=loose关键词,结果也没什么改观)




1.png


PS:计算过程中的全部文件我都放到了附件中  (https://pan.baidu.com/s/1V5lAycdyoHz2PLpWlwgh0A) 提取码:qotl  

2万

帖子

25

威望

3万

eV
积分
62470

管理员

公社社长+计算化学玩家

发表于 2020-12-2 11:11:35 | 显示全部楼层
MPD没有多构象问题,直接gview对称化成C2v,拿Gaussian算得到波函数,然后Multiwfn里直接用两步式RESP计算就完了
如果Multiwfn给你的结果不严格满足对称性,加上等价性设置
北京科音自然科学研究中心http://www.keinsci.com  致力于计算化学的发展和传播,长期开办最高水准的各种量子化学、分子动力学与GROMACS、波函数分析与Multiwfn程序等主题的培训,是快速提升计算化学研究水平的最佳途径。欢迎加入“北京科音”公众号获取培训最新消息和计算化学资讯!培训相关信息见《北京科音办的培训班FAQ》
欢迎加入人气最高、水准最高的综合性理论与计算化学交流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!
您需要登录后才可以回帖 登录 | 现在注册!

本版积分规则

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

GMT+8, 2021-1-16 19:18 , Processed in 0.182195 second(s), 27 queries .

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