| 查看: 1833 | 回复: 3 | ||||
[求助]
郎之万控温法是怎么实现的？ 已有1人参与
|
| 比如说模拟一个高分子链的运动，用郎之万方程描述每个单元的运动的话涉及到随机积分。很多文献提及用文莱特法算，那这样岂不是很麻烦。还有，怎么实现温度的控制，需要对加和每个单元的动能来得到温度吗？ |
回复此楼» 收录本帖的淘帖专辑推荐
 分子动力学模拟 |
» 猜你喜欢
上海电力大学材料防护与新材料重点实验室招收调剂研究生（材料、化学、电化学，环境）
已经有4人回复
-
材料学求调剂
已经有6人回复
-
303求调剂
已经有5人回复
-
一志愿武理085500机械专业总分300求调剂
已经有7人回复
-
考研调剂
已经有4人回复
-
281求调剂
已经有4人回复
-
0805 316求调剂
已经有6人回复
-
085601求调剂总分293英一数二
已经有3人回复
-
08工学调剂
已经有17人回复
-
340求调剂
已经有4人回复
» 本主题相关商家推荐: (我也要在这里推广)
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
