24小时热门版块排行榜    

查看: 1596  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 3xN9SualEUPt 2026-10-24 4/200 2026-10-25 14:42 by XiYP04It3Evb
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 3xN9SualEUPt 2026-10-24 3/150 2026-10-25 14:23 by XiYP04It3Evb
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 3xN9SualEUPt 2026-10-24 4/200 2026-10-25 14:07 by XiYP04It3Evb
[论文投稿] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 3xN9SualEUPt 2026-10-24 4/200 2026-10-25 13:54 by XiYP04It3Evb
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 DI1qFHFK0KbE 2026-10-24 4/200 2026-10-25 13:18 by XiYP04It3Evb
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 7xHjyB3bAiwg 2026-10-24 4/200 2026-10-25 13:15 by XiYP04It3Evb
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 7xHjyB3bAiwg 2026-10-23 8/400 2026-10-25 12:59 by XiYP04It3Evb
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +5 7xHjyB3bAiwg 2026-10-23 10/500 2026-10-25 12:51 by XiYP04It3Evb
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 7xHjyB3bAiwg 2026-10-23 10/500 2026-10-25 12:49 by XiYP04It3Evb
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 7xHjyB3bAiwg 2026-10-23 11/550 2026-10-25 12:45 by XiYP04It3Evb
[找工作] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 7xHjyB3bAiwg 2026-10-23 4/200 2026-10-25 12:36 by XiYP04It3Evb
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 7xHjyB3bAiwg 2026-10-23 4/200 2026-10-25 12:32 by XiYP04It3Evb
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 GVGRDgNfv7P4 2026-10-23 4/200 2026-10-25 12:10 by XiYP04It3Evb
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 GVGRDgNfv7P4 2026-10-23 9/450 2026-10-25 12:10 by XiYP04It3Evb
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 7xHjyB3bAiwg 2026-10-24 4/200 2026-10-25 01:22 by XcOYG53d3x4p
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 7xHjyB3bAiwg 2026-10-23 8/400 2026-10-25 00:15 by XcOYG53d3x4p
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 7xHjyB3bAiwg 2026-10-23 8/400 2026-10-24 22:56 by XcOYG53d3x4p
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 7xHjyB3bAiwg 2026-10-24 3/150 2026-10-24 18:23 by WqtG0Aem8Ufl
[文学芳草园] 亲密关系 +7 宋昙 2026-10-21 13/650 2026-10-23 17:07 by 小牧11
[考博] 郑州大学申博 +5 stars~ 2026-10-19 8/400 2026-10-22 15:06 by stars~
信息提示
请填处理意见