24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 1767  |  回复: 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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[有机交流] 甲基亚磺磺酸钠和甲基磺酸酯反应机理 10+3 kaobao456 2026-03-29 4/200 2026-03-30 23:16 by nBu锂
[考研] 083000学硕274求调剂 +9 Li李鱼 2026-03-26 9/450 2026-03-30 21:06 by dophin1985
[考研] 0703化学321分求调剂 +10 三dd. 2026-03-30 11/550 2026-03-30 19:24 by markhwc
[考研] 329求调剂 +8 星野? 2026-03-26 8/400 2026-03-30 13:41 by chemdavid
[考研] 考研调剂 +7 小蜡新笔 2026-03-29 7/350 2026-03-29 19:00 by 学员8dgXkO
[考研] 329求调剂 +10 钮恩雪 2026-03-25 10/500 2026-03-29 13:32 by peike
[考研] 356求调剂 +4 gysy?s?a 2026-03-28 4/200 2026-03-29 10:32 by 唐沐儿
[考研] 356求调剂 +3 gysy?s?a 2026-03-28 3/150 2026-03-29 00:33 by 544594351
[考研] 315求调剂 +4 akie... 2026-03-28 5/250 2026-03-28 21:05 by zhq0425
[考研] 071000生物学求调剂,初试成绩343 +7 小小甜面团 2026-03-25 7/350 2026-03-28 20:25 by 唐沐儿
[考研] 085602 化工专硕 338分 求调剂 +12 路痴小琪 2026-03-27 12/600 2026-03-28 15:41 by L135790
[考研] 286求调剂 +4 丢掉懒惰 2026-03-27 7/350 2026-03-28 08:07 by baoball
[考研] 265求调剂 +8 小木虫085600 2026-03-27 8/400 2026-03-27 22:16 by 无际的草原
[考研] 266分求材料化工冶金矿业等专业的调剂 +4 哇呼哼呼哼 2026-03-26 4/200 2026-03-27 17:02 by zhyzzh
[考博] 26申博 +3 加油冲啊! 2026-03-26 3/150 2026-03-27 15:38 by cls512
[考研] 考研调剂 +9 小蜡新笔 2026-03-26 9/450 2026-03-27 11:10 by 不吃魚的貓
[考研] 材料求调剂 +5 .m.. 2026-03-25 5/250 2026-03-27 11:08 by 不吃魚的貓
[论文投稿] Journal of Mechanical Science and Technology +3 Russ_ss 2026-03-25 5/250 2026-03-27 10:49 by 陆小果画大饼
[考研] 085602化学工程求调剂。 +4 平乐乐乐 2026-03-26 4/200 2026-03-26 17:57 by fmesaito
[考研] 考研一志愿苏州大学初始315(英一)求调剂 +3 sbdksD 2026-03-24 4/200 2026-03-25 18:16 by xcjcqu
信息提示
请填处理意见