计算化学公社

标题: 求助Gaussian的ADMP的热浴如何设定？ [打印本页]

作者
Author: 一昔a 时间: 2024-2-24 12:53
标题: 求助Gaussian的ADMP的热浴如何设定？
在阅读社长的谈谈温度、压力、同位素设定对量子化学计算结果产生的影响 - 思想家公社的门口：量子化学·分子模拟·二次元 (sobereva.com)文章中，谈及温度对从头算动力学的影响中提到：
“写temperature关键词没有用处，得用IOp来专门地设定热浴。”
想请教下各位热浴如何设定呢，是采用IOp(1/111=N) 来设定吗，为什么采用temperature不能对从头算动力学进行温度设定呢？
我在进行某个反应的从头算动力学的过程中，发现总是跨不过能垒，是不是与温度设定有关，在温度保持不变的同时与NKE有关系吗，NKE是不是设的大一点好？

作者
Author: 一昔a 时间: 2024-2-24 13:21
查找了几个相关帖子，贴在这里了，大家如果有需要的话可以查看
http://sobereva.com/423
http://bbs.keinsci.com/thread-35994-1-1.html
Gaussian IOp手册 https://gaussian.com/iops/
IOp手册原文
IOp(1/80)

L106: 0/1 Cartesian/Normal mode/ Internal coordinate differentiation. 2 is NYI.

L118: = 1 to suppress the 5th order correction after surface hop has been made in Trajectory Surface Hopping calculations. Needs also IOp(10/80 = 1)

L121: Nuclear Kinetic Energy Thermostat Option (currently only velocity scaling is implemented).

0 No Thermostat.
11XXXXX Velocity scaling, but only for the first XXXXX simulation steps. (This option is useful, if thermostating is only required during equilibration.)
1000000 Velocity scaling, all the way through the simulation.

IOp(1/81)

Nuclear KE thermostat in ADMP — temperature is checked and scaled every IOp(81) steps.

IOp(1/82)

L121: Temperature for nuclear KE thermostat.

作者
Author: sobereva 时间: 2024-2-25 07:49
北京科音科音量化高级班专门有讲

(, 下载次数 Times of downloads: 16)

作者
Author: 一昔a 时间: 2024-2-26 16:04

sobereva 发表于 2024-2-25 07:49
北京科音科音量化高级班专门有讲

太谢谢老师了！我就是想了解下这个

作者
Author: 一昔a 时间: 2024-2-26 16:07

sobereva 发表于 2024-2-25 07:49
北京科音科音量化高级班专门有讲

老师，不知道您做过这方面的测试吗，加不加fullscf对计算速度有多大影响呢，我算的一个七八十个原子的体系，加上fullscf后算了2个小时才算了70多步，是不是不加fullscf好些呢。

作者
Author: wzkchem5 时间: 2024-2-26 17:19

一昔a 发表于 2024-2-26 09:07
老师，不知道您做过这方面的测试吗，加不加fullscf对计算速度有多大影响呢，我算的一个七八十个原子的体 ...

最简单的办法是自己不加fullscf跑一下试试，甚至不需要跑2个小时，跑15分钟、然后和fullscf的前15分钟的结果对比就行了，远比发帖等回复快、靠谱。
体系的scf收敛困难程度不同时，fullscf对时间的影响也不同，如果不知道你体系的细节，没有办法精确估计fullscf具体对计算速度影响多大

作者
Author: 一昔a 时间: 2024-2-26 18:11

wzkchem5 发表于 2024-2-26 17:19
最简单的办法是自己不加fullscf跑一下试试，甚至不需要跑2个小时，跑15分钟、然后和fullscf的前15分钟的 ...

好的，谢谢老师！

作者
Author: sobereva 时间: 2024-2-27 01:22

一昔a 发表于 2024-2-26 16:07
老师，不知道您做过这方面的测试吗，加不加fullscf对计算速度有多大影响呢，我算的一个七八十个原子的体 ...

Gaussian里ADMP能比BOMD快关键就在于不用每一步令波函数迭代到收敛。你写了fullscf，ADMP就没任何优势了

作者
Author: 一昔a 时间: 2024-2-27 11:06

sobereva 发表于 2024-2-27 01:22
Gaussian里ADMP能比BOMD快关键就在于不用每一步令波函数迭代到收敛。你写了fullscf，ADMP就没任何优势了

嗯嗯收到，谢谢老师！

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