| 查看: 4943 | 回复: 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 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 动力学强帖 | 气态分子化学动力学 |
» 猜你喜欢
filecode，4个jtjc了
已经有17人回复
-
号外号外
已经有5人回复
-
能否退出参与的面上项目解除限项
已经有22人回复
-
人气不行了
已经有9人回复
-
没有任何消息-是不是就凉了
已经有6人回复
-
如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》
已经有5人回复
-
2026国自然函评费到账
已经有19人回复
-
建议基金发布提前给出明确的时间点
已经有16人回复
-
今日不放榜？网传国自然预计 8 月 27 日可查结果
已经有21人回复
-
估计是周四
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 如何用Matlab曲线拟合一个多输入参数，而且其中一个输入参数与另一个参数有关呢？
已经有6人回复
- 含双键二面角，柔性扫描势能曲线，计算到一半就出错（gaussian）
已经有8人回复
- 求助，如何设置D-P模型
已经有3人回复
- 急求解答，关于力场的电荷参数！！！
已经有4人回复
- E-V曲线拟合后怎么找到能量最低点对应的V值
已经有4人回复
- 复杂网络参数图
已经有6人回复
- 用高斯怎么得到下图中激发态的所有信息（主要是具体的轨道能级数值）
已经有3人回复
- 如何利用GULP拟合分子间势能曲线得出力场参数
已经有10人回复
- 两个自变量怎么拟合求参数值呢？
已经有7人回复
- 求助original曲线拟合。
已经有6人回复
- 紧急求助高斯扫势能曲线（同分子中不同基团之间的距离如何定义）
已经有4人回复
- 求助参数拟合~~~着急
已经有16人回复
- 【求助】氢分子离子的势能曲线
已经有3人回复
- LIBSVM中的svm-train命令后面的-wi参数如何使用？
已经有4人回复
- 非线性最小二乘法编程问题
已经有8人回复
- 用origin里非线性最小二乘法是怎么拟合参数的？？？？
已经有14人回复
- ghcacj与gtolv8688合力推出Code Music之FAQ
已经有53人回复
- 求助非线性最小二乘法
已经有5人回复
- 【求助】分子的LJ势能参数的拟合
已经有8人回复
- 【求助】如何通过理论计算得到分子各态能级以及势能曲线或曲面？
已经有4人回复
- 【求助】蛋白质近表面静电势能模拟
已经有3人回复
- 【求助】使用MATLAB怎么实现拟合力场参数的程序？【已解决】
已经有15人回复
- 【求助】最小二乘法求参数
已经有6人回复
2楼2010-06-12 21:39:48
bay__gulf
金虫 (著名写手)
刘苏州
- 模拟EPI: 8
- 应助: 9 (幼儿园)
- 贵宾: 4.85
- 金币: 2332.8
- 红花: 1
- 帖子: 1344
- 在线: 271小时
- 虫号: 592012
- 注册: 2008-09-03
- 专业: 理论和计算化学
- 管辖: 分子模拟
★ ★
xiaowu787(金币+3): 2010-06-13 09:32:47
xiaowu787(金币+3): 2010-06-13 09:34:02
lei0736(金币+2):谢谢 2010-06-13 09:36:58
xiaowu787(金币+3):Thanks 2010-06-13 10:30:42
xiaowu787(金币+3): 2010-06-13 09:32:47
xiaowu787(金币+3): 2010-06-13 09:34:02
lei0736(金币+2):谢谢 2010-06-13 09:36:58
xiaowu787(金币+3):Thanks 2010-06-13 10:30:42
|
非线性拟合问题, 很多现成软件可以做, 推荐1stOpt http://www.mdbbs.org/viewthread.php?tid=12576 |
3楼2010-06-12 21:56:19
4楼2010-06-12 22:27:17
|
function [a1,b1,a2,b2]=forcefieldfitting(r,e) %使用DFT计算的数据拟合力场参数,e是使用量化计算出现的相互作用能，r是两个分子片段之间的距离 %LJ-12-06公式中的EPSILON(C):a1 %LJ-12-06公式中的SIGMA(C):b1 %LJ-12-06公式中的EPSILON(O):a2 %LJ-12-06公式中的SIGMA(O):b2 %单个的CO2的能量可以认为是两个O的能量和一个C的能量的加和 if(length(e)==length(r)) n=length(e); else disp('e和r的维数不相等！'); return; end %维数检查 A=zeros(4,3) B=zeros(4,2); for i=1:n A(1,1)=A(1,1)+(384*a1^2*b1^23)/r(i)^24; A(1,2)=A(1,2)+(192*a1^2*b1^11)/r(i)^12; A(1,3)=A(1,3)+(48*a1*e*b1^5)/r(i)^6; A(2,1)=A(2,1)+(32*a1*b1^24)/r(i)^24; A(2,2)=A(2,2)+(32*a1*b1^12)/r(i)^12 A(2,3)=A(2,3)+(8*e(i)*b1^6)/r(i)^6; A(3,1)=A(3,1)+(1536*a2^2*b2^23)/r(i)^24; A(3,2)=A(3,2)+(768*a2^2*b2^23)/r(i)^12; A(3,3)=A(3,3)+(96*a2*e(i)*b2^5)/r(i)^6; A(4,1)=A(4,1)+(128*a2*b2^24)/r(i)^24; A(4,2)=A(4,2)+(128*a2*b2^12)/r(i)^12; A(4,3)=A(4,3)+(16*e(i)*b2^6)/r(i)^6; B(1,1)=B(1,1)+(576*a1^2*b1^17)/r(i)^18; B(1,2)=B(1,2)+(96*a1*e(i)*b1^11)/r(i)^12; B(2,1)=B(2,1)+(64*a1*b1^18)/r(i)^18; B(2,2)=B(2,2)+(8*e(i)*a1^12)/r(i)^12; B(3,1)=B(3,1)+(2304*a2^2*b2^17)/r(i)^18; B(3,2)=B(3,2)+(192*a2*e(i)*b2^11)/r(i)^12; B(4,1)=B(4,1)+(256*a2*b2^18)/r(i)^18; B(4,3)=B(4,3)+(16*e(i)*b2^12)/r(i)^12; end s=A\B; a1=s(1); b1=s(2); a2=s(3); b2=s(4); 我这个预先设置UFF的参数，然后这么拟合出来的，不知道行不行？ |
5楼2010-06-12 22:32:50
xiaowu787
木虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 1406.9
- 散金: 95
- 红花: 1
- 帖子: 550
- 在线: 430.4小时
- 虫号: 326816
- 注册: 2007-03-18
- 性别: GG
★ ★ ★
雪狼乖乖(金币+3):谢谢您的参与 2010-06-15 07:48:08
雪狼乖乖(金币+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 是这样吗？请高手指点一下，毕竟对于多参数的函数拟合本身就很麻烦，自己的思路再错了，那就更没指望了。希望有过拟合经验的指点指点。 |
6楼2010-06-14 15:42:29
雪狼乖乖 
木虫 (正式写手)
- 模拟EPI: 1
- 应助: 18 (小学生)
- 贵宾: 0.005
- 金币: 2805.8
- 散金: 567
- 红花: 34
- 帖子: 901
- 在线: 402.2小时
- 虫号: 532524
- 注册: 2008-03-25
- 性别: GG
- 专业: 理论和计算化学
★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
lei0736(金币+5):谢谢 2010-06-15 19:52:35
小木虫(金币+0.5):给个红包，谢谢回帖交流
lei0736(金币+5):谢谢 2010-06-15 19:52:35
|
1stOpt 是国人开发的软件 效果非常好 使用也极其简单 推荐使用 另外 就该势能立场参数的拟合 我个人觉得精度并不是关键 而是物理意义的阐述 否则拟合的精度再高 又有什么用呢？ 推荐您考虑一下物理化学里 Van der Wals方程的建立与改进过程 可以先从一个简单模型入手 不断修改表达的物理意义 以期达到精度和物理意义都统一 拟合的软件很多很多 origin MATLAB mathmatic 都可以实现 关键在初始值的给出和拟合算法的设置 一般的小技巧是 先用遗传算法算出精确解的区域范围 然后再用爬坡法或其他较精确的计算方法计算精确拟合值 个人经验 仅供参考  该图是立体的势能扫描图 其波谷态势可以大致说明势能曲线 请注意波谷值并不是唯一 整个立场中 有多个波谷值 原因是分子的结构并不固定 所以你计算出的势能曲线可能并不唯一 是众多势能曲线的一个 所以不需要过多考虑精度问题 能说明其物理意义即可 [ Last edited by 雪狼乖乖 on 2010-6-15 at 08:13 ] |
7楼2010-06-15 08:04:07
8楼2010-06-15 10:32:36
xiaowu787
木虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 1406.9
- 散金: 95
- 红花: 1
- 帖子: 550
- 在线: 430.4小时
- 虫号: 326816
- 注册: 2007-03-18
- 性别: GG
9楼2010-06-15 12:42:41
戴我飞翔
铁杆木虫 (职业作家)
金虫
- 应助: 0 (幼儿园)
- 金币: 21968.9
- 沙发: 1
- 帖子: 3511
- 在线: 246.5小时
- 虫号: 1031645
- 注册: 2010-05-30
- 专业: 数学物理
10楼2010-06-15 13:13:32
20