应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » LAMMPS/DL_POLY » 关于用Lammps模拟Cu熔点的例子问题
北京石油化工学院2026年研究生招生接收调剂公告
81/1返回列表
查看: 5246  |  回复: 7
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

dd5818

银虫 (正式写手)

[求助] 关于用Lammps模拟Cu熔点的例子问题

输入文件:

# bulk Cu lattice

variable        x index 1
variable        y index 1
variable        z index 1

variable        xx equal 20*$x
variable        yy equal 20*$y
variable        zz equal 20*$z

units           metal
atom_style      atomic

lattice         fcc 3.615
region          box block 0 ${xx} 0 ${yy} 0 ${zz}
create_box      1 box
create_atoms    1 box

pair_style      eam
pair_coeff      1 1 Cu_u3.eam

velocity        all create 1600.0 376847 loop geom

neighbor        1.0 bin
neigh_modify    every 1 delay 5 check yes

fix             1 all nve

timestep        0.005
thermo          50

dump            1 all atom 10 dump.eam

run             1000
输出结果:

LAMMPS (16 Dec 2011)
Lattice spacing in x,y,z = 3.615 3.615 3.615
Created orthogonal box = (0 0 0) to (72.3 72.3 72.3)
  1 by 1 by 2 MPI processor grid
Created 32000 atoms
Setting up run ...
Memory usage per processor = 7.75649 Mbytes
Step Temp E_pair E_mol TotEng Press
       0         1600      -113280            0   -106662.09    18703.573
      50    781.69049   -109873.35            0   -106640.13    52273.088
     100      801.832    -109957.3            0   -106640.77    51322.821
     150    794.67685   -109927.54            0    -106640.6     51669.28
     200    795.75217   -109932.05            0   -106640.66    51628.945
     250    797.58407   -109939.57            0   -106640.61    51564.207
     300    793.23426    -109921.5            0   -106640.53    51741.321
     350    799.19313    -109946.3            0   -106640.68    51461.975
     400    797.16375   -109937.82            0    -106640.6    51585.836
     450    794.86113    -109928.1            0    -106640.4    51651.514
     500    796.77976   -109936.18            0   -106640.54    51565.369
     550    796.68581   -109935.87            0   -106640.62    51581.061
     600    799.54571   -109947.76            0   -106640.68    51434.487
     650    800.84902   -109953.13            0   -106640.66    51407.739
     700    791.74048   -109915.12            0   -106640.33    51729.628
     750    801.43616   -109955.62            0   -106640.72    51340.419
     800    797.91463   -109940.73            0    -106640.4    51527.856
     850    803.71916   -109964.98            0   -106640.64    51321.697
     900    805.40487   -109972.01            0    -106640.7    51275.268
     950    799.22869   -109946.21            0   -106640.44    51489.219
    1000    796.34045   -109934.01            0   -106640.19    51606.217
Loop time of 103.349 on 2 procs for 1000 steps with 32000 atoms

Pair  time (%) = 75.6815 (73.2291)
Neigh time (%) = 7.94726 (7.68973)
Comm  time (%) = 8.33589 (8.06577)
Outpt time (%) = 7.90816 (7.6519)
Other time (%) = 3.4762 (3.36355)

Nlocal:    16000 ave 16044 max 15956 min
Histogram: 1 0 0 0 0 0 0 0 0 1
Nghost:    13820.5 ave 13865 max 13776 min
Histogram: 1 0 0 0 0 0 0 0 0 1
Neighs:    604056 ave 605895 max 602218 min
Histogram: 1 0 0 0 0 0 0 0 0 1

Total # of neighbors = 1208113
Ave neighs/atom = 37.7535
Neighbor list builds = 140
Dangerous builds = 1
问题：
1) 1000步后体系是否达到平衡?
2) 能否说模拟结果给出的Cu的熔点为800K左右?
3) 这个例子是否可以模拟Cu的熔点?
回复此楼

» 收录本帖的淘帖专辑推荐

动力学强帖

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2012-03-08 15:56:29
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

w76990480

荣誉版主 (文坛精英)

super-emuch

优秀版主文献杰出贡献文献杰出贡献

【答案】应助回帖

★ ★
感谢参与，应助指数 +1
dd5818(金币+2): 2012-03-09 08:36:46
西安交通大学有个计算熔点的例子，不知道你看过没有，平不平衡你要画出图给大家看看，估计画出图你自己就能看出来
赞一下 回复此楼
专业文献求助
2楼2012-03-08 23:15:54
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

dd5818

银虫 (正式写手)
引用回帖:
楼: Originally posted by w76990480 at 2012-03-08 23:15:54:
西安交通大学有个计算熔点的例子，不知道你看过没有，平不平衡你要画出图给大家看看，估计画出图你自己就能看出来

