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

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » LAMMPS/DL_POLY » 关于用Lammps模拟Cu熔点的例子问题
81/1返回列表
查看: 5044  |  回复: 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返回列表
信息提示
关闭
请填处理意见
关闭