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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » LAMMPS/DL_POLY » 关于用Lammps模拟Cu熔点的例子问题
51/1返回列表
查看: 5564  |  回复: 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的回帖

dd5818

银虫 (正式写手)

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

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的回帖

dd5818

银虫 (正式写手)
我按照例子运行了, 但是上面的那个Cu的势能我没有, 我用的Cu_u3.eam势能,
绘制出了T-V散点图, 但是没有例子中出现的那个突变点，这又是为什么呢？？
请指教
赞一下 回复此楼
6楼2012-03-09 21:34:29
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +8 医学老男孩 2026-08-13 23/1150 2026-08-13 20:42 by 医学老男孩
[硕博家园] 一作与独作在应聘高校教师时区别大吗 +3 mbygzh 2026-08-08 4/200 2026-08-13 19:31 by 龙-樱
[基金申请] FileCode能看出啥？ +10 要乐观耀哥 2026-08-10 28/1400 2026-08-13 17:51 by 医学老男孩
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +22 Lanmanbaby 2026-08-09 36/1800 2026-08-12 22:50 by sdfapple719
[基金申请] 好奇怪的filecode +5 布布和一二 2026-08-08 6/300 2026-08-12 16:11 by 云上清扬
[基金申请] 奇怪，两个人的filecode固定段从头到尾一模一样 +7 布布和一二 2026-08-10 10/500 2026-08-12 09:12 by 布布和一二
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分？ +11 yufeiwaner 2026-08-09 13/650 2026-08-12 08:17 by yufeiwaner
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[硕博家园] 读博的好处 +3 lnee 2026-08-11 3/150 2026-08-11 18:10 by 希望我好好的
[基金申请] 什么时候出结果，有咨询渠道？？？ +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] 确定了，国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
[基金申请] 2026国自然放榜时间 +9 布布和一二 2026-08-08 9/450 2026-08-10 11:22 by xxxx2020
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
[基金申请] filecode与中标关系的预测 +5 布布和一二 2026-08-07 5/250 2026-08-09 16:15 by 袁向阳007
[基金申请] fileCode有新解读？ +10 Tide man 2026-08-08 18/900 2026-08-09 12:55 by 仁砚薪传
[基金申请] 关于filecode，很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
[基金申请] 国基金的申报应该改成非等额制，评价高的钱多评价低的钱少，但是增加资助率 +7 a089 2026-08-07 7/350 2026-08-08 18:05 by gltch
信息提示
关闭
请填处理意见
关闭