24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 1764  |  回复: 8
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

奉天大元帅

金虫 (小有名气)

[求助] 做金属在氧原子环境下钝化模拟,NVT系综下加热,为什么氧原子一点都不动啊!! 已有2人参与

建模没什么问题,是在盒子中央区域建立了一个球形金属纳米颗粒,周围随机分布氧原子,并在NVT系综下由300K加热到400K,想看金属纳米颗粒的钝化过程。

每次计算后用VMD看结果都是金属原子一直在平衡位置做弛豫运动,这个可以理解,毕竟温度不高,但为什么周围的氧原子一动都不动呢?

势函数我也是参照论文选取的,究竟是哪里出了问题呢????求大神们帮帮忙吧,稍后贴上in文件,为这事已经快一个礼拜没进展了。
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

奉天大元帅

金虫 (小有名气)

这个in文件纯属是为了看 Al和O能不能反应写的一个小试验,结果表明,氧原子根本不动,一点都不像气体,Al原子的位移都比他大得多,O原子基本在原地不动,不会主动和Al原子接触反应,究竟是为什么呢?
3楼2016-05-24 18:34:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 9 个回答

奉天大元帅

金虫 (小有名气)

# 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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

奉天大元帅

金虫 (小有名气)

而且我用了 fix all nvt 300.0 400.0 0.1命令为什么算完后看log文件 发现除了第一步以外的每一步温度都是0呢???
4楼2016-05-24 18:40:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

abinitio

版主 (著名写手)

【答案】应助回帖

感谢参与,应助指数 +1
引用回帖:
3楼: Originally posted by 奉天大元帅 at 2016-05-24 18:34:34
这个in文件纯属是为了看 Al和O能不能反应写的一个小试验,结果表明,氧原子根本不动,一点都不像气体,Al原子的位移都比他大得多,O原子基本在原地不动,不会主动和Al原子接触反应,究竟是为什么呢?

你先不要al,用这个势函数做一个纯o的,看看是不是势的问题?
5楼2016-05-25 09:41:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 291求调剂 +8 Y-cap 2026-03-29 11/550 2026-03-30 21:03 by Y-cap
[考研] 一志愿郑大材料工程290求调剂 +11 Youth_ 2026-03-30 11/550 2026-03-30 21:01 by peike
[考研] 359求调剂 +5 王了个楠 2026-03-25 5/250 2026-03-30 19:36 by 源_2020
[考研] 293求调剂 +3 末未mm 2026-03-30 5/250 2026-03-30 17:23 by 王保杰33
[考研] 材料专硕调剂 +11 椰椰。 2026-03-29 11/550 2026-03-30 16:21 by wangjy2002
[考研] 0703化学/290求调剂/本科经历丰富/工科也可 +13 丹青奶盖 2026-03-26 15/750 2026-03-30 12:35 by fangnagu
[考研] 085602 化学工程专硕 340分求调剂 +4 qianbai11 2026-03-29 4/200 2026-03-30 11:34 by 唐沐儿
[考研] 材料与化工272求调剂 +21 阿斯蒂芬2004 2026-03-28 21/1050 2026-03-30 10:52 by 晴空210210
[考研] 318求调剂 +7 陈晨79 2026-03-30 7/350 2026-03-30 10:49 by 探123
[考研] 348求调剂 +6 小懒虫不懒了 2026-03-28 6/300 2026-03-30 10:29 by Evan_Liu
[考研] 一志愿211,335分,0856,求调剂院校和导师 +7 倾____萧 2026-03-27 8/400 2026-03-30 09:37 by longlotian
[考研] 275求调剂 +15 Micky11223 2026-03-25 20/1000 2026-03-29 20:44 by 唐沐儿
[考研] 352分-085602-一志愿985 +5 海纳百川Ly 2026-03-29 5/250 2026-03-29 09:57 by Sjndkwm
[硕博家园] 招收生物学/细胞生物学调剂 +4 IceGuo 2026-03-26 5/250 2026-03-29 01:25 by griffith2014
[考研] 286求调剂 +4 丢掉懒惰 2026-03-27 7/350 2026-03-28 08:07 by baoball
[考研] 086000调剂 +3 7901117076 2026-03-26 3/150 2026-03-27 21:34 by Jianing_Mi
[考研] 266分求材料化工冶金矿业等专业的调剂 +4 哇呼哼呼哼 2026-03-26 4/200 2026-03-27 17:02 by zhyzzh
[考研] 求调剂 一志愿 本科 北科大 化学 343 +6 13831862839 2026-03-24 7/350 2026-03-26 22:57 by 不吃魚的貓
[考研] 347求调剂 +4 L when 2026-03-25 4/200 2026-03-25 13:37 by cocolv
[考研] 0854人工智能方向招收调剂 +4 章小鱼567 2026-03-24 4/200 2026-03-25 13:29 by 2177681040
信息提示
请填处理意见