| 查看: 2044 | 回复: 3 | ||||
[求助]
郎之万控温法是怎么实现的？ 已有1人参与
|
| 比如说模拟一个高分子链的运动，用郎之万方程描述每个单元的运动的话涉及到随机积分。很多文献提及用文莱特法算，那这样岂不是很麻烦。还有，怎么实现温度的控制，需要对加和每个单元的动能来得到温度吗？ |
回复此楼» 收录本帖的淘帖专辑推荐
 分子动力学模拟 |
» 猜你喜欢
极限拉扯，如何抉择
已经有15人回复
-
申请科研项目纯靠实力能成功吗？
已经有8人回复
-
中科院海洋所招收海洋材料相关领域博士研究生及研究助理
已经有4人回复
-
40+中年人失业后，还能找到好工作吗？
已经有3人回复
-
二审4个多月了还是Awaiting Reviewer Scores
已经有4人回复
-
OpenAI又攻克数学难题了，离AI做有机合成还有多远？
已经有13人回复
-
课题组会帮你通过考核么？
已经有7人回复
-
入职后第一次管课题组的专利缴费，对账才发现官费比别人多掏了一倍
已经有6人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有9人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有10人回复
» 本主题相关商家推荐: (我也要在这里推广)
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
200