干货分子动力学 干货详情

lammps中文教程3

春风得意马蹄疾30
分子动力学LAMMPS分子模拟计算材料学NVT系综热力学量初始化温度控制系综调控运行步数控制升温模拟

1.3 run命令使⽤技巧,多次run正确写法

本⽂介绍run命令的⼏种不同⽤法。

run命令是lammps模拟中必不可缺的⼀个命令,主要功能是设定运算步数并启动运算。

最常⻅的命令是run+步数,其实,run命令还有更⾼级的应⽤,下⾯分类介绍run命令的⽤法。

(1)运⾏N步:run N

run+步数是⽤的最多的⼀条命令,主要功能是运⾏N步。

例如,run 1000 指在当前状态下继续运⾏1000步。

还有⼀个⽐较特殊的⽤法是:run 0

run 0功能是在不启动积分运算的情况下,计算热⼒学信息,此功能可⽤来进⾏热⼒学变量的初始化,如压⼒、温度等。

(2)运⾏到N步:run N upto

加关键词upto,表示从当前步数开始,运⾏到第N步结束。

例如,当前步数为100步,run 300 upto 的意思是从第100步开始,运⾏到第300步,这个run语句实际运⾏步数为200步。

(3)run N start N1 stop N2

这个命令表示从N1到N2步数内,运算N步。例如

fix          1 all nvt 200.0 300.0 1.0

run        1000 start 0 stop 10000

run        1000 start 0 stop 10000

...

run        1000 start 0 stop 10000

上⾯这段代码每⼀个run命令运⾏1000步,10个run共运⾏10000步,在这10000步运⾏过程中,温度从200K到300K线性变化。

与之类似的命令:

fix 1 all nvt 200.0 300.0 1.0

run 1000

run 1000

...

run 1000


这种写法,每运⾏1000步,都要重新开始,按照系综的设置从200->300K进⾏升温,运⾏的结果类似下图:

(4) run N every N1 command

共运⾏N步,每隔N1步执⾏⼀⼀次command命令,例:

run 6000 every 2000 "print 'Coord = $q'"

共运⾏6000步,每2000步在屏幕上输出⼀⾏提示。

这个命令和下⾯的写法效果是⼀样的:

variable q equal x[100]

run 2000

print "Coord = $q"

run 2000

print "Coord = $q"

run 2000

print "Coord = $q"

评论

0/1000发布评论
全部评论

春风得意马蹄疾

关注
TA的主页

干货推荐