查看: 3845  |  回复: 10
本帖产生 1 个 模拟EPI ,点击这里进行查看
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

xk6891

至尊木虫 (著名写手)

[交流] 【讨论】关于md中系统稳定的问题 已有3人参与

小虫初涉md,问一个弱弱的问题
作完结构优化之后,看很多网上的资料都要求系统在热浴的环境中(是不是准确的应该说是在一定的系综下)跑一会儿,意思是说希望系统能在这个环境中稳定,如何判断系统在这个环境中是否达到了稳定?稳定之后,在做以下的性质计算时,是否还需要这个系综(譬如Nose Hoover热浴)的关键词?
急盼各位高手拍砖,谢谢
回复此楼
人生自古多险阻,勤者自得天酬助。试问否泰何所依,枯藤老枝待新抽。临渊踌躇终迈步,振翅鹏起云霄冲。似是前程甚堪忧,他日振臂揽苍穹。
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

manutdhk

木虫 (小有名气)

★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
zh1987hs(金币+3): 谢谢 2010-10-13 23:19:52
xk6891(金币+5): 2010-10-14 00:34:40
引用回帖:
Originally posted by xk6891 at 2010-10-13 22:48:18:

系统跑md的过程就是一个各态历经的过程?
热统的知识忘的太多了,还要补起来,我自己这两天再看看。
再次谢谢前辈指教,感激不尽。学到很多知识。

客气了。。。
在达到平衡以后 是各态历经。各态历经只是一个假说,严格地说并不成立。只是目前还没有证据证明不成立。有人提出,系统有可能会陷入某些局部极小附近出不来,之类的说辞。

不过在run MD的时候,就不管这个了,反正系综平均等价于时间平均就可以了。

多说两句:统计物理学的基础是 微正则系综假设,也就是说 微观的状态点是等几率的。而在热力学极限下,宏观的观测量等同于 对应的微观量的系综平均。MD和MC都是基于这样一个原则来实现的。当然,对于MD而言,还有一个各态历经假说。 如果从微观体系(也就是原子的自由度)得出宏观量,可以说是 所有做计算模拟的人就是靠这个吃饭的。

[ Last edited by manutdhk on 2010-10-13 at 23:33 ]
7楼2010-10-13 23:13:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 11 个回答

xk6891

至尊木虫 (著名写手)

不好意思,自己顶一下
人生自古多险阻,勤者自得天酬助。试问否泰何所依,枯藤老枝待新抽。临渊踌躇终迈步,振翅鹏起云霄冲。似是前程甚堪忧,他日振臂揽苍穹。
2楼2010-10-13 13:12:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

manutdhk

木虫 (小有名气)

★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
ghcacj(金币+4):谢谢 2010-10-13 17:02:14
xk6891(金币+5): 2010-10-13 18:18:06
引用回帖:
Originally posted by xk6891 at 2010-10-12 15:51:15:
小虫初涉md,问一个弱弱的问题
作完结构优化之后,看很多网上的资料都要求系统在热浴的环境中(是不是准确的应该说是在一定的系综下)跑一会儿,意思是说希望系统能在这个环境中稳定,如何判断系统在这个环境中是 ...

首先,不知道你要做什么性质的计算,平衡态还是非平衡态?
如果是平衡态的计算,所谓的稳定,确切地说应该是 热力学平衡。对于一个热力学系统,平衡是指 系统达到 吉布斯分布。NVT系综下,系统的温度与 热浴温度一致,能量在一定平衡值附近振荡……等等。

至于计算平衡态性质,是否需要继续保持该系综,理论上来说 其实问题不大。根据统计物理学原理,NVE系综和NVT系综,在热力学极限下,得到的结果是相一致的。
3楼2010-10-13 15:36:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xk6891

至尊木虫 (著名写手)

★ ★ ★
zh1987hs(金币+3):鼓励交流 2010-10-13 18:59:19
引用回帖:
Originally posted by manutdhk at 2010-10-13 15:36:56:


首先,不知道你要做什么性质的计算,平衡态还是非平衡态?
如果是平衡态的计算,所谓的稳定,确切地说应该是 热力学平衡。对于一个热力学系统,平衡是指 系统达到 吉布斯分布。NVT系综下,系统的温度与 热浴 ...

先谢谢指导,呵呵,提的问题太不专业了。是平衡态,这一点应该说明。可以这样说么,譬如NVT系综,md结果当中ions的温度和electrons的动能应该在nose hoover设定的温度和动能值附近震荡,且粒子数目和所计算系统体积(晶胞大小)不变,这就可以说这个系统符合了NVT的系综,这个系统就稳定了。而如果是NPT系综就看输出文件中系统的压力和温度值的变化,如果其在设定值周围震荡,就可以说系统平衡了。
同时由于系统已经平衡了,可以认为系统已经满足了热力学极限的条件,此时在什么系综下已经不重要了,因为各个系综的情况已经相同。按照在某一系综下稳定的系统,再计算性质,加不加关键词以不重要。
不知这样理解是否正确?还请赐教
人生自古多险阻,勤者自得天酬助。试问否泰何所依,枯藤老枝待新抽。临渊踌躇终迈步,振翅鹏起云霄冲。似是前程甚堪忧,他日振臂揽苍穹。
4楼2010-10-13 18:40:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 课题组会帮你通过考核么? +3 好运来~9527 2026-10-10 6/300 2026-10-12 14:55 by 简单化xn
[有机交流] OpenAI又攻克数学难题了,离AI做有机合成还有多远? +4 wangyl_123 2026-10-10 7/350 2026-10-12 14:47 by chuans0923
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:34 by XQNyDCa1Tk5n
[找工作] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 22:26 by XQNyDCa1Tk5n
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:26 by XQNyDCa1Tk5n
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:15 by XQNyDCa1Tk5n
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 22:11 by XQNyDCa1Tk5n
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:06 by XQNyDCa1Tk5n
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:06 by XQNyDCa1Tk5n
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 14:25 by VOMexwnDzOFZ
[教师之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 (EPI+-1)(金币-50) +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 14:05 by VOMexwnDzOFZ
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:51 by VOMexwnDzOFZ
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:45 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:54 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:34 by VOMexwnDzOFZ
[考博] 招收博士生一名 +3 mandolin 2026-10-10 3/150 2026-10-11 10:21 by 不是山谷!
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 qR2AdCWTVtkw 2026-10-09 4/200 2026-10-11 03:14 by M0zrYVSAnxpE
[有机交流] 各位大佬,路线设计求助 +6 zapen 2026-10-08 8/400 2026-10-10 15:59 by zapen
[基金申请] 入职后第一次管课题组的专利缴费,对账才发现官费比别人多掏了一倍 +5 13108017953 2026-10-08 5/250 2026-10-10 11:15 by wdwd123
[有机交流] 我应该选哪条合成路线让实验成功率高一些? +5 jing816y 2026-10-08 9/450 2026-10-10 09:05 by zapen
信息提示
请填处理意见