24小时热门版块排行榜    

查看: 4288  |  回复: 17

清泥-3

金虫 (小有名气)

[求助] Kinetic Monte carlo动力学蒙特卡罗模拟化学反应。。 已有3人参与

各位虫友们,有谁对做过动力学蒙特卡罗模拟化学反应? 求指导呀,最近老师让我研究这个,真心看不懂?有哪位大神能够指导下,感激不尽。。。做模拟的时候初始模型是随意给定的还是?如果某个事件能够发生,比如反应A能生成B, 在模拟的时候是不是只在某个格点上发生一次这样的事情,然后进入到下一个循环,还是所有的满足的格子上面的A都能生成B。?
回复此楼

» 收录本帖的淘帖专辑推荐

LAMMPS

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

刚刚开始接触,以后还要继续努力学习。
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

清泥-3

金虫 (小有名气)

怎么都没有人搭理?
刚刚开始接触,以后还要继续努力学习。
2楼2015-04-27 09:55:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

vanceeasleaf

新虫 (初入文坛)

【答案】应助回帖

KMC的code有lammps官方的spparks,以及VTST官方的eon和aKMC,spparks主要用于Lattice KMC,功能比较有限,aKMC是perl写的控制程序,调用vasp用NEB寻找过渡态然后执行KMC,并且是adaptive的,eon是python写的,我比较喜欢这个,功能也是最强的
8楼2015-07-11 13:19:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wolffporter

新虫 (初入文坛)

引用回帖:
10楼: Originally posted by 吕京毛衣 at 2016-03-15 10:12:43
请问针对高分子化学反应的Monte Carlo用什么软件实现

我用matlab做的,详见我给本主题楼主的回复,应该在楼上。
12楼2016-03-18 14:55:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

kjsun

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
清泥-3: 金币+30, ★★★很有帮助 2015-05-04 11:13:10
模拟的时候初始模型是随意给定的还是?如果某个事件能够发生,比如反应A能生成B, 在模拟的时候是不是只在某个格点上发生一次这样的事情,然后进入到下一个循环,还是所有的满足的格子上面的A都能生成B。?

答:初始模型是在满足一些条件的情况下,随便定的。
   一般是:某个格子上的A生成B,然后进入下一循环。




特别提醒:我没做过KMC模拟反应。
我做的是模拟脱附过程,
所以仅供参考。
3楼2015-04-27 10:27:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

清泥-3

金虫 (小有名气)

引用回帖:
3楼: Originally posted by kjsun at 2015-04-27 10:27:22
模拟的时候初始模型是随意给定的还是?如果某个事件能够发生,比如反应A能生成B, 在模拟的时候是不是只在某个格点上发生一次这样的事情,然后进入到下一个循环,还是所有的满足的格子上面的A都能生成B。?

答:初 ...

脱附过程,是用什么算法,BKL,还是第一反应法,还是别的,这些算法的使用是不是看自己的选择,还是说根据自己要模拟什么来定?
刚刚开始接触,以后还要继续努力学习。
4楼2015-04-27 16:14:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

清泥-3

金虫 (小有名气)

引用回帖:
3楼: Originally posted by kjsun at 2015-04-27 10:27:22
模拟的时候初始模型是随意给定的还是?如果某个事件能够发生,比如反应A能生成B, 在模拟的时候是不是只在某个格点上发生一次这样的事情,然后进入到下一个循环,还是所有的满足的格子上面的A都能生成B。?

答:初 ...

非常感谢大神的回答,。。。
刚刚开始接触,以后还要继续努力学习。
5楼2015-04-27 16:15:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kjsun

木虫 (正式写手)

引用回帖:
4楼: Originally posted by 清泥-3 at 2015-04-27 16:14:50
脱附过程,是用什么算法,BKL,还是第一反应法,还是别的,这些算法的使用是不是看自己的选择,还是说根据自己要模拟什么来定?...

