24小时热门版块排行榜    

查看: 1957  |  回复: 8
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

奉天大元帅

金虫 (小有名气)

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

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

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

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

» 猜你喜欢

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

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

S201405010

新虫 (初入文坛)

你好,请教一下,我做有机分子的动力学模拟,观察X原子在有机分子周围的迁移情况,想得到X原子迁移位点图(大概每隔10个点取一个坐标),请问你有没有小的脚本程序
8楼2016-05-26 17:02:16
已阅   回复此楼   关注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的回帖

奉天大元帅

金虫 (小有名气)

这个in文件纯属是为了看 Al和O能不能反应写的一个小试验,结果表明,氧原子根本不动,一点都不像气体,Al原子的位移都比他大得多,O原子基本在原地不动,不会主动和Al原子接触反应,究竟是为什么呢?
3楼2016-05-24 18:34:34
已阅   回复此楼   关注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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +6 yuleib84 2026-08-26 7/350 2026-08-31 19:46 by 鱼翔浅底1
[基金申请] 面上意见出来了 +12 黄鸟于飞Chao 2026-08-29 23/1150 2026-08-31 18:57 by 黄鸟于飞Chao
[基金申请] 中青基了要发朋友圈吗? +7 349506619 2026-08-28 7/350 2026-08-31 13:39 by 冼亮淀粉酶
[基金申请] 投票:  有多少人是今天查系统知道结果的? +17 爱看书的可乐 2026-08-26 19/950 2026-08-31 11:51 by xiangy672
[基金申请] 29号明天会评吗 +4 笨笨唐 2026-08-28 4/200 2026-08-31 09:30 by huixian257
[基金申请] 为什么到现在没收到通知? +5 tannykie 2026-08-29 5/250 2026-08-30 21:05 by purplejack
[考博] 找导师 +6 yuanjiabao 2026-08-29 7/350 2026-08-30 14:40 by 生科新手
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +12 Kittylucky 2026-08-27 13/650 2026-08-29 00:04 by superceng
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 基金未中,这种答复是模板吗? +5 zhaosm1982 2026-08-27 6/300 2026-08-27 16:00 by lfy8008
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
信息提示
请填处理意见