求助CP2K中Gibbs自由能的计算

桂の猿 · 发表于 Post on 2024-5-5 22:18:10

各位老师好

我最近在计算完过渡态和中间体后，希望计算不同温度下的反应物，过渡态，中间体和产物的吉布斯自由能，从而根据TST理论推算不同温度下的反应速率常数。
我使用的软件是CP2K.在振动频率分析后，利用Shermo软件处理out文件，得到结果。

对于不同温度下的吉布斯自由能的计算，是否只需要将振动频率计算中的TC_TEMPERATURE 298.15 修改就可以？

&VIBRATIONAL_ANALYSIS
  DX 0.001
  NPROC_REP 1
  TC_PRESSURE 100000
  TC_TEMPERATURE 298.15
  THERMOCHEMISTRY
  INTENSITIES T
  FULLY_PERIODIC T
&END VIBRATIONAL_ANALYSIS

谢谢

桂の猿 · 发表于 Post on 2024-5-6 00:31:15

sobereva 发表于 2024-5-5 23:33
完全不需要改TC_TEMPERATURE，根本不影响输出的振动频率、Shermo的计算结果。TC_TEMPERATURE仅仅影响CP2K直 ...

谢谢卢老师，我以前算的时候，发现273和1473的结果一样，所以，一直以为是自己的问题，非常感谢。

sobereva · 发表于 Post on 2024-5-5 23:33:11

完全不需要改TC_TEMPERATURE，根本不影响输出的振动频率、Shermo的计算结果。TC_TEMPERATURE仅仅影响CP2K直接给出的热力学量
算什么温度下的热力学量，就在Shermo的settings.ini里把T改成什么就完了

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[CP2K] 求助CP2K中Gibbs自由能的计算