| 查看: 1951 | 回复: 8 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
奉天大元帅金虫 (小有名气)
|
[求助]
做金属在氧原子环境下钝化模拟，NVT系综下加热，为什么氧原子一点都不动啊！！ 已有2人参与
|
|
|
建模没什么问题，是在盒子中央区域建立了一个球形金属纳米颗粒，周围随机分布氧原子，并在NVT系综下由300K加热到400K，想看金属纳米颗粒的钝化过程。 每次计算后用VMD看结果都是金属原子一直在平衡位置做弛豫运动，这个可以理解，毕竟温度不高，但为什么周围的氧原子一动都不动呢？ 势函数我也是参照论文选取的，究竟是哪里出了问题呢？？？？求大神们帮帮忙吧，稍后贴上in文件，为这事已经快一个礼拜没进展了。 |
回复此楼» 猜你喜欢
为什么到现在没收到通知？
已经有3人回复
-
有没有仍没收到信息的
已经有5人回复
-
面上意见出来了
已经有10人回复
-
导师吐槽：我怎么摊上了这么个极品研究生！
已经有9人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有5人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有7人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有4人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有10人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有9人回复
» 本主题相关商家推荐: (我也要在这里推广)
abinitio
版主 (著名写手)
- 模拟EPI: 2
- 应助: 221 (大学生)
- 贵宾: 0.093
- 金币: 7365.2
- 散金: 3917
- 红花: 57
- 帖子: 1317
- 在线: 485.2小时
- 虫号: 368545
- 注册: 2007-05-11
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
- 管辖: 分子模拟
5楼2016-05-25 09:41:20
奉天大元帅
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 659.8
- 散金: 13
- 帖子: 76
- 在线: 178.6小时
- 虫号: 3745685
- 注册: 2015-03-17
- 性别: GG
- 专业: 工程热力学
|
# Al dunhua 20nm chiyu echo screen units metal atom_style charge dimension 3 boundary f f f neighbor 2.0 bin neigh_modify every 10 delay 0 check yes # create geometry # create Feball region box block 0 800 0 800 0 800 units box create_box 2 box lattice fcc 3.61 region lvqiu sphere 400 400 400 10 units box create_atoms 1 region lvqiu region lvqiuwai sphere 400 400 400 15 units box side out region yangqiu sphere 400 400 400 30 units box region yangqi intersect 2 yangqiu lvqiuwai create_atoms 2 random 1000 591659 yangqi group yang region yangqi group Al region lvqiu pair_style hybrid/overlay coul/streitz 12.0 wolf 0.31 eam/alloy pair_coeff * * coul/streitz AlO.streitz Al O pair_coeff * * eam/alloy AlO.eam.alloy Al O # compute compute 2 Al pe/atom velocity Al create 300.0 56944 dist gaussian velocity yang create 300.0 591037 dist gaussian fix 1 all nvt temp 300.0 300.0 0.1 #fix 1 all nve fix 2 all wall/reflect xlo EDGE xhi EDGE ylo EDGE yhi EDGE zlo EDGE zhi EDGE fix 3 all qeq/slater 1 12.0 1.0e-6 100 coul/streitz # thermo thermo 1000 thermo_style custom step temp etotal pe ke thermo_modify lost ignore flush yes timestep 0.001 #dump 5 all cfg 1000 dump.h523.*.cfg mass type xs ys zs dump 1 all atom 1000 h523.lammpstrj dump 2 Al custom 1000 h523.*.dump c_2 run 50000 |
2楼2016-05-24 18:31:25
奉天大元帅
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 659.8
- 散金: 13
- 帖子: 76
- 在线: 178.6小时
- 虫号: 3745685
- 注册: 2015-03-17
- 性别: GG
- 专业: 工程热力学
3楼2016-05-24 18:34:34
奉天大元帅
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 659.8
- 散金: 13
- 帖子: 76
- 在线: 178.6小时
- 虫号: 3745685
- 注册: 2015-03-17
- 性别: GG
- 专业: 工程热力学
4楼2016-05-24 18:40:27
投票:
10