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

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » 其他 » 【讨论】分子模拟中力场参数的调节讨论
51/1返回列表
查看: 2674  |  回复: 7
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

abin007

银虫 (初入文坛)

[交流] 【讨论】分子模拟中力场参数的调节讨论已有3人参与

经验力场在分子模拟中的重要性不用我多说了，但困惑我们的是，当我们的结构在现有力场中找不到参数时，或者是即使找到了力场参数，模拟结果并不准确，这时候调整立场参数成了不可避免的问题。
    问题的复杂性在于，这么多的力场参数，从哪儿下手？怎么调整？而且有时候是牵一发而动全身啊。我最近涉及到了这一个问题，就我的经验跟大家交流交流，希望得到大家的响应。
    就我的例子而言，如果想要修改二面角参数，首先第一步是取消二面角中1-4原子间的静电力和范德华力，可以这样比方：1-4原子间的非键力像是一个大小可变的盒子的尺寸，如果盒子里面的系统始终小于盒子，那么该盒子大小的变化不会带来多大的系统变化；一旦盒子里面的系统要主动变大时（就像我们调整二面角的参数时），这时该盒子的尺寸就会发挥较明显的限制作用了。因此，如果我们希望调整二面角的力场参数并能看到明显地效果时，应该取消1-4原子间的非键力。不知我这个比喻是否恰当。
   接下来，我的办法就比较笨了，逐一调整二面角的力参数，来比较最终结果，有时候会添加或删除一些二面角参数。我的原则是，二面角的个数从少到多，不要刚开始就一股脑儿地把所有的二面角都用上。
   以上就是我跟大家分享的浅薄经验，希望能跟大家继续交流这个话题。
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐，对您同样有帮助:
1楼2011-02-13 19:00:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaowu759

铁杆木虫 (著名写手)
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
ghcacj(金币+2): 谢谢 2011-02-16 13:04:48
我的理解是，做微小范围的调节是允许的，否则需要重新参数化
赞一下(2人) 回复此楼
5楼2011-02-15 09:59:23
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答

043114076

木虫 (正式写手)
★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
ghcacj(金币+3): 谢谢 2011-02-16 13:04:28
不太理解你那个盒子的比喻，能不能再详细解释一下。
你说的取消1-4原子间的静电力和范德华力，据我了解，很多力场本身就将1-4原子间的非键作用定义为0。不明白你说的取消的意思
赞一下(2人) 回复此楼
2楼2011-02-14 10:26:04
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

abin007

银虫 (初入文坛)
★ ★ ★ ★
zh1987hs(金币+4): 谢谢 2011-02-18 22:53:52
引用回帖:
Originally posted by 043114076 at 2011-02-14 10:26:04:
不太理解你那个盒子的比喻，能不能再详细解释一下。
你说的取消1-4原子间的静电力和范德华力，据我了解，很多力场本身就将1-4原子间的非键作用定义为0。不明白你说的取消的意思

谢谢回复，是的，这可能跟我用的软件有关，在NAMD中使用charmm力场时，是可以调节1-4原子间的非键力的，比如说对1-4原子间的静电力的调节，可以设定一个0-1之间的标定参数，我就是在设定这个标定参数和调节二面角的力场参数上卡壳了，我刚开始使用的标定参数为1，就是完全考虑静电力，在这样的状况下怎么调节二面角的力场参数都不能达到理想值，当我取消1-4原子间的所有非键力时，发现结构对二面角力场参数的响应变得比较明显了。当找到最佳的力场参数时，再加上1-4原子间的非键力，结构依然保持较好。我上述的盒子的比喻不太成熟，大家不用在意也罢。
赞一下(1人) 回复此楼
3楼2011-02-14 21:40:25
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

043114076

木虫 (正式写手)
★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
ghcacj(金币+3): 谢谢 2011-02-16 13:04:38
我上面犯了一个错误。很多力场将1-4原子间的非键作用因子定义为0.5 而不是0。
1-4原子间的非键作用因子应该从力场文献中找，而不是自己去优化。力场在拟合时已经规定了1-4原子间的非键作用因子。
赞一下(2人) 回复此楼
4楼2011-02-15 09:21:12
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答
普通表情 龙 兔 虎 猫 高级回复(可上传附件)
信息提示
关闭
请填处理意见
关闭