24小时热门版块排行榜    

Znn3bq.jpeg
汕头大学海洋科学接受调剂
查看: 1493  |  回复: 1

anviacc

新虫 (初入文坛)

[求助] 求助 用lammps模拟铜纳米熔化过程, 不同温度下的rdf不变? 已有1人参与

初学lammps,求铜纳米熔化过程不同温度下的rdf函数? 不知哪里出了错误,求大神指点



clear
units metal
dimension 3
boundary p p p
atom_style atomic
atom_modify map array

# ---------- Create Atoms ---------------------

lattice fcc 3.61
region box block 0 200 0 200 0 200 units lattice
create_box 1 box


lattice fcc 3.61  
region my1 sphere 100.0 100.0 100.0 2.8
create_atoms 1 region my1


# ---------- Define Interatomic Potential ---------------------

pair_style eam/alloy
pair_coeff * * Cu_mishin1.eam.alloy Cu




# ---------- Run Minimization ---------------------
reset_timestep 0
velocity all create 700.0 4928459
#fix 1 all npt temp 900 900 1 iso 0 0 1000

neighbor 1.0 bin
neigh_modify delay 5 check yes

compute myRDF all rdf 50
fix 1 all ave/time 200 1 200 c_myRDF
file tmp.rdf mode vector

thermo 1000
thermo_style custom step temp pe lx ly lz
dump my all movie 10 my.avi type type
dump my1 all xyz 1000 file.xtc
run 50000

print "Total energy (eV) = ${teng};"
print "Number of atoms = ${natoms};"
print "Lattice constant (Angstoms) = ${length};"
print "Cohesive energy (eV) = ${ecoh};"

print "All done!"
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

brucefan

专家顾问 (著名写手)

【答案】应助回帖

感谢参与,应助指数 +1
你为什么把fix npt 注释掉?没它怎么积分?

[ 发自手机版 http://muchong.com/3g ]
2楼2018-03-24 20:37:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 anviacc 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 考研调剂 +7 长弓傲 2026-04-13 8/400 2026-04-13 20:46 by 长弓傲
[考研] 273求调剂 +3 白居不易. 2026-04-09 5/250 2026-04-13 19:27 by or7ro
[考研] 280求调剂 +13 wzzz王 2026-04-09 13/650 2026-04-12 00:31 by 勇攀高峰0126
[考研] 331求调剂 +5 王国帅 2026-04-11 5/250 2026-04-11 22:56 by 溪涧流水
[考研] 359求调剂 +5 胃痉挛累了 2026-04-11 5/250 2026-04-11 19:55 by lbsjt
[考研] 求调剂,一志愿大连理工大学354分 +5 雨声余生 2026-04-11 6/300 2026-04-11 16:12 by 雨声余生
[考研] 农业管理302分求调剂 +3 xuening1 2026-04-10 3/150 2026-04-11 10:18 by zhq0425
[考研] 298求调剂 +9 钉叮咚冬瓜 2026-04-07 11/550 2026-04-11 09:35 by zhq0425
[考研] 一志愿北理工298英一数二已上岸,感谢各位老师 +14 Reframe 2026-04-10 16/800 2026-04-10 23:07 by caotw2020
[考研] 263能源动力专硕求调剂 +3 加大号饭盒袋 2026-04-10 3/150 2026-04-10 22:23 by 286640313
[考研] 调剂 +19 小张ZA 2026-04-10 20/1000 2026-04-10 22:08 by 猪会飞
[考研] 085800 能源动力求调剂 +6 阿biu啊啊啊啊啊 2026-04-10 6/300 2026-04-10 15:03 by hemengdong
[考研] 314求调剂 +14 weltZeng 2026-04-09 14/700 2026-04-09 23:14 by wolf97
[考研] 材料化工总分334求调剂 +16 Riot2025 2026-04-08 17/850 2026-04-09 20:19 by maddjdld
[考研] 考研调剂 +13 冰冰,,, 2026-04-07 13/650 2026-04-09 17:01 by Lilly_Li
[考研] 求调剂材料科学与工程一志愿985初试365分 +5 材化李可 2026-04-08 5/250 2026-04-09 17:00 by Lilly_Li
[考研] 二次调剂求老师收留 +3 笑笑袁 2026-04-08 3/150 2026-04-08 23:50 by 醉在风里
[考研] 机械专硕273请求调剂 +6 庚申壬申 2026-04-07 6/300 2026-04-08 22:41 by bljnqdcc
[考研] 263分B区求调剂 +6 李nihao 2026-04-08 6/300 2026-04-08 09:38 by 南开小綦
[考研] 325 调剂 +6 QQ小虾 2026-04-07 6/300 2026-04-07 15:17 by Ccclqqq
信息提示
请填处理意见