但是用什么数据绘图呢
赞一下 回复此楼
3楼2012-03-09 08:37:04
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

w76990480

荣誉版主 (文坛精英)

super-emuch

优秀版主文献杰出贡献文献杰出贡献

【答案】应助回帖

★ ★
dd5818(金币+2): 2012-03-09 13:26:26
内容已删除
赞一下(1人) 回复此楼
专业文献求助
4楼2012-03-09 08:42:33
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

dd5818

银虫 (正式写手)
内容已删除
回复此楼
5楼2012-03-09 13:26:14
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

dd5818

银虫 (正式写手)
我按照例子运行了, 但是上面的那个Cu的势能我没有, 我用的Cu_u3.eam势能,
绘制出了T-V散点图, 但是没有例子中出现的那个突变点，这又是为什么呢？？
请指教
赞一下 回复此楼
6楼2012-03-09 21:34:29
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

dd5818

银虫 (正式写手)

如图上图所示
赞一下 回复此楼
7楼2012-03-09 21:43:44
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

abinitio

版主 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
chaizhm: 金币+3, 谢谢~ 2012-04-27 19:44:08
dd5818: 金币+3 2012-04-27 20:17:47
正好这两天也在用lammps练习Cu熔化的例子，对楼主的问题回答如下：
1) 1000步后体系是否达到平衡?
   A:从最后的温度值看，远未达到平衡。1000步太少了，建议加到10w或100w再试。
2) 能否说模拟结果给出的Cu的熔点为800K左右?
   A: 找熔点的可以用npt升温的方法，找V-t或pe-t的突变点，如果用nvt系综，一般是在理论熔化附近的温度内每个温度都进行充分的nvt驰誉，对比不同温度下的rdf或pe值。
3) 这个例子是否可以模拟Cu的熔点?
   A：个人认为，在nvt下只使用一下温度，不能模拟Cu的熔点。
赞一下(1人) 回复此楼
8楼2012-04-27 16:03:44
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 dd5818 的主题更新
81/1返回列表
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 320分人工智能调剂 +8 振—TZ 2026-04-03 8/400 2026-04-05 22:33 by 范式思维
[考研] 277求调剂 数一104分 +6 瓶子PZ 2026-04-05 6/300 2026-04-05 20:38 by 啵啵啵0119
[考研] 288求调剂 +7 没有答案_ 2026-04-05 7/350 2026-04-05 17:01 by yulian1987
[考研] 0703化学调剂325分 +9 15771691647 2026-04-04 9/450 2026-04-05 11:39 by 猪会飞
[考研] 材料调剂 +7 dxy调剂 2026-04-04 7/350 2026-04-05 09:15 by 陌秋26
[考研] 一志愿211，化学学硕，310分，本科重点双非，求调剂 +11 努力奋斗112 2026-04-04 11/550 2026-04-04 20:51 by 蓝云思雨
[考研] 324求调剂 +14 想上学求调 2026-04-02 15/750 2026-04-04 20:31 by 无际的草原
[考研] 265求调剂 +20 梁梁校校 2026-04-01 21/1050 2026-04-04 00:38 by userper
[考研] 320调剂 +4 农业工程与信息� 2026-04-03 4/200 2026-04-03 21:40 by lbsjt
[考研] 070300一志愿211，312分求调剂院校 +16 小黄鸭宝 2026-03-30 16/800 2026-04-03 19:53 by lijunpoly
[考研] 学硕机械工程303求调剂 +6 无名所以叫吴明 2026-03-30 7/350 2026-04-03 16:48 by asdfzly
[考研] 319求调剂 +18 太容易1018 2026-04-01 18/900 2026-04-03 11:18 by linyelide
[考研] 285求调剂 +14 AZMK 2026-04-02 14/700 2026-04-02 15:54 by 上九天揽月（好�
[考研] 271求调剂 +15 勒布朗@ 2026-03-31 20/1000 2026-04-02 11:24 by Sammy2
[考研] 考研调剂 +12 Amber00 2026-03-31 12/600 2026-04-02 09:04 by sanrepian
[考研] 0703一志愿南师大334求调剂 +4 seven7yu 2026-03-30 4/200 2026-04-01 16:10 by oooqiao
[硕博家园] 博一被送出联培感觉不适应怎么办 +3 全村的狗 2026-03-31 3/150 2026-04-01 10:44 by 328838485
[考研] 267求调剂 +13 uiybh 2026-03-31 13/650 2026-04-01 10:25 by 探123
[考研] 一志愿大连理工大学，机械工程学硕，341 +3 西瓜田的守望者 2026-03-30 3/150 2026-03-31 11:08 by asdfzly
[考研] 083000环境科学与工程调剂，总分281 +4 橙子（胜意） 2026-03-30 4/200 2026-03-31 00:44 by Linzejun
信息提示
关闭
请填处理意见
关闭