24小时热门版块排行榜    

查看: 1738  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 放逐孤岛 的主题更新
信息提示
请填处理意见