我被你吓到了,因为我根本不知道BKL和第一反应法是什么区别。
好吧,我承认我不是这个领域的人,我是外行。

我用KMC模拟的是吸附分子的程序升温脱附(TPD)过程
表面上的吸附的分子有两种行为:扩散和脱附 (忽略再吸附的可能)
扩散远远快于脱附,所以:
1. 让分子扩散,直到扩散平衡。
2. 随机选择一分子,让其脱附,计算时间变化 t= -ln(随机数)/脱附率的加和, 和新的温度值 T=前一步温度+t*升温速度
3 如果T未达到结束温度, 跳转1。否则,中止,结束。
把每个T下的值画出来,就是TPD的结果。


简单到我都不好意思说了。
6楼2015-04-28 12:29:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

清泥-3

金虫 (小有名气)

引用回帖:
6楼: Originally posted by kjsun at 2015-04-28 12:29:01
我被你吓到了,因为我根本不知道BKL和第一反应法是什么区别。
好吧,我承认我不是这个领域的人,我是外行。

我用KMC模拟的是吸附分子的程序升温脱附(TPD)过程
表面上的吸附的分子有两种行为:扩散和脱附 (忽 ...

非常感谢你的回答,你是自己写程序的吗? 能够给一个你的程序代码让我学习下吗?
刚刚开始接触,以后还要继续努力学习。
7楼2015-05-04 11:12:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

清泥-3

金虫 (小有名气)

引用回帖:
8楼: Originally posted by vanceeasleaf at 2015-07-11 13:19:51
KMC的code有lammps官方的spparks,以及VTST官方的eon和aKMC,spparks主要用于Lattice KMC,功能比较有限,aKMC是perl写的控制程序,调用vasp用NEB寻找过渡态然后执行KMC,并且是adaptive的,eon是python写的,我比较喜 ...

你好,你用KMC主要是用来做哪方面的研究?
刚刚开始接触,以后还要继续努力学习。
9楼2015-10-20 15:34:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

吕京毛衣

银虫 (小有名气)

请问针对高分子化学反应的Monte Carlo用什么软件实现
10楼2016-03-15 10:12:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 清泥-3 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 今天系统多次维护,明天很可能放榜! +5 zju2000 2026-08-16 5/250 2026-08-16 23:23 by Haru815
[论文投稿] 岩土工程学报什么时候才能终审完呐 25+3 yeager111 2026-08-10 4/200 2026-08-16 20:42 by tfang
[基金申请] 2027广东省杰青 +3 奶牛小黑 2026-08-15 6/300 2026-08-16 20:07 by 奶牛小黑
[基金申请] 我的国基提前知道中了,可是同事的操作让我实在接受不了,怎么会有这样的人 +11 家与远方 2026-08-10 16/800 2026-08-16 10:28 by ray43
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +10 majunge000 2026-08-11 12/600 2026-08-16 08:18 by xli1984
[基金申请] 小木虫上这么多卖论文的,真有人买论文么?感觉没必要啊 +12 Tide man 2026-08-10 13/650 2026-08-15 16:34 by 氺木
[基金申请] 关于Filecode分析方法 +10 majunge000 2026-08-10 13/650 2026-08-15 12:26 by hanpeng972
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 奇怪,两个人的filecode固定段从头到尾一模一样 +8 布布和一二 2026-08-10 11/550 2026-08-14 14:58 by Equinoxhua
[基金申请] filecode +15 documentary 2026-08-10 17/850 2026-08-14 10:08 by kissu88
[基金申请] FileCode能看出啥? +10 要乐观耀哥 2026-08-10 32/1600 2026-08-14 09:37 by 要乐观耀哥
[基金申请] 静等基金结果 +8 gjjjzhong 2026-08-10 21/1050 2026-08-13 17:56 by 且听虎啸
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[基金申请] 什么时候出结果,有咨询渠道??? +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
信息提示
请填处理意见