| 查看: 1963 | 回复: 8 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
奉天大元帅金虫 (小有名气)
|
[求助]
做金属在氧原子环境下钝化模拟，NVT系综下加热，为什么氧原子一点都不动啊！！ 已有2人参与
|
||
|
建模没什么问题，是在盒子中央区域建立了一个球形金属纳米颗粒，周围随机分布氧原子，并在NVT系综下由300K加热到400K，想看金属纳米颗粒的钝化过程。 每次计算后用VMD看结果都是金属原子一直在平衡位置做弛豫运动，这个可以理解，毕竟温度不高，但为什么周围的氧原子一动都不动呢？ 势函数我也是参照论文选取的，究竟是哪里出了问题呢？？？？求大神们帮帮忙吧，稍后贴上in文件，为这事已经快一个礼拜没进展了。 |
回复此楼» 猜你喜欢
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
梦想
已经有9人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有15人回复
-
科研人应该花精力去思考如何解决问题，而不是去凝练问题
已经有14人回复
-
要骂人了，新模版改版就是要淡化问题凝练这种虚的东西，结果有个评委还在说凝练得不够
已经有17人回复
-
基金系统什么内容也没有
已经有10人回复
-
学科评审组评审是指会评吗？
已经有5人回复
-
面上合作单位盖章
已经有9人回复
-
面上函评意见出来了，像什么等级？
已经有17人回复
-
小白求助 投论文要求的highlights应该如何写
已经有4人回复
» 本主题相关商家推荐: (我也要在这里推广)
★ ★
奉天大元帅(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
奉天大元帅
金虫 (小有名气)
- 应助: 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
30