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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » LAMMPS/DL_POLY » 关于用Lammps模拟Cu熔点的例子问题
81/1返回列表
查看: 5196  |  回复: 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返回列表
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0703化学调剂 290分有科研经历，论文在投 +3 腻腻gk 2026-03-14 3/150 2026-03-15 17:28 by 小物理化学
[考研] 材料专硕326求调剂 +4 墨煜姒莘 2026-03-15 4/200 2026-03-15 11:02 by dyw
[考研] 085601材料工程315分求调剂 +3 yang_0104 2026-03-15 3/150 2026-03-15 10:58 by peike
[考研] 268求调剂 +5 一定有学上- 2026-03-14 6/300 2026-03-14 22:20 by 运气yunqi
[考研] 085600材料与化工 求调剂 +11 enenenhui 2026-03-13 12/600 2026-03-14 21:48 by laoshidan
[考研] 本科南京大学一志愿川大药学327 +3 麦田耕者 2026-03-14 3/150 2026-03-14 20:04 by 外星文明
[考研] 【0703化学调剂】-一志愿华中师范大学-六级475 +5 Becho359 2026-03-11 5/250 2026-03-14 11:35 by 哦哦123
[考研] 求调剂 +5 鹤遨予卿 2026-03-09 5/250 2026-03-14 02:44 by JourneyLucky
[考研] 288求调剂 +14 王晓阳- 2026-03-09 19/950 2026-03-14 02:05 by JourneyLucky
[考研] 311求调剂 +5 牛乳糖的卡卡 2026-03-10 5/250 2026-03-14 00:05 by JourneyLucky
[考研] 材料与化工304求B区调剂 +5 邱gl 2026-03-11 6/300 2026-03-13 22:37 by JourneyLucky
[考研] 308求调剂 +5 是Lupa啊 2026-03-11 5/250 2026-03-13 22:13 by JourneyLucky
[考研] 0856材料与化工301求调剂 +5 奕束光 2026-03-13 5/250 2026-03-13 22:00 by 星空星月
[考研] 333求调剂 +3 球球古力 2026-03-11 3/150 2026-03-13 21:27 by JourneyLucky
[考研] 311求调剂 +3 冬十三 2026-03-13 3/150 2026-03-13 20:41 by JourneyLucky
[考研] 295求调剂 +3 小匕仔汁 2026-03-12 3/150 2026-03-13 15:17 by vgtyfty
[论文投稿] 投稿问题 5+4 星光灿烂xt 2026-03-12 6/300 2026-03-13 14:17 by god_tian
[考研] 321求调剂（食品/专硕） +3 xc321 2026-03-12 6/300 2026-03-13 08:45 by xc321
[考研] 081200-11408-276学硕求调剂 +3 崔wj 2026-03-12 4/200 2026-03-12 19:33 by 求调剂zz
[考研] 0857 资源与环境 285分 +6 未名考生 2026-03-09 6/300 2026-03-11 21:08 by 30660438
信息提示
关闭
请填处理意见
关闭