计算化学公社

标题: lammps升温问题求教 [打印本页]

作者
Author: 灰天皮匠 时间: 2019-4-2 11:51
标题: lammps升温问题求教
本帖最后由灰天皮匠于 2019-4-2 11:52 编辑

接触分子模拟时间短，基础比较差，还在学习中。有一个问题想请教各位老师，如果一个金属体系在0K时处于优化好的稳定状态，考虑它在某个温度的热稳定性时：
方法一是直接建立一个高温（例如1200K）的速度分布，然后恒温至1200K；方法二从0K逐步升温至1200K，这两个过程区别在哪里？

我在别的论坛中看到升温退火过程，是要找到体系最稳定的状态，是不是第二种方法需要先做升温退火，第一种方法不用？

作者
Author: 少年爱吃地三鲜 时间: 2019-4-2 19:10
我理解的退火过程就是逐步的降温升温

作者
Author: 灰天皮匠 时间: 2019-4-2 21:00

少年爱吃地三鲜发表于 2019-4-2 19:10
我理解的退火过程就是逐步的降温升温

直接高温，再稳定至所需温度来查看体系的热稳定性，是不是拿0K时优化好的结构就能做啊？

作者
Author: 少年爱吃地三鲜 时间: 2019-4-3 08:16
你为啥要1200K那么高？

作者
Author: 灰天皮匠 时间: 2019-4-3 11:13

少年爱吃地三鲜发表于 2019-4-3 08:16
你为啥要1200K那么高？

想重复出一篇文献的工作，文献中是这么设定的

作者
Author: Andie 时间: 2020-4-6 10:22

灰天皮匠发表于 2019-4-3 11:13
想重复出一篇文献的工作，文献中是这么设定的

遇到了同样的困惑，审稿人质疑我用退火过程说明材料热稳定性的可靠性，麻烦楼主把要重复的文献挂出来，看看能不能提供一点理论支撑，谢谢啦

作者
Author: loading 时间: 2021-9-2 16:10
请问楼主，逐渐升温的命令应该怎么写呢？是fix 1 all nvt temp 300.0 2000.0 100.0么？怎么设定升温速率呢？谢谢您

作者
Author: xiaoj 时间: 2021-9-13 10:23

loading 发表于 2021-9-2 16:10
请问楼主，逐渐升温的命令应该怎么写呢？是fix 1 all nvt temp 300.0 2000.0 100.0么？怎么设定升温速率呢 ...

升温速率通过步数来控制，2000-300=1700 K 1700/步数=每步升温多少K（把步数换算成时间也可以）

作者
Author: loading 时间: 2022-10-5 09:27

xiaoj 发表于 2021-9-13 10:23
升温速率通过步数来控制，2000-300=1700 K 1700/步数=每步升温多少K（把步数换算成时间也可以）

懂了，谢谢楼主

欢迎光临计算化学公社 (http://bbs.keinsci.com/)