| 查看: 4948 | 回复: 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 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 动力学强帖 | 气态分子化学动力学 |
» 猜你喜欢
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有8人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有3人回复
-
明天应该可查了！？
已经有5人回复
-
如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》
已经有7人回复
-
没有任何消息-是不是就凉了
已经有8人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
人气不行了
已经有10人回复
-
我面上完蛋了
已经有10人回复
-
能否退出参与的面上项目解除限项
已经有24人回复
-
filecode，4个jtjc了
已经有17人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 如何用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人回复
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
2楼2010-06-12 21:39:48
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
20