24小时热门版块排行榜    

查看: 9025  |  回复: 16
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

xiaorongwei

金虫 (初入文坛)

[交流] lammps的时间步长问题 已有10人参与

我的in文件中units 是metal 默认的时间步长是0.001ps,而我现在模拟的时间比较长,所以我将时间步长timestep提高了10倍,即0.01。建立的模型也可以跑,暂时没有出现什么异常。  不知道各位对时间步长这么大有没有什么看法,会不会对计算结果有影响?
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chaizhm

木虫 (著名写手)

★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖
御剑江湖(金币+3): 谢谢 2011-12-10 20:55:38
引用回帖:
4楼: Originally posted by xiaorongwei at 2011-12-09 19:30:18:
那在提高速度方面还有其他比较好的方法吗?

你看看在体系大小方面能不能做点文章,但是体系也不能太小,因为分子模拟好多参数都是通过统计平均得到,系统太小结果也不准确
壁立万仞,无欲则刚
5楼2011-12-10 13:04:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 17 个回答

rencl1018

铜虫 (初入文坛)

★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖
chaizhm(金币+2): 谢谢交流 2011-12-09 10:31:21
chaizhm(金币+2): 谢谢交流 2011-12-09 10:31:24
Lammps里面使用的是Velocity Verlet 积分时间步长算法,较大的时间步长会使能量的涨落比较大,10fs有点大了。
2楼2011-12-09 10:00:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chaizhm

木虫 (著名写手)

★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖
御剑江湖(金币+3): 谢谢 2011-12-10 20:55:41
jiaoyixiong(金币+2): 鼓励交流 2011-12-11 22:49:46
引用回帖:
1楼: Originally posted by xiaorongwei at 2011-12-09 08:17:45:
我的in文件中units 是metal 默认的时间步长是0.001ps,而我现在模拟的时间比较长,所以我将时间步长timestep提高了10倍,即0.01。建立的模型也可以跑,暂时没有出现什么异常。  不知道各位对时间步长这么大有没有 ...

不建议你通过增加时间步长的方法来提高速度。因为分子模拟用到很多数值方法对微分方程进行离散化,这样如果取的间隔太大,微积分就不满足了,系统很容易发散,对结果也会有比较大的影响
壁立万仞,无欲则刚
3楼2011-12-09 10:33:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaorongwei

金虫 (初入文坛)

引用回帖:
3楼: Originally posted by chaizhm at 2011-12-09 10:33:05:
不建议你通过增加时间步长的方法来提高速度。因为分子模拟用到很多数值方法对微分方程进行离散化,这样如果取的间隔太大,微积分就不满足了,系统很容易发散,对结果也会有比较大的影响

那在提高速度方面还有其他比较好的方法吗?
4楼2011-12-09 19:30:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 08:14 by C79jjtjAKjEn
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-15 08:02 by C79jjtjAKjEn
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-15 08:02 by C79jjtjAKjEn
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 s3fFTmArrBt6 2026-09-14 4/200 2026-09-15 07:50 by C79jjtjAKjEn
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +5 s3fFTmArrBt6 2026-09-14 6/300 2026-09-15 07:38 by C79jjtjAKjEn
[考博] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +6 s3fFTmArrBt6 2026-09-14 6/300 2026-09-15 07:14 by C79jjtjAKjEn
[教师之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 s3fFTmArrBt6 2026-09-13 4/200 2026-09-15 07:02 by 5BDX0d0WFp7t
[公派出国] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 s3fFTmArrBt6 2026-09-13 5/250 2026-09-15 06:38 by 5BDX0d0WFp7t
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 LwdutQ8HoqWP 2026-09-13 4/200 2026-09-15 06:14 by 5BDX0d0WFp7t
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 LwdutQ8HoqWP 2026-09-13 4/200 2026-09-15 06:02 by 5BDX0d0WFp7t
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 3n8v2C8RimXI 2026-09-13 4/200 2026-09-15 05:50 by 5BDX0d0WFp7t
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +6 3n8v2C8RimXI 2026-09-13 6/300 2026-09-15 05:50 by 5BDX0d0WFp7t
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 3n8v2C8RimXI 2026-09-13 3/150 2026-09-15 05:49 by 5BDX0d0WFp7t
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +6 QUjhNVAcOSff 2026-09-13 7/350 2026-09-15 05:22 by 5BDX0d0WFp7t
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +5 QUjhNVAcOSff 2026-09-13 6/300 2026-09-15 05:13 by 5BDX0d0WFp7t
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +5 QUjhNVAcOSff 2026-09-13 6/300 2026-09-15 05:01 by mibUvS8DDCwf
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 23jxep3nCNZb 2026-09-14 3/150 2026-09-15 01:21 by mibUvS8DDCwf
[考博] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +5 QUjhNVAcOSff 2026-09-13 5/250 2026-09-14 21:36 by vZfe6xYu34yj
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 LwdutQ8HoqWP 2026-09-13 4/200 2026-09-14 16:13 by Aj1rhIDL5ixY
[考博] 27秋季材料博士申请求助 10+3 Cahal1 2026-09-11 4/200 2026-09-14 08:50 by 北京莱茵编辑
信息提示
请填处理意见