24小时热门版块排行榜    

查看: 982  |  回复: 7

chuntao118

新虫 (著名写手)

[求助] Noneqilibrium dynamics和 equilibrium dynamics 有什么差别? 已有1人参与

Noneqilibrium dynamics和 equilibrium dynamics 有什么差别?在进行MD模拟时又有哪些具体的不同?有谁能够具体详细的解释一下,或者提供个比较详细可靠的资料!谢谢!当然,大家也可以发表自己的见解,和研究心得!谢谢!谢谢!金币是不是问题。只要回答的可信,有道理。都会奉上金币!
回复此楼

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lsloneil

专家顾问 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
chuntao118: 金币+20, 有帮助, 有一定的帮助。 2016-01-15 15:49:35

» 本帖已获得的红花(最新10朵)

2楼2016-01-15 11:26:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chuntao118

新虫 (著名写手)

送红花一朵
引用回帖:
2楼: Originally posted by lsloneil at 2016-01-15 11:26:23
http://muchong.com/bbs/viewthread.php?tid=7698181

Noneqilibrium dynamics和 equilibrium dynamics 与系综有什么联系那?是不是说Noneqilibrium dynamics MD 与系综就没有关系?equilibrium dynamics  MD 依赖于系综那?
3楼2016-01-15 15:49:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lsloneil

专家顾问 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
chuntao118: 金币+20, ★★★很有帮助, 很专业,解释的比较到位。 2016-01-16 15:54:50
引用回帖:
3楼: Originally posted by chuntao118 at 2016-01-14 19:49:02
Noneqilibrium dynamics和 equilibrium dynamics 与系综有什么联系那?是不是说Noneqilibrium dynamics MD 与系综就没有关系?equilibrium dynamics  MD 依赖于系综那?...

equilibrium dynamics  MD 依赖于系综那?
是的。只要在平衡态MD里计算的性质都是系综平均。

是不是说Noneqilibrium dynamics MD 与系综就没有关系?
这个不好说。只能说NEMD不是在一个系综里进行采样计算,但其计算的性质还是和之前没有被扰动的系综是有关系的。

举个例子,为了用NEMD计算体系在某个温度下的剪切粘度,对系统添加一个扰动来产生剪切力,进而通过形变计算粘度。整个模拟的过程是非平衡态的(系统一直在演化),但模拟得到的粘度,对应于扰动之前的系统的粘度。很显然,扰动之前的处于平衡态的系统,对应一个系综。

系综这个概念本来就是和平衡态捆绑在一起的,因为根据统计物理的遍历假设,物理量的系综平均=时间平均=平衡态的热力学值。平衡态意味着体系各个热力学性质不随时间变化,而对于随时间演化的非平衡态过程,没有一个与之对应的系综。

有时对于NEMD模拟,会存在所谓的稳态(steady state),比如体系中存在一个稳定的热流,或者粒子流,但体系的其他性质像体积什么的都不变。这种处于稳态的非平衡态,和与之相对应的平衡态的系综有什么关联,是当前非平衡态统计物理研究的一项课题之一, 一个比较典型的例子是处于稳态的self-propelled particle。

» 本帖已获得的红花(最新10朵)

4楼2016-01-16 06:36:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chuntao118

新虫 (著名写手)

送红花一朵
引用回帖:
4楼: Originally posted by lsloneil at 2016-01-16 06:36:14
equilibrium dynamics  MD 依赖于系综那?
是的。只要在平衡态MD里计算的性质都是系综平均。

是不是说Noneqilibrium dynamics MD 与系综就没有关系?
这个不好说。只能说NEMD不是在一个系综里进行采样计算,但 ...

"举个例子,为了用NEMD计算体系在某个温度下的剪切粘度,对系统添加一个扰动来产生剪切力,进而通过形变计算粘度。整个模拟的过程是非平衡态的(系统一直在演化),但模拟得到的粘度,对应于扰动之前的系统的粘度。很显然,扰动之前的处于平衡态的系统,对应一个系综“
回答的很专业,但是还有点问题想请教一下:就一上面您举的例子吧?要是计算体系的在某个温度下的剪切粘度吧?以Poiseuille flow来说吧!在采用力驱动时,应该是NEMD吧(原因:因为体系受外力驱动)?但是要保证体系的温度是某个确定的温度是不是就要采用恒温器来调整温度那?谢谢!
5楼2016-01-16 15:53:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lsloneil

专家顾问 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
chuntao118: 金币+160, ★★★很有帮助, 40 2016-01-17 09:21:19
引用回帖:
5楼: Originally posted by chuntao118 at 2016-01-15 19:53:34
"举个例子,为了用NEMD计算体系在某个温度下的剪切粘度,对系统添加一个扰动来产生剪切力,进而通过形变计算粘度。整个模拟的过程是非平衡态的(系统一直在演化),但模拟得到的粘度,对应于扰动之前的系统的 ...

我没有用NEMD计算粘度的实际操作经验,只是对原理略知一二,以下回答可能不准确。

的确,在lammps用NEMD计算粘度的说明中,当使用fix deform命令施加扰动时,需要同时把体系耦合一个恒温器来防止体系被加热。

类似的,在某些施加外加电场计算电导率的NEMD模拟中,系统也会耦合一个恒温器,来防止外加电场对离子的不断加速产生的升温。

所以,采用恒温器在NEMD模拟中控制温度是说得通的。

需要说明一点,像温度还有压力这样的热力学量,它们的定义是只和平衡态相关的:

对于温度,定义就是平衡态能量对熵的偏导数,在分子模拟中瞬时温度的计算具体体现为对粒子平均动能的计算,粒子动能的系综平均就是温度乘以玻尔兹曼常数。

对于压力,定义就是平衡态能量对体积的偏导数,在分子模拟中瞬时压力的计算具体体现为对应力张量(stress tensor)或者Virial的计算,应力张量或者Virial的系综平均就对应于压力。

对于非平衡态,你没法把不断演化的系统和一个系综联系起来(因为做不到遍历),所以在这种情况下平衡态下定义的温度和压力对于非平衡态严格来讲是失效的。虽然在NEMD模拟中,仍然可以通过计算平均动能和应力张量来计算所谓的瞬时温度和瞬时压力,但这样得到的所谓的压力和温度的物理意义,目前从理论上来讲还不是很清晰,有一些人在做这方面的工作。
6楼2016-01-17 04:58:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chuntao118

新虫 (著名写手)

引用回帖:
6楼: Originally posted by lsloneil at 2016-01-17 04:58:16
我没有用NEMD计算粘度的实际操作经验,只是对原理略知一二,以下回答可能不准确。

的确,在lammps用NEMD计算粘度的说明中,当使用fix deform命令施加扰动时,需要同时把体系耦合一个恒温器来防止体系被加热。
...

十分感谢。我个人也是这样理解的就是怕不对,昨天又查了一大推文献,现在又有您这么详尽的回答和分析!谢谢!遗憾的是我没有把金币全给您,主要是想有其他的人也能够参与进来讨论一下。也许还会有其他的见解!不过还是很感谢您的回答!谢谢!
7楼2016-01-17 09:25:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lvxinyuabc

木虫 (正式写手)

材料大数据

一句话,平衡是通过计算自相关函数;非平衡是通过外力,根据线性响应关系求得物理量。我的论坛(创意杂谈论坛)也曾经发表过相关发面的帖子。
专心科研,分享成长心路历程。关注个人博客www.yffslxy.com,共同鉴证成长之路.
8楼2016-03-02 16:17:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 chuntao118 的主题更新
信息提示
请填处理意见