| 查看: 359 | 回复: 2 | ||
[求助]
请大家帮我看看，为什么不能增加原子
|
|
# 3d EAM V surface diffusion units metal boundary p p s atom_style atomic #-------------------create simulation lattice-------------------# lattice bcc 3.04139 orient x 1 0 0 orient y 0 1 -1 orient z 0 1 1 region box block 0 8 0 8 0 16 #To define 5 lattice create_box 2 box create_atoms 1 box #----------------to define the potential--------------# pair_style meam pair_coeff * * library.meam V V.meam V V #---------------define the top layers and construct the vacuum-----------# region TopAtom block INF INF INF INF 16 INF #To define a delete region group Top region TopAtom delete_atoms group Top region BotAtom block INF INF INF INF INF 0 #To define a delete region group Bot region BotAtom delete_atoms group Bot #--------------constraint the adatom diffuse in the surface-------------# #fix constraint1 adAtom1 setforce NULL NULL 0.0 #velocity adAtom1 set NULL NULL 0.0 #fix constraint2 adAtom2 setforce NULL NULL 0.0 #velocity adAtom2 set NULL NULL 0.0 #group surf subtract all Bottom timestep 0.001 #dump id1 all xyz 1 initial_coor.xyz run 0 variable E equal pe variable N equal atoms print "the number of atoms & system energy now are $N $E" neighbor 2 bin neigh_modify every 20 delay 0 check yes #velocity surf create 300.0 4928459 rot yes dist gaussian print "nvt performed, temperature up:$N atoms,total energy is $E" compute 4 all temp thermo_style custom step temp etotal press thermo_modify temp 4 thermo 100 dump id3 all xyz 1000 snap_all300k.xyz #---------------set the NVT resemble and relax ---------------# min_style sd minimize 1.0e-15 1.0e-10 1000 10000 fix myRelax all nvt temp 10 100.0 0.1 run 10000 minimize 1.0e-20 1.0e-15 1000 10000 unfix myRelax fix myRelax1 all nvt temp 100 100 0.1 run 5000 unfix myRelax1 fix nextRelax all nvt temp 100.0 300.0 0.1 run 10000 minimize 1.0e-30 1.0e-20 1000 10000 unfix nextRelax fix nextRelax1 all nvt temp 300 300 0.1 run 20000 unfix nextRelax1 fix relax all nve run 10000 unfix relax #----------------add adatom----------------# create_atoms 2 single 4 4 16 units lattice region adAtom1 block 4 4 4 4 16 16 units lattice group adAtom1 region adAtom1 create_atoms 2 single 4 4 0 units lattice region adAtom2 block 4 4 4 4 0 0 units lattice group adAtom2 region adAtom2 .... log文件中提示增加0个原子, ??????? |
回复此楼» 猜你喜欢
297，工科调剂?河南农业大学本科
已经有12人回复
-
化学070300 求调剂
已经有24人回复
-
0854求调剂
已经有19人回复
-
294求调剂
已经有11人回复
-
一志愿华中农业071010，320求调剂
已经有14人回复
-
药学求调剂
已经有16人回复
-
291求调剂
已经有8人回复
-
新一代电子信息294求调剂 不挑学校
已经有9人回复
-
求调剂
已经有23人回复
-
290调剂生物0860
已经有44人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 请大家帮忙看看这个XRD图谱，谢谢
已经有11人回复
- 大家帮我看看我这样的青年基金梯队合理不？
已经有29人回复
- 请大家帮我看看这个计算结果，多谢！
已经有7人回复
- 投稿有答复了，些许激动请大家帮忙看看，谢谢大家！（在13楼增加第二封审稿意见）
已经有32人回复
- 大家来帮我看看这两个INCAR 文件有什么区别哦。。为啥计算结果不一样呢
已经有4人回复
- 第一次写文章，大家帮我看看这种情况能不能发表
已经有10人回复
- 请大侠帮看一下74个原子ENCUT测试应该选哪个值？
已经有9人回复
- 请大家帮我看看我的电猫怎么回事
已经有4人回复
- 【求助】大家帮我看看这多出来的10个原子是从哪里来的？
已经有4人回复
- 请大家帮我看看，这个博后读还是不读？
已经有34人回复
- 【求助/交流】请大家帮忙看看这个柱子为什么这样？
已经有13人回复
phq323
至尊木虫 (职业作家)
- 应助: 13 (小学生)
- 金币: 12622.1
- 红花: 1
- 帖子: 3487
- 在线: 348.1小时
- 虫号: 703929
- 注册: 2009-02-19
- 性别: GG
- 专业: 金属材料的微观结构
2楼2012-01-30 20:13:47
3楼2012-01-30 21:29:49
