| 查看: 1378 | 回复: 14 | |||
| 当前主题已经存档。 | |||
[交流]
【求助】求一段monte carlo 程序 已有1人参与
|
|||
|
求一段monte carlo 程序，不管是实现什么的，我主要是想看看这个运算法则的样子，以及写法 [ Last edited by wuli8 on 2009-11-13 at 21:40 ] |
回复此楼» 猜你喜欢
081700，311，求调剂
已经有15人回复
-
一志愿北京化工085600 310分求调剂
已经有18人回复
-
085600材料与化工301分求调剂院校
已经有4人回复
-
材料与化工371求调剂
已经有14人回复
-
336材料与化工085600求调剂
已经有7人回复
-
材料334求调剂
已经有18人回复
-
331求调剂
已经有8人回复
-
332求调剂
已经有17人回复
-
一志愿南京航空航天大学 材料与化工329分求调剂
已经有4人回复
-
材料专硕322
已经有7人回复
pfzhangphy
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1607.1
- 帖子: 61
- 在线: 173.4小时
- 虫号: 832213
- 注册: 2009-08-22
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
2楼2009-11-16 14:28:20
bay__gulf
金虫 (著名写手)
刘苏州
- 模拟EPI: 8
- 应助: 9 (幼儿园)
- 贵宾: 4.85
- 金币: 2332.8
- 红花: 1
- 帖子: 1344
- 在线: 271小时
- 虫号: 592012
- 注册: 2008-09-03
- 专业: 理论和计算化学
- 管辖: 分子模拟
★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zeoliters(金币+4,VIP+0):多谢回帖交流！ 11-17 15:04
小木虫(金币+0.5):给个红包，谢谢回帖交流
zeoliters(金币+4,VIP+0):多谢回帖交流！ 11-17 15:04
|
MC方法很好,只是形式复杂的力场不好加进去 当初学习MC时候写一个例子 用于逃离曲线上 local minimum, 找到的global minimum 实验函数y = cos (x) + cos (3x), 图像见图top 最低点-2.0, 收敛目标-1.99999 代码 需要手工输入三个参数, a: x初值, d:尝试移动的最大范围, t: 相当于kT 图middle 和buttom 分别在(a,d,t)=(0.2, 0.1, 0.8) 时候, x 和y 随运行步数的演化 可见在有限的时间内, 通过简单的算法可以逃离曲线上 local minimum, 找到的global minimum 适合分子复杂, 不能一一列举的各个构象时候的构象分析 |
3楼2009-11-16 14:58:51
wuchenwf
荣誉版主 (职业作家)
- 应助: 0 (幼儿园)
- 贵宾: 3.433
- 金币: 19419.2
- 散金: 10
- 红花: 4
- 帖子: 3560
- 在线: 1035.7小时
- 虫号: 398569
- 注册: 2007-06-10
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
- 管辖: 第一性原理
4楼2009-11-16 21:52:01
zhiping780
木虫 (著名写手)
天鹰帮帮主
- 应助: 0 (幼儿园)
- 贵宾: 0.007
- 金币: 1584.7
- 散金: 20
- 沙发: 2
- 帖子: 1396
- 在线: 9.7小时
- 虫号: 589991
- 注册: 2008-08-30
- 性别: GG
- 专业: 工业催化
5楼2009-11-22 19:25:54
6楼2010-03-07 00:13:07
bay__gulf
金虫 (著名写手)
刘苏州
- 模拟EPI: 8
- 应助: 9 (幼儿园)
- 贵宾: 4.85
- 金币: 2332.8
- 红花: 1
- 帖子: 1344
- 在线: 271小时
- 虫号: 592012
- 注册: 2008-09-03
- 专业: 理论和计算化学
- 管辖: 分子模拟
7楼2010-03-07 11:25:40
8楼2010-03-07 14:11:12
bay__gulf
金虫 (著名写手)
刘苏州
- 模拟EPI: 8
- 应助: 9 (幼儿园)
- 贵宾: 4.85
- 金币: 2332.8
- 红花: 1
- 帖子: 1344
- 在线: 271小时
- 虫号: 592012
- 注册: 2008-09-03
- 专业: 理论和计算化学
- 管辖: 分子模拟
★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
lei0736(金币+2):谢谢 2010-03-07 16:18
ghcacj(金币+2):谢谢 2010-03-09 09:20
小木虫(金币+0.5):给个红包，谢谢回帖交流
lei0736(金币+2):谢谢 2010-03-07 16:18
ghcacj(金币+2):谢谢 2010-03-09 09:20
|
是你把我问住了 -- 话说我只是学了点MC的皮毛,从未真正用过MC做化学方面的研究 就这几点我说一下我的看法吧 1 关于力场 随着计算能力的的告诉发展,研究的体系也越来越复杂. 比如模拟膜蛋白,就要把蛋白的,水的,lipid uint的,以及水中离子 就连那个常常被人抱怨参数不全的charmm力场,其原子类型也有200余个 键/角/二面角的类型更是200的几次方的数目 而据我所致,做MC的大多是自编程序.通用软件的数目也非常不乐观 要把这样复杂的力场写到自己的程序里面,这样的通才数目有限,不能普及 2 问题好像不能这样问,缺了哪一点都不行 3 metropolis 取样实现了波尔兹曼分布,所以才在分子模拟中有广泛应用吧 |
9楼2010-03-07 14:44:11
10楼2010-03-08 21:01:41
