版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

返回列表

wolffporter

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 12.8
帖子: 31
在线: 13小时
虫号: 280092
注册: 2006-09-17
专业: 高分子合成化学

【答案】应助回帖

我博士课题期间自学并用matlab写了个MC的代码，不过是用于高分子聚合反应的，参考了杨玉良与张红东写的 “高分子科学中的 Monte Carlo 方法”。化学反应的MC模拟大体来说就是用到两组随机数：r1与r2。

r1是一个1×1的数组，也就是一个数，用于产生每一步反应的时间；
r2是一个N×1的数组，用于选择反应类型，N等于反应物种的数目，也就是说每一个反应物种会分配得到一个随机数值。

设体系有 M 种反应，各反应的反应频率为 ai （i=1,2,...,M）。那么总的反应频率就是 SIGMA（ai）。如果一些反应物种所分配到的 r2 值的范围在（a1,a2]/SIGMA(ai) ，那么这些反应物种选择进行第一种反应；以此类推， r2 值的范围在（a2,a3]/SIGMA(ai) 的那些反应物种，选择进行第二种反应........
反应类型的排序不影响结果，也就是说随你喜欢把哪种反应放第几个。

每当遍历了所有类型的反应，反应时间增加 tao，tao 值由 r1 计算产生 tao=ln(1/r1)/SIGMA(ai)。当时间达到你要的反应时间，终止程序。

赞一下

回复此楼

11楼2016-03-18 14:51:35

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

wolffporter

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 12.8
帖子: 31
在线: 13小时
虫号: 280092
注册: 2006-09-17
专业: 高分子合成化学

引用回帖:

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

我用matlab做的，详见我给本主题楼主的回复，应该在楼上。

赞一下(1人)

回复此楼

12楼2016-03-18 14:55:33

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

wolffporter

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 12.8
帖子: 31
在线: 13小时
虫号: 280092
注册: 2006-09-17
专业: 高分子合成化学

【答案】应助回帖

这是我写的计算流程。

图片1.png

赞一下

回复此楼

13楼2016-03-18 14:58:18

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

吕京毛衣

银虫 (小有名气)

应助: 0 (幼儿园)
金币: 461.3
散金: 7
帖子: 121
在线: 15.1小时
虫号: 3245265
注册: 2014-05-30
专业: 高分子合成化学

引用回帖:

13楼: Originally posted by wolffporter at 2016-03-18 14:58:18
这是我写的计算流程。

图片1.png

嗯嗯非常感谢您的帮助

发自小木虫IOS客户端

赞一下

回复此楼

14楼2016-03-18 22:32:25

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

吕京毛衣

银虫 (小有名气)

应助: 0 (幼儿园)
金币: 461.3
散金: 7
帖子: 121
在线: 15.1小时
虫号: 3245265
注册: 2014-05-30
专业: 高分子合成化学

引用回帖:

11楼: Originally posted by wolffporter at 2016-03-18 14:51:35
我博士课题期间自学并用matlab写了个MC的代码，不过是用于高分子聚合反应的，参考了杨玉良与张红东写的 “高分子科学中的 Monte Carlo 方法”。化学反应的MC模拟大体来说就是用到两组随机数：r1与r2。

r1是一个1 ...

想请教你一下，你开始有没有设定反应体积，如果有要怎么设定呢？

发自小木虫IOS客户端

赞一下

回复此楼

15楼2016-03-31 21:01:12

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