24小时热门版块排行榜    

查看: 1840  |  回复: 3
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

我是飞锅

新虫 (初入文坛)

[求助] 郎之万控温法是怎么实现的? 已有1人参与

比如说模拟一个高分子链的运动,用郎之万方程描述每个单元的运动的话涉及到随机积分。很多文献提及用文莱特法算,那这样岂不是很麻烦。还有,怎么实现温度的控制,需要对加和每个单元的动能来得到温度吗?
回复此楼

» 收录本帖的淘帖专辑推荐

分子动力学模拟

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

致力于搞科研,就是这么简单
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

我是飞锅

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by lsloneil at 2015-11-22 15:03:21
你如果想比较analytically地分析Langevin equation当然要用到stochastic calculus。但是对于做分子模拟的人来说,需要的只是数值解,Langevin Equation在他们看来和牛顿方程无异,所以依然可以用Verlet算法来进行数 ...

首先,谢谢了。我觉得郎之万控温没有这么简单,并不是直接把随机力加到确定力(比如说LJ力和键力)上就可以了,因为按涨落耗散定理得出的随机力在量级上和确定力要小很多,随机力几乎对整体的受力不起作用。
致力于搞科研,就是这么简单
3楼2015-11-24 18:36:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 4 个回答

lsloneil

专家顾问 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lsloneil

专家顾问 (正式写手)

引用回帖:
3楼: Originally posted by 我是飞锅 at 2015-11-23 22:36:31
首先,谢谢了。我觉得郎之万控温没有这么简单,并不是直接把随机力加到确定力(比如说LJ力和键力)上就可以了,因为按涨落耗散定理得出的随机力在量级上和确定力要小很多,随机力几乎对整体的受力不起作用。...

how do you know "按涨落耗散定理得出的随机力在量级上和确定力要小很多"? Do you have the data?

Could you give an estimate of the random force as well as the viscous force in your system?
4楼2015-11-24 23:09:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 上海电力大学材料防护与新材料重点实验室招收调剂研究生(材料、化学、电化学,环境) +4 我爱学电池 2026-03-23 4/200 2026-03-25 00:59 by 1027_324
[考研] 085601求调剂总分293英一数二 +3 钢铁大炮 2026-03-24 3/150 2026-03-24 22:03 by bingxueer79
[考研] 化工专硕求调剂 +3 question挽风 2026-03-24 3/150 2026-03-24 18:48 by jhhcooi
[考研] 307求调剂 +3 余意卿 2026-03-21 6/300 2026-03-24 15:03 by 余意卿
[考研] 081700 调剂 267分 +9 迷人的哈哈 2026-03-23 9/450 2026-03-24 11:58 by 544594351
[考研] 384求调剂 +3 子系博 2026-03-22 6/300 2026-03-23 21:45 by 子系博
[考研] 291求调剂 +5 孅華 2026-03-22 5/250 2026-03-23 09:20 by haoshis
[考研] 求调剂一志愿海大,0703化学学硕304分,有大创项目,四级已过 +6 幸运哩哩 2026-03-22 10/500 2026-03-22 20:10 by edmund7
[考研] 一志愿 西北大学 ,070300化学学硕,总分287,双非一本,求调剂。 +3 晨昏线与星海 2026-03-20 3/150 2026-03-22 16:00 by ColorlessPI
[考研] 考研调剂 +3 呼呼?~+123456 2026-03-21 3/150 2026-03-21 20:04 by 无际的草原
[考研] 一志愿东华大学控制学硕320求调剂 +3 Grand777 2026-03-21 3/150 2026-03-21 19:23 by 简之-
[考研] 0703化学297求调剂 +3 Daisy☆ 2026-03-20 3/150 2026-03-21 17:45 by ColorlessPI
[考研] 材料 271求调剂 +5 展信悦_ 2026-03-21 5/250 2026-03-21 17:29 by 学员8dgXkO
[考研] 一志愿华南师大 070300(化学)304分求调剂 +3 0703武芊慧雪304 2026-03-18 3/150 2026-03-21 00:48 by JourneyLucky
[考研] 中南大学化学学硕337求调剂 +3 niko- 2026-03-19 6/300 2026-03-20 21:58 by luoyongfeng
[考研] A区线材料学调剂 +5 周周无极 2026-03-20 5/250 2026-03-20 21:33 by laoshidan
[考研] 295复试调剂 +8 简木ChuFront 2026-03-19 8/400 2026-03-20 20:44 by zhukairuo
[考研] 295材料求调剂,一志愿武汉理工085601专硕 +5 Charlieyq 2026-03-19 5/250 2026-03-20 20:35 by JourneyLucky
[考研] 本科郑州大学物理学院,一志愿华科070200学硕,346求调剂 +4 我不是一根葱 2026-03-18 4/200 2026-03-19 09:11 by 浮云166
[考研] 收复试调剂生 +4 雨后秋荷 2026-03-18 4/200 2026-03-18 14:16 by elevennnne
信息提示
请填处理意见