24小时热门版块排行榜    

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

奉天大元帅

金虫 (小有名气)

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

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

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

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

S201405010

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by 奉天大元帅 at 2016-05-24 18:31:25
# 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 ...

你好,请教一下,我做有机分子的动力学模拟,观察X原子在有机分子周围的迁移情况,想得到X原子迁移位点图(大概每隔10个点取一个坐标),类似下面的图,请问你有没有小的脚本程序
做金属在氧原子环境下钝化模拟,NVT系综下加热,为什么氧原子一点都不动啊!!
图片2.png

9楼2016-05-26 17:03:32
已阅   回复此楼   关注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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 投票:  有多少人是今天查系统知道结果的? +16 爱看书的可乐 2026-08-26 18/900 2026-08-28 17:09 by winsaint
[基金申请] 有没有仍没收到信息的 +3 德尚中行 2026-08-27 3/150 2026-08-28 15:24 by rongshuxia
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +8 苏东坡二世 2026-08-23 8/400 2026-08-28 14:41 by lkforward
[基金申请] 基金不中,共勉 +11 eulota 2026-08-26 11/550 2026-08-28 14:22 by 火星超人xi
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +4 gdd2018 2026-08-28 9/450 2026-08-28 14:22 by jklily
[基金申请] 2026年叶企孙基金 +3 bud_bud 2026-08-27 6/300 2026-08-28 11:53 by bud_bud
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +12 晴天加油 2026-08-26 13/650 2026-08-28 11:32 by skyjudy
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +5 yuleib84 2026-08-26 6/300 2026-08-28 00:02 by yudaoqian88
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 申请删除本帖 +6 lyz123lyz 2026-08-27 7/350 2026-08-27 17:31 by 宁静致远sy
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 国合现在查不到了吗? +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
信息提示
请填处理意见