| 查看: 1871 | 回复: 3 | ||||
[求助]
郎之万控温法是怎么实现的？ 已有1人参与
|
| 比如说模拟一个高分子链的运动，用郎之万方程描述每个单元的运动的话涉及到随机积分。很多文献提及用文莱特法算，那这样岂不是很麻烦。还有，怎么实现温度的控制，需要对加和每个单元的动能来得到温度吗？ |
回复此楼» 收录本帖的淘帖专辑推荐
 分子动力学模拟 |
» 猜你喜欢
297，工科调剂?
已经有4人回复
-
恳请有学校收留
已经有7人回复
-
291求调剂
已经有9人回复
-
300求调剂
已经有11人回复
-
22专硕求调剂
已经有12人回复
-
材料相关专业344求调剂双非工科学校或课题组
已经有25人回复
-
急需调剂
已经有7人回复
-
求调剂
已经有10人回复
-
一志愿华中农业071010，320求调剂
已经有16人回复
-
304求调剂
已经有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
5