| 查看: 1974 | 回复: 3 | |||||
[求助]
郎之万控温法是怎么实现的？ 已有1人参与
|
| 比如说模拟一个高分子链的运动，用郎之万方程描述每个单元的运动的话涉及到随机积分。很多文献提及用文莱特法算，那这样岂不是很麻烦。还有，怎么实现温度的控制，需要对加和每个单元的动能来得到温度吗？ |
回复此楼» 收录本帖的淘帖专辑推荐
 分子动力学模拟 |
» 猜你喜欢
filecode
已经有3人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有10人回复
-
静等基金结果
已经有13人回复
-
关于Filecode分析方法
已经有3人回复
-
国自然结果
已经有8人回复
-
FileCode能看出啥？
已经有8人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有22人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有4人回复
-
8月时间戳变的，举个手。玩一下，释放压力
已经有13人回复
-
好奇怪的filecode
已经有5人回复
» 本主题相关商家推荐: (我也要在这里推广)
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
100