| 查看: 3688 | 回复: 10 | |||
| 本帖产生 1 个 模拟EPI ，点击这里进行查看 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
xk6891至尊木虫 (著名写手)
|
[交流]
【讨论】关于md中系统稳定的问题 已有3人参与
|
||
|
小虫初涉md，问一个弱弱的问题 作完结构优化之后，看很多网上的资料都要求系统在热浴的环境中（是不是准确的应该说是在一定的系综下）跑一会儿，意思是说希望系统能在这个环境中稳定，如何判断系统在这个环境中是否达到了稳定？稳定之后，在做以下的性质计算时，是否还需要这个系综（譬如Nose Hoover热浴）的关键词？ 急盼各位高手拍砖，谢谢  |
回复此楼» 猜你喜欢
人气不行了
已经有10人回复
-
如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》
已经有6人回复
-
我面上完蛋了
已经有10人回复
-
能否退出参与的面上项目解除限项
已经有24人回复
-
filecode，4个jtjc了
已经有17人回复
-
没有任何消息-是不是就凉了
已经有6人回复
-
2026国自然函评费到账
已经有19人回复
-
建议基金发布提前给出明确的时间点
已经有16人回复
-
今日不放榜？网传国自然预计 8 月 27 日可查结果
已经有21人回复
-
估计是周四
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 分子模拟时间越长越好吗？
已经有44人回复
- 利用siesta做verlet MD时出现问题siesta: WARNING: Atoms
已经有3人回复
- 【求助】关于玻尔兹曼公式
已经有29人回复
manutdhk
木虫 (小有名气)
- 模拟EPI: 1
- 应助: 5 (幼儿园)
- 金币: 5939.4
- 散金: 70
- 红花: 5
- 帖子: 94
- 在线: 97.3小时
- 虫号: 983544
- 注册: 2010-03-27
- 专业: 核技术及其应用
★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
ghcacj(金币+4):谢谢 2010-10-13 17:02:14
xk6891(金币+5): 2010-10-13 18:18:06
小木虫(金币+0.5):给个红包，谢谢回帖交流
ghcacj(金币+4):谢谢 2010-10-13 17:02:14
xk6891(金币+5): 2010-10-13 18:18:06
|
首先，不知道你要做什么性质的计算，平衡态还是非平衡态？ 如果是平衡态的计算，所谓的稳定，确切地说应该是 热力学平衡。对于一个热力学系统，平衡是指 系统达到 吉布斯分布。NVT系综下，系统的温度与 热浴温度一致，能量在一定平衡值附近振荡……等等。 至于计算平衡态性质，是否需要继续保持该系综，理论上来说 其实问题不大。根据统计物理学原理，NVE系综和NVT系综，在热力学极限下，得到的结果是相一致的。 |
3楼2010-10-13 15:36:56
xk6891
至尊木虫 (著名写手)
- 模拟EPI: 1
- 应助: 0 (幼儿园)
- 金币: 13735.8
- 散金: 379
- 红花: 9
- 帖子: 1514
- 在线: 546小时
- 虫号: 890342
- 注册: 2009-11-01
- 专业: 理论和计算化学
不好意思，自己顶一下 2楼2010-10-13 13:12:04
xk6891
至尊木虫 (著名写手)
- 模拟EPI: 1
- 应助: 0 (幼儿园)
- 金币: 13735.8
- 散金: 379
- 红花: 9
- 帖子: 1514
- 在线: 546小时
- 虫号: 890342
- 注册: 2009-11-01
- 专业: 理论和计算化学
★ ★ ★
zh1987hs(金币+3):鼓励交流 2010-10-13 18:59:19
zh1987hs(金币+3):鼓励交流 2010-10-13 18:59:19
|
先谢谢指导，呵呵，提的问题太不专业了。是平衡态，这一点应该说明。可以这样说么，譬如NVT系综，md结果当中ions的温度和electrons的动能应该在nose hoover设定的温度和动能值附近震荡，且粒子数目和所计算系统体积（晶胞大小）不变，这就可以说这个系统符合了NVT的系综，这个系统就稳定了。而如果是NPT系综就看输出文件中系统的压力和温度值的变化，如果其在设定值周围震荡，就可以说系统平衡了。 同时由于系统已经平衡了，可以认为系统已经满足了热力学极限的条件，此时在什么系综下已经不重要了，因为各个系综的情况已经相同。按照在某一系综下稳定的系统，再计算性质，加不加关键词以不重要。 不知这样理解是否正确？还请赐教  |
4楼2010-10-13 18:40:35
manutdhk
木虫 (小有名气)
- 模拟EPI: 1
- 应助: 5 (幼儿园)
- 金币: 5939.4
- 散金: 70
- 红花: 5
- 帖子: 94
- 在线: 97.3小时
- 虫号: 983544
- 注册: 2010-03-27
- 专业: 核技术及其应用
★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zh1987hs(金币+4):谢谢 2010-10-13 18:59:32
xk6891(金币+5):再次感谢 2010-10-13 22:49:12
小木虫(金币+0.5):给个红包，谢谢回帖交流
zh1987hs(金币+4):谢谢 2010-10-13 18:59:32
xk6891(金币+5):再次感谢 2010-10-13 22:49:12
|
对，大体是这样的。 有点小错误：NVT系综，要看 压强随时间的振荡； NPT系综，要看体积随时间的振荡。 我上个帖子提到的 热力学平衡的判据 只是理论上的。 在MD 过程中，很难跟踪整个体系是不是处于 吉布斯分布。而温度，能量的判据是常用的。因为温度或者能量，原本就是吉布斯分布中的参数，或者说是 可观测量，或者说是 一阶矩。更细致一点的，可能还有 二阶矩，譬如 能量的振荡幅度，对应的是比热之类的物理量。 以上这些判据，我个人认为，实质上就是将 吉布斯分布做一个泰勒式的展开 的函数逼近。 一阶矩对应这个展开的零阶项，二阶矩对应的是二阶项……等等。 |
5楼2010-10-13 18:53:22
20