计算化学公社

 找回密码 Forget password
 注册 Register
Views: 15917|回复 Reply: 17
打印 Print 上一主题 Last thread 下一主题 Next thread

[综合交流] EPON862单体oplsaa力场选择问题

[复制链接 Copy URL]

155

帖子

0

威望

2879

eV
积分
3034

Level 5 (御坂)

跳转到指定楼层 Go to specific reply
楼主
EPON862化学结构如图所示;对其选择oplsaa力场时候遇到几个问题:
1、与苯环相连的氧原子给什么type比较合适,另外与环氧基团相连的两个C原子选择那种type比较合适;
2、所有原子type确定之后,很有可能总电荷不为0,该如何处理
3、还有就是这个oplsaa力场文件的前面给了部分oplsua的力场type,能不能混合来用。
这是oplsaa力场文件的位置:http://dasher.wustl.edu/tinker/distribution/params/oplsaa.prm



Screenshot from 2016-10-14 13-35-20.png (63.04 KB, 下载次数 Times of downloads: 57)

EPON862化学结构

EPON862化学结构

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

2#
发表于 Post on 2016-10-14 18:30:48 | 只看该作者 Only view this author
可以尝试借助http://www.aribeiro.net.br/mktop/http://erg.biophys.msu.ru/tpp/来确定原子类型
原子电荷一般是需要自己计算的,加和肯定为0。用RESP、AM1-BCC、ADCH这些对静电势重现性好的为宜。
没太大必要的话还是用OPLS-AA为宜。UA也省不了多少时间。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

155

帖子

0

威望

2879

eV
积分
3034

Level 5 (御坂)

3#
 楼主 Author| 发表于 Post on 2016-10-14 20:26:26 | 只看该作者 Only view this author
sobereva 发表于 2016-10-14 18:30
可以尝试借助http://www.aribeiro.net.br/mktop/和http://erg.biophys.msu.ru/tpp/来确定原子类型
原子电 ...

谢谢Sob老师,我看到96年关于oplsaa力场描述的文章《Development and Testing of the OPLS All-Atom Force Field on Conformational Energetics and Properties of Organic Liquids》中提到——That is, if a set of partial charges had to be modified owing to unsatisfactory reproduction of liquid properties, the associated torsional terms had to be refit owing to changes in the other terms in eq 5, especially the nonbonded energy.
问题就在于不使用他们的charge,而使用自己计算得到charge后,torsional参数还需要修正?

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

4#
发表于 Post on 2016-10-14 20:54:52 | 只看该作者 Only view this author
dark_cosmos 发表于 2016-10-14 20:26
谢谢Sob老师,我看到96年关于oplsaa力场描述的文章《Development and Testing of the OPLS All-Atom Forc ...


原子电荷是显著依赖于化学环境的,同一个原子类型在不同分子里的电荷也是不同的。
OPLSAA弄各种参数,也包括电荷,目的是尽可能好重现实验数据,所以这些参数和电荷之间难免有一定耦合。但如果照它们这么弄,对付一个新的分子还得自己将模拟数据和实验反复对比调整,显然不现实。所以相对简单的做法就是用一个对静电势重现性不错的原子电荷计算方法。或者若OPLSAA库里已有的一些分子和当前体系有类似性,可以将原子电荷挪到当前体系,让后适当人为处理一下使净电荷为0。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

155

帖子

0

威望

2879

eV
积分
3034

Level 5 (御坂)

5#
 楼主 Author| 发表于 Post on 2016-10-15 15:29:25 | 只看该作者 Only view this author
sobereva 发表于 2016-10-14 20:54
原子电荷是显著依赖于化学环境的,同一个原子类型在不同分子里的电荷也是不同的。
OPLSAA弄各种参数, ...

谢谢Sob老师,TPPmktop给的atom type给了我很大的参考。MKTOP给的atom type不太科学。最终加在一块电荷为0.03e,我觉得已经很好了,打算先计算看看结果如何。
请问sob老师,您所说的“然后适当人为处理一下使净电荷为0”具体怎么做呢,对于这个分子差的0.03e能不能给处理掉

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

6#
发表于 Post on 2016-10-15 16:58:22 | 只看该作者 Only view this author
dark_cosmos 发表于 2016-10-15 15:29
谢谢Sob老师,TPPmktop给的atom type给了我很大的参考。MKTOP给的atom type不太科学。最终加在一块电荷为 ...

这点电荷不怎么会影响结果,直接平摊到几个原子上抹掉就完了
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

155

帖子

0

威望

2879

eV
积分
3034

Level 5 (御坂)

7#
 楼主 Author| 发表于 Post on 2016-10-15 18:36:42 | 只看该作者 Only view this author
sobereva 发表于 2016-10-15 16:58
这点电荷不怎么会影响结果,直接平摊到几个原子上抹掉就完了

好的,谢谢sob老师

155

帖子

0

威望

2879

eV
积分
3034

Level 5 (御坂)

8#
 楼主 Author| 发表于 Post on 2016-11-1 22:50:25 | 只看该作者 Only view this author
