| 查看: 1979 | 回复: 3 | ||||
[求助]
郎之万控温法是怎么实现的？ 已有1人参与
|
| 比如说模拟一个高分子链的运动，用郎之万方程描述每个单元的运动的话涉及到随机积分。很多文献提及用文莱特法算，那这样岂不是很麻烦。还有，怎么实现温度的控制，需要对加和每个单元的动能来得到温度吗？ |
回复此楼» 收录本帖的淘帖专辑推荐
 分子动力学模拟 |
» 猜你喜欢
帮忙看看fileCode
已经有8人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有3人回复
-
FileCode能看出啥？
已经有15人回复
-
filecode
已经有8人回复
-
确定了，国自然21号放榜
已经有7人回复
-
综述论文作为代表作会不会影响评审专家的印象分？
已经有11人回复
-
关于代码变化问题，想知道的进来
已经有24人回复
-
为什么网上很多人说本周 12号出结果
已经有3人回复
-
静等基金结果
已经有16人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有11人回复
» 本主题相关商家推荐: (我也要在这里推广)
lsloneil
专家顾问 (正式写手)
-
专家经验: +49 - 模拟EPI: 2
- 应助: 148 (高中生)
- 金币: 8640.9
- 散金: 507
- 红花: 30
- 帖子: 505
- 在线: 598.4小时
- 虫号: 556471
- 注册: 2008-05-11
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 分子模拟
【答案】应助回帖
感谢参与，应助指数 +1
|
你如果想比较analytically地分析Langevin equation当然要用到stochastic calculus。但是对于做分子模拟的人来说，需要的只是数值解，Langevin Equation在他们看来和牛顿方程无异，所以依然可以用Verlet算法来进行数值积分。 用Verlet算法积分运动方程，无非就是需要知道粒子的坐标，速度和受力，Langevin equation和牛顿方程唯一不同的地方就是多了个随机力，那就用随机数生成器生成一个力呗。知道了粒子的受力，下一步的速度和坐标用Verlet算法都可以算出来了。 至于Langevin Equation怎么实现温控，核心就在那个随机力里。那个随机力要符合一个高斯分布，那个高斯分布的峰宽（或者说标准差）是精心选择的，其物理依据是涨落耗散定理(Fluctuation-Dissipation Theorem），一旦选好了高斯分布，随机力和Langevin Equation里的粘滞力配合起来，就能得到想要的温度，无须对粒子动能另作额外处理。 |
2楼2015-11-22 15:03:21
3楼2015-11-24 18:36:31
lsloneil
专家顾问 (正式写手)
-
专家经验: +49 - 模拟EPI: 2
- 应助: 148 (高中生)
- 金币: 8640.9
- 散金: 507
- 红花: 30
- 帖子: 505
- 在线: 598.4小时
- 虫号: 556471
- 注册: 2008-05-11
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 分子模拟
4楼2015-11-24 23:09:22
