计算化学公社

标题: 求助CP2K中Gibbs自由能的计算 [打印本页]

作者
Author: 桂の猿 时间: 2024-5-5 22:18
标题: 求助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

谢谢

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

作者
Author: 桂の猿 时间: 2024-5-6 00:31

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

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

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