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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » LAMMPS/DL_POLY » 求助：如何选取morse势截断半径
41/1返回列表
查看: 2177  |  回复: 3
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
【悬赏金币】回答本帖问题，作者半寸日光将赠送您 20 个金币

半寸日光

新虫 (小有名气)

[求助] 求助：如何选取morse势截断半径

lammps模拟界面问题，界面之间作用力用morse势，查了文献，只查到morse势的三个基本参数，没找到cutoff，这个cutoff选取有什么经验吗？

@smutao @oxox6085 发自小木虫Android客户端
回复此楼

» 猜你喜欢
1楼 2021-06-21 17:11:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

HJieY

新虫 (初入文坛)
前辈
      您好！我目前在用分子动力学做单晶铜的纳米切削模拟，在计算原子间相互作用力时，我也用到了Morse势函数，但不知道如何选取截断半径，而且我经过测试发现截断半斤的选取对于原子间的相互作用能的计算确实是非常重要的，但不知道您现在有什么好的方法了吗？
期待您的回复！
                                                                                                                                                                                                                                                         模拟瞌睡虫
赞一下 回复此楼
人在塔在
2楼2021-11-17 10:37:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

半寸日光

新虫 (小有名气)
引用回帖:
2楼: Originally posted by HJieY at 2021-11-17 10:37:40
前辈
      您好！我目前在用分子动力学做单晶铜的纳米切削模拟，在计算原子间相互作用力时，我也用到了Morse势函数，但不知道如何选取截断半径，而且我经过测试发现截断半斤的选取对于原子间的相互作用能的计算确 ...

我之前请教过别人，得到的有两个回复一是Morse势函数是二元的，可以根据公式画图的，用Matlab画出来就知道怎么选了，随便选，趋近于0的数，越大越准计算也越慢；二是没什么特别要求，12A左右

发自小木虫Android客户端
赞一下 回复此楼
3楼2021-11-17 20:57:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

乐奏球

新虫 (著名写手)
引用回帖:
3楼: Originally posted by 半寸日光 at 2021-11-17 20:57:59
我之前请教过别人，得到的有两个回复一是Morse势函数是二元的，可以根据公式画图的，用Matlab画出来就知道怎么选了，随便选，趋近于0的数，越大越准计算也越慢；二是没什么特别要求，12A左右
...

你好，按照你说的截断半径越大越好的说法，我试了一下，但是不符合，我的截断半径只在某些值时系统是正常的，截断半径过大或过小系统都会不正常，是因为参数有问题吗
赞一下 回复此楼
4楼2021-12-22 17:44:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 半寸日光 的主题更新
41/1返回列表
不应助 确定回帖应助 (注意：应助才可能被奖励，但不允许灌水，必须填写15个字符以上)
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[博后之家] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +4 2JOx3r2CYEgw 2026-08-22 6/300 2026-08-23 19:34 by lXgU0RvNTDC9
[博后之家] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 2JOx3r2CYEgw 2026-08-22 6/300 2026-08-23 19:31 by lXgU0RvNTDC9
[公派出国] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 2JOx3r2CYEgw 2026-08-22 4/200 2026-08-23 19:22 by lXgU0RvNTDC9
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 2JOx3r2CYEgw 2026-08-21 5/250 2026-08-23 19:10 by lXgU0RvNTDC9
[论文投稿] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-21 8/400 2026-08-23 18:58 by lXgU0RvNTDC9
[论文投稿] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +4 2JOx3r2CYEgw 2026-08-21 9/450 2026-08-23 18:46 by lXgU0RvNTDC9
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 2JOx3r2CYEgw 2026-08-21 5/250 2026-08-23 14:33 by KM5EcsNQRBPn
[教师之家] 跳槽后在研项目怎么办？ +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 今天放榜吗？ +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[博后之家] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 DpSrDtM079iu 2026-08-22 5/250 2026-08-23 04:07 by OEbVnUOu01ol
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 2JOx3r2CYEgw 2026-08-22 7/350 2026-08-23 03:55 by OEbVnUOu01ol
[硕博家园] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 7/350 2026-08-23 03:43 by OEbVnUOu01ol
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +4 QTy3jDtz1uLt 2026-08-21 11/550 2026-08-23 02:55 by OEbVnUOu01ol
[考博] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 DpSrDtM079iu 2026-08-22 3/150 2026-08-23 00:55 by OEbVnUOu01ol
[教师之家] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +4 2JOx3r2CYEgw 2026-08-21 4/200 2026-08-22 16:52 by sunzitan
[基金申请] 放榜前的不淡定 40+4 snowwithsea 2026-08-19 14/700 2026-08-21 23:51 by cratir
[基金申请] 科研孤儿太难了 +17 我4大白菜 2026-08-20 18/900 2026-08-21 20:57 by zhangev
[基金申请] 应该是下周三26日公布了吧？ +4 哈哈蛤？ 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 重要消息，中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
信息提示
关闭
请填处理意见
关闭