求助CP2K中Gibbs自由能的计算

桂の猿

各位老师好


我最近在计算完过渡态和中间体后，希望计算不同温度下的反应物，过渡态，中间体和产物的吉布斯自由能，从而根据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



谢谢

sobereva

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

桂の猿

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

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