24小时热门版块排行榜    

查看: 4264  |  回复: 17

清泥-3

金虫 (小有名气)

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

各位虫友们,有谁对做过动力学蒙特卡罗模拟化学反应? 求指导呀,最近老师让我研究这个,真心看不懂?有哪位大神能够指导下,感激不尽。。。做模拟的时候初始模型是随意给定的还是?如果某个事件能够发生,比如反应A能生成B, 在模拟的时候是不是只在某个格点上发生一次这样的事情,然后进入到下一个循环,还是所有的满足的格子上面的A都能生成B。?
回复此楼
刚刚开始接触,以后还要继续努力学习。
已阅   回复此楼   关注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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] FileCode能看出啥? +7 要乐观耀哥 2026-08-10 16/800 2026-08-11 08:18 by wmfsnow
[基金申请] 帮忙看看fileCode +5 wwncly 2026-08-10 9/450 2026-08-11 08:10 by li9527
[基金申请] 基金中了 +15 laoda193707 2026-08-06 15/750 2026-08-11 00:11 by jiafei2190
[基金申请] filecode +8 documentary 2026-08-10 8/400 2026-08-10 19:20 by loufangrui
[基金申请] 确定了,国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分? +10 yufeiwaner 2026-08-09 11/550 2026-08-10 18:47 by yufeiwaner
[基金申请] 关于代码变化问题,想知道的进来 +17 且听虎啸 2026-08-07 24/1200 2026-08-10 18:35 by zhangduo2008
[基金申请] 我的国基提前知道中了,可是同事的操作让我实在接受不了,怎么会有这样的人 +8 家与远方 2026-08-10 11/550 2026-08-10 17:09 by 六两废铜
[基金申请] 奇怪,两个人的filecode固定段从头到尾一模一样 +4 布布和一二 2026-08-10 4/200 2026-08-10 14:50 by 医学老男孩
[基金申请] 好奇怪的filecode +4 布布和一二 2026-08-08 5/250 2026-08-10 14:20 by 冰心玉壶晴
[基金申请] 2026国自然放榜时间 +9 布布和一二 2026-08-08 9/450 2026-08-10 11:22 by xxxx2020
[基金申请] 这样的filecode谁见过 +11 布布和一二 2026-08-08 22/1100 2026-08-10 11:10 by wmfsnow
[基金申请] filecode与中标关系的预测 +5 布布和一二 2026-08-07 5/250 2026-08-09 16:15 by 袁向阳007
[基金申请] fileCode有新解读? +10 Tide man 2026-08-08 18/900 2026-08-09 12:55 by 仁砚薪传
[基金申请] 化学口download_prp&fileCode的固定段好像这几天一直没变,有变的大神么? +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] 固定端突然变了,今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] 听说今天filecode变了 +24 布布和一二 2026-08-06 47/2350 2026-08-07 16:02 by zhiyanjiang
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
[基金申请] 求各位大神看下 100+6 hpkpkpkp 2026-08-05 33/1650 2026-08-06 14:49 by zhiyanjiang
[基金申请] 有没有H口的?有收到消息的吗? +3 超级海虾 2026-08-04 3/150 2026-08-04 17:26 by 学教育滴
信息提示
请填处理意见