24小时热门版块排行榜    

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

xk6891

至尊木虫 (著名写手)

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

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

xk6891

至尊木虫 (著名写手)

引用回帖:
Originally posted by manutdhk at 2010-10-13 23:13:14:


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

不过在run MD的时 ...

前辈,模拟退火就是为了防止原子陷入局部极小化而进一步降低其能量的一种方法吧。那是不是所有的系统做完几何优化之后,都还需要再进行模拟退火的操作呢?
哦,明白各态历经的大概意思了,这个理论企图用少量个体(原子,分子)的力学规律来等同于大系统(同种原子,分子)的统计规律,就是说少量个体的某个物理量的时间平均值等于这个物理量在某一系综下所有状态的平均值。但是统计规律无法等同力学规律,少量个体在模拟的一定时间内也无法跑完所有的状态,所以准各态历经的说法也有可疑之处,更别说各态历经了。
那前辈提到“微观体系(也就是原子的自由度)”看这个原子自由度很多说明书中都提及,不理解它在计算中指的是一个物理量还是什么东西?
人生自古多险阻,勤者自得天酬助。试问否泰何所依,枯藤老枝待新抽。临渊踌躇终迈步,振翅鹏起云霄冲。似是前程甚堪忧,他日振臂揽苍穹。
8楼2010-10-14 00:33:37
已阅   回复此楼   关注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的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 08:14 by C79jjtjAKjEn
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-15 08:02 by C79jjtjAKjEn
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-15 07:50 by C79jjtjAKjEn
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +5 s3fFTmArrBt6 2026-09-14 6/300 2026-09-15 07:38 by C79jjtjAKjEn
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 s3fFTmArrBt6 2026-09-13 4/200 2026-09-15 06:50 by 5BDX0d0WFp7t
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +6 3n8v2C8RimXI 2026-09-13 6/300 2026-09-15 05:50 by 5BDX0d0WFp7t
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 3n8v2C8RimXI 2026-09-13 3/150 2026-09-15 05:49 by 5BDX0d0WFp7t
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +6 3n8v2C8RimXI 2026-09-13 6/300 2026-09-15 05:38 by 5BDX0d0WFp7t
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 3n8v2C8RimXI 2026-09-13 5/250 2026-09-15 05:34 by 5BDX0d0WFp7t
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 3n8v2C8RimXI 2026-09-13 3/150 2026-09-15 05:26 by 5BDX0d0WFp7t
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 QUjhNVAcOSff 2026-09-13 4/200 2026-09-15 05:00 by mibUvS8DDCwf
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 23jxep3nCNZb 2026-09-14 3/150 2026-09-15 01:21 by mibUvS8DDCwf
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 LwdutQ8HoqWP 2026-09-13 6/300 2026-09-14 22:37 by vZfe6xYu34yj
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 QUjhNVAcOSff 2026-09-13 4/200 2026-09-14 21:37 by vZfe6xYu34yj
[考博] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +5 QUjhNVAcOSff 2026-09-13 5/250 2026-09-14 21:36 by vZfe6xYu34yj
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 QUjhNVAcOSff 2026-09-13 5/250 2026-09-14 21:25 by vZfe6xYu34yj
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 QUjhNVAcOSff 2026-09-13 3/150 2026-09-14 15:09 by Aj1rhIDL5ixY
[基金申请] 求教各位大神:2026教育部人文社科青年基金项目何时公示呀? +5 云宴山人 2026-09-10 7/350 2026-09-14 12:11 by ecnu2013
[考博] 西北工业大学材料学院电化学传感与催化课题组招收2027年推免研究生、博士生、博士后 +3 马方园 2026-09-10 4/200 2026-09-14 09:29 by Chem张zz
[考博] 27秋季材料博士申请求助 10+3 Cahal1 2026-09-11 4/200 2026-09-14 08:50 by 北京莱茵编辑
信息提示
请填处理意见