| 查看: 1785 | 回复: 8 | ||
奉天大元帅金虫 (小有名气)
|
[求助]
做金属在氧原子环境下钝化模拟，NVT系综下加热，为什么氧原子一点都不动啊！！ 已有2人参与
|
|
建模没什么问题，是在盒子中央区域建立了一个球形金属纳米颗粒，周围随机分布氧原子，并在NVT系综下由300K加热到400K，想看金属纳米颗粒的钝化过程。 每次计算后用VMD看结果都是金属原子一直在平衡位置做弛豫运动，这个可以理解，毕竟温度不高，但为什么周围的氧原子一动都不动呢？ 势函数我也是参照论文选取的，究竟是哪里出了问题呢？？？？求大神们帮帮忙吧，稍后贴上in文件，为这事已经快一个礼拜没进展了。 |
回复此楼» 猜你喜欢
考研英一数一338分
已经有3人回复
-
一志愿085502，267分求调剂
已经有16人回复
-
085801电气专硕272求调剂
已经有3人回复
-
366求调剂
已经有9人回复
-
材料工程085601，270求调剂
已经有37人回复
-
279学硕食品专业求调剂院校
已经有18人回复
-
290调剂生物0860
已经有31人回复
-
一志愿085802 323分求调剂
已经有13人回复
-
277求调剂
已经有23人回复
-
322求调剂，08工科
已经有4人回复
» 本主题相关商家推荐: (我也要在这里推广)
奉天大元帅
金虫 (小有名气)
- 应助: 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
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
★ ★
奉天大元帅(abinitio代发): 金币+2, 感谢应助 2016-05-27 08:48:47
奉天大元帅(abinitio代发): 金币+2, 感谢应助 2016-05-27 08:48:47
|
按我的理解没有O-Al之间以及O-O之间的lj势，O原子肯定是不会动的 发自小木虫IOS客户端 |
6楼2016-05-26 00:37:50
7楼2016-05-26 15:02:16
8楼2016-05-26 17:02:16
9楼2016-05-26 17:03:32
5