24小时热门版块排行榜    

查看: 1953  |  回复: 2

放逐孤岛

铜虫 (小有名气)

[求助] lammps中fix_npt命令详解

我在使用fix_npt命令时,产生了以下疑问:

fix                     1 all npt temp 298 298 $(100.0*dt) iso 1.0 -200 1000.0
timestep           0.5
run                   100

原型即:
fix                     ID group-ID npt temp Tstart Tstop Tdamp iso Pstart Pstop Pdamp

这个是我的代码段,我对Tdamp与Pdamp产生了疑问。比如压力控制时,我在100个时间步内让压力从1atm降到-200atm,后面的1000起什么样的作用呢?或者说,有没有一个公式包含Tstart,Tstop,Tdamp能够描述压力降低的过程呢?感谢各位虫友@smutao@smutao
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Nater1ver

木虫 (正式写手)

是时间间隔,它的倒数是频率,影响的是温度/压强变化的快慢,一般自己估摸着取就行了,根据你的单位决定,如果变化太慢了就要调小,每个人习惯不一样,1000可能太大了,有时候我在real单位会考虑25左右,视不同的体系和单位,这个值是可以设置到很小的,比如1.0甚至更小,要看你的需求和调试过程,而且你100步跑不到-200atm的话,步数还要加,但是要注意一点的是,MD中压强的涨落是很大的,那么你很有可能体系的压强一下变成-4000或者3000,这个1~-200如果只看压强可能不是那么明显,那就需要取多个点的平衡涨落做系综平均了

发自小木虫Android客户端

» 本帖已获得的红花(最新10朵)

2楼2019-03-27 23:57:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

放逐孤岛

铜虫 (小有名气)

送红花一朵
引用回帖:
2楼: Originally posted by Nater1ver at 2019-03-27 23:57:10
是时间间隔,它的倒数是频率,影响的是温度/压强变化的快慢,一般自己估摸着取就行了,根据你的单位决定,如果变化太慢了就要调小,每个人习惯不一样,1000可能太大了,有时候我在real单位会考虑25左右,视不同的体 ...

哇,感谢,您回答的太详了,谢谢

发自小木虫Android客户端
3楼2019-03-28 08:56:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 放逐孤岛 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 明天E口面上会评 +4 布布和一二 2026-07-12 5/250 2026-07-13 10:44 by zzybjr
[公派出国] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 obzynn3tuv 2026-07-11 4/200 2026-07-13 10:35 by emlgkvoijw
[考博] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 obzynn3tuv 2026-07-11 4/200 2026-07-13 10:25 by emlgkvoijw
[论文投稿] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 obzynn3tuv 2026-07-11 4/200 2026-07-13 10:15 by emlgkvoijw
[硕博家园] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 obzynn3tuv 2026-07-11 3/150 2026-07-13 09:56 by emlgkvoijw
[考博] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 obzynn3tuv 2026-07-11 4/200 2026-07-13 00:24 by v2i90kmrcu
[考研] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 obzynn3tuv 2026-07-11 4/200 2026-07-13 00:15 by v2i90kmrcu
[基金申请] 不要再数国自然申请书的 filecode 的分隔符个数了 +10 sadpencil 2026-07-12 15/750 2026-07-12 21:50 by sadpencil
[基金申请] 祈祷青基必中 50+6 hadengry 2026-07-09 16/800 2026-07-12 20:31 by 登山虫虫
[基金申请] 生命口会评 +11 wuke100666 2026-07-07 12/600 2026-07-11 19:30 by hhs666
[考博] 27届辽宁大学应届毕业生申博 +3 可乐微风大海 2026-07-09 3/150 2026-07-10 21:23 by 贪玩的研究僧
[基金申请] 关于如何从代码看上不上会 +17 且听虎啸 2026-07-09 23/1150 2026-07-10 19:51 by sdfapple719
[有机交流] chemdraw 20+5 ?慕鱼? 2026-07-08 6/300 2026-07-10 15:37 by 紫枫星域
[基金申请] b口会评 +6 当空照 2026-07-08 6/300 2026-07-09 20:49 by tonyliu0401
[基金申请] 信息学部会评时间 +3 xuel2011 2026-07-06 3/150 2026-07-09 19:05 by mollyzhang_2003
[考研] 在职考研 +6 wu5d 2026-07-08 7/350 2026-07-09 16:40 by 化工小霸王呀
[论文投稿] MSER送审了还被拒稿 +3 S778950906 2026-07-08 5/250 2026-07-09 15:46 by S778950906
[有机交流] 硝基还原 20+3 暮年飞雪 2026-07-07 5/250 2026-07-09 13:30 by tfang
[基金申请] 化学口会评?应用基础研究和基础研究,希望能分开 +3 Tide man 2026-07-08 3/150 2026-07-09 11:36 by dede5556
[有机交流] Suzuki偶联中硼酸酯原料掉下的硼酸酯是否可以自身偶联,生成B2Pin2杂质 500+3 天下大同 2026-07-06 3/150 2026-07-07 10:31 by martyr_oyx
信息提示
请填处理意见