计算化学公社

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

[综合交流] 在AMBER19SB力场中CMAP项的参数含义是什么?

[复制链接 Copy URL]

8

帖子

0

威望

259

eV
积分
267

Level 3 能力者

本帖最后由 hfc 于 2023-6-26 16:20 编辑

1. AMBER在19SB力场中添加了CMAP项,那么在CMAP相共有24*24=576相参数,分别是什么意思?如何使用这些参数进行计算呢?
PS:原始文献截图如下,使用双三次差值,不应该是一个grid有16个参数吗?利用grid四个角的信息求解grid内部的值?所以总共16*24*24个参数。
2. amber当中如何实现给存在耦合的二面角添加CMAP项呢?(非主链)如何指定CMAP涉及的index呢?

202306261617302009..png (43.51 KB, 下载次数 Times of downloads: 16)

202306261617302009..png

6万

帖子

99

威望

5万

eV
积分
124668

管理员

公社社长

2#
发表于 Post on 2023-6-26 18:31:16 | 只看该作者 Only view this author
发帖的时候看清楚板块再发,这在首页的公告栏里已经明确提醒了,求助帖别发到资源分享子版块,这次给你移动了,以后发错地方直接删帖处理。
北京科音自然科学研究中心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

8

帖子

0

威望

259

eV
积分
267

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2023-6-26 20:01:10 | 只看该作者 Only view this author
sobereva 发表于 2023-6-26 18:31
发帖的时候看清楚板块再发,这在首页的公告栏里已经明确提醒了,求助帖别发到资源分享子版块,这次给你移动 ...

不好意思 下次一定注意 sobereva老师

6万

帖子

99

威望

5万

eV
积分
124668

管理员

公社社长

4#
发表于 Post on 2023-6-26 23:08:51 | 只看该作者 Only view this author
这就是把phi,psi二维的校正项势能面用多项式展开,a_ij是i,j阶数的项的展开系数。具体计算是实际程序负责的,不用自己去添加,AMBER里直接用AMBER19SB自然就考虑了cmap了
北京科音自然科学研究中心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

8

帖子

0

威望

259

eV
积分
267

Level 3 能力者

5#
 楼主 Author| 发表于 Post on 2023-6-27 00:13:59 | 只看该作者 Only view this author
sobereva 发表于 2023-6-26 23:08
这就是把phi,psi二维的校正项势能面用多项式展开,a_ij是i,j阶数的项的展开系数。具体计算是实际程序负责的 ...

那老师请问TOP文件中的CMAP参数是24*24 每个具体的参数是意思呢 和这个aij 对应不起来?
第二 我想要不仅仅是在主链上保留CMAP势 我还想要在侧链上也添加CMAP 请问在amber中如何实现?

1664

帖子

5

威望

4768

eV
积分
6532

Level 6 (一方通行)

喵星人

6#
发表于 Post on 2023-6-27 01:29:26 | 只看该作者 Only view this author
本帖最后由 喵星大佬 于 2023-6-27 01:36 编辑


cmap对应的是两个二面角的耦合项,从第一个点开始就是24×24的点阵,但是由于二面角是周期性的,所以转了一圈之后要回到原点,即头=尾

cmap的参数是势能点阵,但是实际计算中角度不可能正对着某个点,所以使用的时候程序是先通过带约束条件的双三次插值得到完整的修正是能面(即让每个点之间空隙上也能计算数值),再在计算中使用

8

帖子

0

威望

259

eV
积分
267

Level 3 能力者

7#
 楼主 Author| 发表于 Post on 2023-6-27 09:33:11 | 只看该作者 Only view this author
本帖最后由 hfc 于 2023-6-27 09:34 编辑
喵星大佬 发表于 2023-6-27 01:29
cmap对应的是两个二面角的耦合项,从第一个点开始就是24×24的点阵,但是由于二面角是周期性的,所以转了 ...

谢谢老师我了解CMAP的含义了 这个应该是gromacs下面的TOP文件 请问amber中如何实现自定义CMAP呢?

140

帖子

0

威望

858

eV
积分
998

Level 4 (黑子)

8#
发表于 Post on 2025-11-20 04:31:29 | 只看该作者 Only view this author
sobereva 发表于 2023-6-26 23:08
这就是把phi,psi二维的校正项势能面用多项式展开,a_ij是i,j阶数的项的展开系数。具体计算是实际程序负责的 ...

就是因为ff19SB多了这个CMAP项,导致使用Amber自带的ParmEd构造TI的初始结构的命令"TiMerge"报错,使用ff14SB力场则不会有这样的问题。在github上最早2021年就有人注意到了这个问题,但是直到Amber24版本,这个问题似乎仍未解决,目前如果想要Amber跑TI依然无法使用ff19SB力场。

1664

帖子

5

威望

4768

eV
积分
6532

Level 6 (一方通行)

喵星人

9#
发表于 Post on 2025-11-21 02:21:30 | 只看该作者 Only view this author
ZetaFunction 发表于 2025-11-20 04:31
就是因为ff19SB多了这个CMAP项,导致使用Amber自带的ParmEd构造TI的初始结构的命令"TiMerge"报错,使用ff ...

我已经给了gmx直接可以用的Amber19SB力场包,TI直接用gmx跑然后用那个python包算就行

本版积分规则 Credits rule

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

GMT+8, 2026-1-23 20:34 , Processed in 0.235100 second(s), 23 queries , Gzip On.

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