应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » LAMMPS/DL_POLY » 【求助】如何根据势能曲线拟合力场参数
51/1返回列表
查看: 4952  |  回复: 19
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
本帖产生 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 ]
回复此楼

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

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

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2010-06-12 21:26:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

雪狼乖乖

木虫 (正式写手)
★ ★ ★ ★ ★ ★
小木虫(金币+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 ]
赞一下(2人) 回复此楼
7楼2010-06-15 08:04:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 20 个回答

zyj8119

木虫 (著名写手)
xiaowu787(金币+1): 2010-06-13 09:32:30
这个也许得编程，我也同问高人。。
赞一下 回复此楼
好好学习，天天向上。
2楼2010-06-12 21:39:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bay__gulf

金虫 (著名写手)

刘苏州
★ ★
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
赞一下(1人) 回复此楼
3楼2010-06-12 21:56:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zyj8119

木虫 (著名写手)
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
雪狼乖乖(金币+5):感谢您的参与 2010-06-15 07:46:37
雪狼乖乖(金币+5):感谢您的参与 2010-06-15 07:47:00
雪狼乖乖(金币+2, 模拟EPI+1):您的计算平台是不是MATLAB？您的算法收敛否？精度如何？能否给予说明？ 2010-06-15 07:52:35
不知道GULP能做这个不？
赞一下(3人) 回复此楼
好好学习，天天向上。
4楼2010-06-12 22:27:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 20 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 国合现在查不到了吗？ +8 chengyan1220 2026-08-24 16/800 2026-08-26 08:50 by paperzjh
[基金申请] 牛来！米来！面来！ +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 7K1CJE38xLG4 2026-08-25 5/250 2026-08-26 07:17 by cNXvBfCpiZOM
[论文投稿] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 7K1CJE38xLG4 2026-08-25 4/200 2026-08-26 02:04 by cNXvBfCpiZOM
[基金申请] 在坚冰还盖着北海的时候，我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 明天应该可查了！？ +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 有没有大神帮我看看基金代码 28+4 1234567wang 2026-08-24 10/500 2026-08-25 19:15 by lfy8008
[基金申请] 放榜前的不淡定 20+4 snowwithsea 2026-08-19 19/950 2026-08-25 17:57 by chick875
[基金申请] 某些机构，以效率低为荣，以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 今日不放榜？网传国自然预计 8 月 27 日可查结果 +17 医学老男孩 2026-08-20 22/1100 2026-08-25 15:36 by 医学老男孩
[基金申请] 如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[基金申请] 人气不行了 +11 fansofjerry 2026-08-21 11/550 2026-08-25 11:04 by 孤独的英雄6
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +5 process2012 2026-08-23 7/350 2026-08-25 09:42 by huixian257
[基金申请] 能否退出参与的面上项目解除限项 +21 koalala 2026-08-24 24/1200 2026-08-24 19:25 by 家与远方
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[基金申请] 让我中一个面上吧！ +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 by 太傻了
[教师之家] 跳槽后在研项目怎么办？ +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 时间戳今天，20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 看来今天不会放榜了？ +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
信息提示
关闭
请填处理意见
关闭