计算化学公社

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

[其它程序] 使用ztop拟合二面角时曲线很差

[复制链接 Copy URL]

35

帖子

0

威望

235

eV
积分
270

Level 3 能力者

使用高斯进行二面角柔性扫描后,将能量与轨迹文件调整成ztop的输入格式,结果怎么也拟合不好,曲线非常的曲折
有坛友在使用ztop拟合的时候曲线也出现非常崎岖的现象吗
求教有没有别的拟合参数的方法



202311231113074451..png (96.03 KB, 下载次数 Times of downloads: 11)

202311231113074451..png

202311231112435168..png (72.51 KB, 下载次数 Times of downloads: 15)

周期1234

周期1234

截图 2023-11-23 11-13-54.png (61.16 KB, 下载次数 Times of downloads: 15)

周期1,2,4

周期1,2,4

截图 2023-11-23 11-14-36.png (64.4 KB, 下载次数 Times of downloads: 11)

周期24

周期24

截图 2023-11-23 11-14-16.png (58.59 KB, 下载次数 Times of downloads: 17)

周期 234

周期 234

35

帖子

0

威望

235

eV
积分
270

Level 3 能力者

2#
 楼主 Author| 发表于 Post on 2023-11-23 12:35:20 | 只看该作者 Only view this author
当我在top文件中手动更改力常数后,参数曲线也并不是平滑的高斯峰,这是什么原理呀

202311231234408328..png (17.31 KB, 下载次数 Times of downloads: 14)

202311231234408328..png

202311231235182885..png (41.05 KB, 下载次数 Times of downloads: 12)

202311231235182885..png

909

帖子

37

威望

5518

eV
积分
7167

Level 6 (一方通行)

3#
发表于 Post on 2023-11-23 14:35:20 | 只看该作者 Only view this author
渡边京介 发表于 2023-11-23 12:35
当我在top文件中手动更改力常数后,参数曲线也并不是平滑的高斯峰,这是什么原理呀

柔性扫描,导致其他键长键角二面角的改变,可能会使得这些改变对应的分子力学的势能面比较崎岖。你先试试刚性扫描。另外你这PM6明显是有位阻的。gaff2力场通常会高估这个位阻,导致二面角参数必须变得很大去抗衡这个位阻效应。而产生奇怪的势能面。可能的办法是调节产生位阻的原子的原子类型,使得他们之间的vdw排斥变小。

35

帖子

0

威望

235

eV
积分
270

Level 3 能力者

4#
 楼主 Author| 发表于 Post on 2023-11-23 14:57:48 | 只看该作者 Only view this author
ggdh 发表于 2023-11-23 14:35
柔性扫描,导致其他键长键角二面角的改变,可能会使得这些改变对应的分子力学的势能面比较崎岖。你先试试 ...

也就是说initial parameter这条曲线实际上不止有二面角的能量,还通过力场计算了其他的部分,所以会变得很崎岖?

35

帖子

0

威望

235

eV
积分
270

Level 3 能力者

5#
 楼主 Author| 发表于 Post on 2023-11-23 17:09:15 | 只看该作者 Only view this author
换成刚性扫描就拟合的很完美了  谢谢老师

202311231708497764..png (184.81 KB, 下载次数 Times of downloads: 17)

202311231708497764..png

本版积分规则 Credits rule

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

GMT+8, 2026-1-25 17:47 , Processed in 0.740681 second(s), 23 queries , Gzip On.

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