sobereva 发表于 2016-10-14 20:54
原子电荷是显著依赖于化学环境的,同一个原子类型在不同分子里的电荷也是不同的。
OPLSAA弄各种参数, ...

请问sob老师,为什么OPLSAA力场的K4为0,也就是力场文件中(Ryckaert-Bellemans形式)最后两项为0。而且始终为0。这样要这项的作用在哪
还有就是看到oplsaa的文献提到“The charges for functional groups are taken to be transferable between molecules and the use of neutral subunits makes the derivation of charges for
large molecules straightforward. This represents a major difference with the AMBER94 force field 3c for which charges are obtained on a case-by-case basis from fitting to electrostatic potential surfaces from ab initio 6-31G* calculations.”。后面改进的力场参数有使用CM1A-BCC charge mode的,这就有点郁闷了,是不是只能利用oplsaa力场自身的charge了,计算得到的resp或者时AM1-BCC电荷用起来就不太合适了。谢谢老师

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

9#
发表于 Post on 2016-11-2 00:13:16 | 只看该作者 Only view this author
dark_cosmos 发表于 2016-11-1 22:50
请问sob老师,为什么OPLSAA力场的K4为0,也就是力场文件中(Ryckaert-Bellemans形式)最后两项为0。而且 ...

后两项始终为0就为0,说明当前情况用不着。但留出那两项可以为实现特殊目的之用
如果OPLSAA分子库里有能借到你的体系来的片段的电荷定义可以借过来,如果没有合适的能借来的就只能自己直接算。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

155

帖子

0

威望

2879

eV
积分
3034

Level 5 (御坂)

10#
 楼主 Author| 发表于 Post on 2016-11-2 09:51:00 | 只看该作者 Only view this author
sobereva 发表于 2016-11-2 00:13
后两项始终为0就为0,说明当前情况用不着。但留出那两项可以为实现特殊目的之用
如果OPLSAA分子库里有能 ...

谢谢sob老师,我现在计算的另外一个分子,电荷找不全。计算的resp电荷又与已有的电荷差别太大(甚至是符号上的差别),请问sob老师这种情况下该如何处理。

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

11#
发表于 Post on 2016-11-2 11:10:14 | 只看该作者 Only view this author
dark_cosmos 发表于 2016-11-2 09:51
谢谢sob老师,我现在计算的另外一个分子,电荷找不全。计算的resp电荷又与已有的电荷差别太大(甚至是符 ...

如果RESP电荷符合化学直觉,就用这个
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

155

帖子

0

威望

2879

eV
积分
3034

Level 5 (御坂)

12#
 楼主 Author| 发表于 Post on 2016-11-2 14:25:23 | 只看该作者 Only view this author
sobereva 发表于 2016-11-2 11:10
如果RESP电荷符合化学直觉,就用这个

谢谢sob老师,讲道理是更加相信resp的。但是opls原始文献里面说他们的很多都是精心调整过的,比如甲基乙基里面的的C和H,这种就与具体分子里面计算得到的resp差别很大了。以前没有什么这方面的经验,所以不太清楚究竟该如何处理这种情况。

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

13#
发表于 Post on 2019-5-15 16:14:18 | 只看该作者 Only view this author
dark_cosmos 发表于 2016-11-2 14:25
谢谢sob老师,讲道理是更加相信resp的。但是opls原始文献里面说他们的很多都是精心调整过的,比如甲基乙 ...

在J. Phys. Chem. B, 121, 3864 (2017)中证明了CM5电荷乘以1.2很适合结合OPLS-AA力场模拟,所以建议用1.2*CM5电荷。CM5电荷可以用Multiwfn非常容易地计算,载入诸如Gaussian的fch文件,进入布居分析功能(主功能7)后选CM5即可。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

155

帖子

0

威望

2879

eV
积分
3034

Level 5 (御坂)

14#
 楼主 Author| 发表于 Post on 2019-5-16 09:54:46 | 只看该作者 Only view this author
sobereva 发表于 2019-5-15 16:14
在J. Phys. Chem. B, 121, 3864 (2017)中证明了CM5电荷乘以1.2很适合结合OPLS-AA力场模拟,所以建议用1.2 ...

感谢sob老师,去年参加过gromacs的培训,目前对于这些力场电荷算是搞清楚了

12

帖子

0

威望

1497

eV
积分
1509

Level 5 (御坂)

15#
发表于 Post on 2020-8-23 21:20:52 | 只看该作者 Only view this author
dark_cosmos 发表于 2016-10-15 15:29
谢谢Sob老师,TPPmktop给的atom type给了我很大的参考。MKTOP给的atom type不太科学。最终加在一块电荷为 ...

请问你现在还能正常使用TPPmktop吗,我submit后,显示这个:Calculation was run by you or another user. Starting two simultaneous calculations is forbidden. Wait a little and update page - your request wil be restarted!

本版积分规则 Credits rule

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

GMT+8, 2026-2-21 07:00 , Processed in 0.185777 second(s), 23 queries , Gzip On.

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