24小时热门版块排行榜    

Znn3bq.jpeg
汕头大学海洋科学、生物学、生物与医药接受调剂
查看: 1492  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0856专硕求调剂 希望是a区院校 +23 好好休息好不好 2026-04-09 26/1300 2026-04-13 14:22 by 张zhihao
[考研] 人工智能320调剂08工类还有机会吗 +15 振—TZ 2026-04-10 16/800 2026-04-13 14:11 by 张zhihao
[考研] 一志愿厦大生物学332求调剂 +11 池池池池池池 2026-04-08 11/550 2026-04-13 14:10 by 科研论
[考研] 生物学调剂,一志愿西南大学348,Top期刊一区二作、二区三作,三等奖学金三次 +5 candyyyi 2026-04-09 5/250 2026-04-13 09:02 by 可淡不可忘
[考研] 268分085602化学工程调剂 +31 月照花林。 2026-04-09 31/1550 2026-04-12 23:37 by 学zh
[考研] 326求调剂 +6 Shansyn 2026-04-10 6/300 2026-04-12 09:46 by hammer3
[考研] 327求调剂 +5 拾光任染 2026-04-06 5/250 2026-04-12 02:29 by 秋豆菜芽
[考研] 269求调剂 +11 啊啊我我 2026-04-07 11/550 2026-04-11 16:45 by vgtyfty
[考研] 调剂 +4 电气300求调剂不 2026-04-08 7/350 2026-04-11 10:44 by 紫曦紫棋
[考研] 22408调剂求助 +7 毂12 2026-04-09 9/450 2026-04-11 09:23 by 哦哦123
[考研] 284求调剂 +9 让我上岸吧阿西 2026-04-09 11/550 2026-04-10 19:18 by 靖jing
[考研] 0858求调剂 5+5 Gky09300550, 2026-04-10 8/400 2026-04-10 19:13 by chemisry
[考研] 273求调剂 +51 麦小叮当 2026-04-06 58/2900 2026-04-10 15:54 by jiajinhpu
[考研] 材料化工总分334求调剂 +16 Riot2025 2026-04-08 17/850 2026-04-09 20:19 by maddjdld
[考研] 328求调剂 +17 lftmya 2026-04-07 18/900 2026-04-09 08:05 by 5268321
[考研] 270求调剂 +3 031127 2026-04-06 4/200 2026-04-08 21:00 by 逆水乘风
[考研] 277求调剂 +4 考研调剂lxh 2026-04-06 6/300 2026-04-08 10:40 by 逆水乘风
[考研] 323求调剂 +3 林zlu 2026-04-07 4/200 2026-04-07 23:21 by lbsjt
[考研] 307求调剂 +3 Youth@@ 2026-04-07 3/150 2026-04-07 22:00 by hemengdong
[考研] 软工学硕299求调剂 +6 useryy 2026-04-07 6/300 2026-04-07 09:50 by vgtyfty
信息提示
请填处理意见