24小时热门版块排行榜    

查看: 4243  |  回复: 19
本帖产生 1 个 模拟EPI ,点击这里进行查看
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

xiaowu787

木虫 (正式写手)

[交流] 【求助】如何根据势能曲线拟合力场参数已有8人参与

下图是Au-Au的Dimer二聚体的解离曲线


势能曲线有了,现在需要拟合下面三个公式里的多个参数,具体该怎么做?有直接能用的软件吗?请大家帮忙
E=E1+E2+E3


参数a,b,c,d,e,f


参数p,q,s


参数g,h,i,j,k,m


谢谢帮忙!!!

[ Last edited by xiaowu787 on 2010-6-12 at 21:31 ]
回复此楼

» 收录本帖的淘帖专辑推荐

动力学强帖 气态分子化学动力学

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaowu787

木虫 (正式写手)

★ ★ ★
雪狼乖乖(金币+3):谢谢您的参与 2010-06-15 07:48:08
如果势函数形式比较复杂
以求角度部分参数为例:
E(total)=E(bond)+E(angle)+E(nobond)
用第一性原理扫描出角度变化的势能曲线,与角度变化相关的项有E(total),E(angle),E(nobond)三项,拟合的时候是不是利用
E(total)2-E(total)1=[E(angle)+E(nobond)]2-[E(angle)+E(nobond)]1
是这样吗?请高手指点一下,毕竟对于多参数的函数拟合本身就很麻烦,自己的思路再错了,那就更没指望了。希望有过拟合经验的指点指点。
引用回帖:
Originally posted by bay__gulf at 2010-06-12 21:56:19:
非线性拟合问题, 很多现成软件可以做, 推荐1stOpt
http://www.mdbbs.org/viewthread.php?tid=12576

6楼2010-06-14 15:42:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaowu787

木虫 (正式写手)

我现在要用的是ReaxFF力场,发现里面的许多参数的物理意义并不是十分明确,因此许多力场数据使用起来确实不踏实。不过这个力场不着重阐述参数的物理意义,他们注重的势函数的意义计算结果的描述。现在刚开始学习,感觉很多东西都对不上号。
引用回帖:
Originally posted by 雪狼乖乖 at 2010-06-15 08:04:07:
1stOpt 是国人开发的软件 效果非常好 使用也极其简单 推荐使用
另外 就该势能立场参数的拟合 我个人觉得精度并不是关键 而是物理意义的阐述 否则拟合的精度再高 又有什么用呢?
推荐您考虑一下物理化学里 Van d ...

9楼2010-06-15 12:42:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaowu787

木虫 (正式写手)

引用回帖:
Originally posted by 老虎大王 at 2010-06-15 16:17:44:
既然有数据,有公式,就直接找软件拟合就是了。你有15个参数,数据点要搞多一点才行。

现在正在尝试,感觉参数比较多,总不理想。我现在做的是别人做过的体系,希望能把他们的参数重现一次。但是把三部分加在一起时,由于函数形式比较复杂在lstOpt里面直接不能识别,我想把一部分参数值先代进去一部分一部分拟合,不知这样的思路对不?
12楼2010-06-15 22:50:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaowu787

木虫 (正式写手)

引用回帖:
Originally posted by 老虎大王 at 2010-06-16 08:31:41:
你可以根据经验,先固定一部分参数,来拟合另一部分参数。再在此基础上,换一些固定与非固定的参数,反复拟合,直到目标函数达到足够的精度。当然了,这样复杂的函数,拟合起来一定是很困难的,这个要有充分的思想 ...

谢谢您的建议和鼓励!
14楼2010-06-16 08:52:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaowu787

木虫 (正式写手)

引用回帖:
Originally posted by zyj8119 at 2010-06-18 12:15:20:
gulp拟合的还行,然后还有 孙淮教授他们组里面编写的一个 工具 direct force field

不同的力场,拟合出来的参数很多情况下是不能通用的吧,毕竟势函数形式一般都不一样,里面参数的物理意义可能也不同呀
19楼2010-06-18 14:55:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xiaowu787 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见