24小时热门版块排行榜    

查看: 507  |  回复: 0

我是飞锅

新虫 (初入文坛)

[求助] 高分子链的模拟——积分的时间步长

我现在在做一些简单的高分子链的模拟,用珠链模型,单元与单元之间用LJ势,键用FENE势,用langevin热浴积分。问题是,我用来积分的时间步长只能取10-4的量级,否则的话键会被拉断得出,不物理的结果。很多文献上用类似的模拟方法,时间步长可以取到10-3甚至10-2的量级。模拟中时间步长的大小很关键,因为这直接影响了计算量。我仔细检查了自己的程序,目前为止没有找到错误。大家有做这方面的模拟的吗?这问题大概是出在什么地方了,望指点迷津,不甚感激!
回复此楼
致力于搞科研,就是这么简单
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 我是飞锅 的主题更新
信息提示
请填处理意见