| 查看: 5576 | 回复: 7 | |||
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的熔点? |
回复此楼» 收录本帖的淘帖专辑推荐
 动力学强帖 |
» 猜你喜欢
2026国自然函评费到账
已经有17人回复
-
让我中一个面上吧！
已经有15人回复
-
2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战
已经有5人回复
-
跳槽后在研项目怎么办？
已经有10人回复
-
什么时候开奖？
已经有10人回复
-
今天放榜吗？
已经有16人回复
-
93BebMhtakh前后11位开头都是大写
已经有8人回复
-
范进中举一文的中心思想
已经有4人回复
-
只有每年这种时候来逛逛小木虫
已经有26人回复
-
filecode，4个jtjc了
已经有16人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- lammps的时间步长问题
已经有16人回复
- 寻求lammps的应用 欢迎大家畅所欲言
已经有8人回复
- 【500金币本版活动】“经验力场分子模拟+第一性原理计算问题讨论互助及资源共享”
已经有82人回复
- 【求助】lammps安装过程中，在安装fftw的过程中遇到问题
已经有8人回复
- 【讨论】想学lammps是否应该先学习C++
已经有28人回复
- 【讨论】lammps里的boundary的位置能够随意更改？
已经有14人回复
- 【讨论】lammps 用sw势计算硅的熔点怎么这么高（2400k),附in文件
已经有13人回复
- 【讨论】lammps中利用velocity命令无法控制住速度问题
已经有11人回复
- 【讨论】lammps中关于fix ave/spatial命令的输出问题
已经有4人回复
- 【求助】lammps中多元化合物的晶格创建问题
已经有14人回复
- 【求助】MS建模导入lammps问题求救
已经有14人回复
- 【求助】Lammps 多机并行的问题
已经有11人回复
w76990480
荣誉版主 (文坛精英)
super-emuch
- 应助: 47 (小学生)
- 贵宾: 13.729
- 金币: 109764
- 散金: 29386
- 红花: 337
- 沙发: 50
- 帖子: 23156
- 在线: 3166.8小时
- 虫号: 942733
- 注册: 2010-01-15
- 专业: 固体力学
- 管辖: 文献求助
2楼2012-03-08 23:15:54
dd5818
银虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 24.7
- 散金: 23
- 帖子: 361
- 在线: 239.7小时
- 虫号: 774387
- 注册: 2009-05-19
- 性别: GG
- 专业: 热力学
3楼2012-03-09 08:37:04
w76990480
荣誉版主 (文坛精英)
super-emuch
- 应助: 47 (小学生)
- 贵宾: 13.729
- 金币: 109764
- 散金: 29386
- 红花: 337
- 沙发: 50
- 帖子: 23156
- 在线: 3166.8小时
- 虫号: 942733
- 注册: 2010-01-15
- 专业: 固体力学
- 管辖: 文献求助
4楼2012-03-09 08:42:33
dd5818
银虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 24.7
- 散金: 23
- 帖子: 361
- 在线: 239.7小时
- 虫号: 774387
- 注册: 2009-05-19
- 性别: GG
- 专业: 热力学
5楼2012-03-09 13:26:14
dd5818
银虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 24.7
- 散金: 23
- 帖子: 361
- 在线: 239.7小时
- 虫号: 774387
- 注册: 2009-05-19
- 性别: GG
- 专业: 热力学
6楼2012-03-09 21:34:29
dd5818
银虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 24.7
- 散金: 23
- 帖子: 361
- 在线: 239.7小时
- 虫号: 774387
- 注册: 2009-05-19
- 性别: GG
- 专业: 热力学
7楼2012-03-09 21:43:44
abinitio
版主 (著名写手)
- 模拟EPI: 2
- 应助: 221 (大学生)
- 贵宾: 0.093
- 金币: 7365.2
- 散金: 3917
- 红花: 57
- 帖子: 1317
- 在线: 485.2小时
- 虫号: 368545
- 注册: 2007-05-11
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
- 管辖: 分子模拟
【答案】应助回帖
★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
chaizhm: 金币+3, 谢谢~ 2012-04-27 19:44:08
dd5818: 金币+3 2012-04-27 20:17:47
感谢参与，应助指数 +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的熔点。 |
8楼2012-04-27 16:03:44
