| 查看: 1812 | 回复: 1 | |||
liuchenhan新虫 (小有名气)
|
[求助]
关于石墨加压的lammps模拟
|
|
小弟希望做lammps加压下石墨晶格常数的变化。 要求：c-axis 方向加1GPa，面内不加压，自由变化，温度要设在300K.以下是我的脚本。 # simulation # Initialization variable dt equal 0.0005 variable Tstart equal 1000.0 variable Tstop equal 300.0 variable Tdamp equal 100.0*${dt} variable Pstart equal 10.0 variable Pstop equal 10.0 variable Pdamp equal 1000.0*${dt} units metal dimension 3 newton on boundary p p p atom_style atomic neighbor 0.3 nsq neigh_modify check yes read_data nanonew.xyz # tersoff pair_style hybrid tersoff lj/cut 8.5 pair_coeff * * tersoff Cnew.tersoff C C C C C C C C pair_coeff 1 2 lj/cut 0.0024 3.4 #8.5 3.4 pair_coeff 1 3 lj/cut 0.0024 3.4 pair_coeff 1 4 lj/cut 0.0024 3.4 pair_coeff 1 5 lj/cut 0.0024 3.4 pair_coeff 1 6 lj/cut 0.0024 3.4 pair_coeff 1 7 lj/cut 0.0024 3.4 pair_coeff 1 8 lj/cut 0.0024 3.4 pair_coeff 2 3 lj/cut 0.0024 3.4 pair_coeff 2 4 lj/cut 0.0024 3.4 pair_coeff 2 5 lj/cut 0.0024 3.4 pair_coeff 2 6 lj/cut 0.0024 3.4 pair_coeff 2 7 lj/cut 0.0024 3.4 pair_coeff 2 8 lj/cut 0.0024 3.4 pair_coeff 3 4 lj/cut 0.0024 3.4 pair_coeff 3 5 lj/cut 0.0024 3.4 pair_coeff 3 6 lj/cut 0.0024 3.4 pair_coeff 3 7 lj/cut 0.0024 3.4 pair_coeff 3 8 lj/cut 0.0024 3.4 pair_coeff 4 5 lj/cut 0.0024 3.4 pair_coeff 4 6 lj/cut 0.0024 3.4 pair_coeff 4 7 lj/cut 0.0024 3.4 pair_coeff 4 8 lj/cut 0.0024 3.4 pair_coeff 5 6 lj/cut 0.0024 3.4 pair_coeff 5 7 lj/cut 0.0024 3.4 pair_coeff 5 8 lj/cut 0.0024 3.4 pair_coeff 6 7 lj/cut 0.0024 3.4 pair_coeff 6 8 lj/cut 0.0024 3.4 pair_coeff 7 8 lj/cut 0.0024 3.4 velocity all create ${Tstart} 1000001 mom yes rot yes dist gaussian fix NPT all npt temp ${Tstart} ${Tstop} ${Tdamp} x 10.0 10.0 1000.0*${dt} y 10.0 10.0 1000.0*${dt} z 1000.0 1000.0 1000.0*${dt} variable V equal vol variable LX equal lx variable LY equal ly variable LZ equal lz variable T equal temp variable PE equal pe variable E equal etotal fix VOL all print 100 "$V ${LX} ${LY} ${LZ}" file volume.txt screen no fix TE all print 100 "$T ${PE} $E" file TE.txt screen no run 1000000 问题：1、周期性边界条件加的对不对，不对的话应该怎么改; 2、压强加的对不对，不对的话应该怎么改；加压强时需要周期性边界条件吗？ 3、脚本还有哪些错误的地方，麻烦指出 非常感谢大神能为小弟指点迷津，刚刚入门，还望海涵，拜托啦。 |
回复此楼» 猜你喜欢
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有3人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有3人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有4人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- MS修改键长，MS画石墨烯，MS文件导入LAMMPS
已经有3人回复
- 石墨烯 lammps
已经有11人回复
liuchenhan
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1536.5
- 散金: 4
- 帖子: 170
- 在线: 108小时
- 虫号: 1508714
- 注册: 2011-11-25
- 专业: 传热传质学
2楼2014-02-25 17:38:48
20