求助cp2k计算此反应OOH*+H→O*+H2O 吉布斯，这样算是否正确？

qzm

1式子： *+O2+H→OOH*     
2式子：OOH*+H→O*+H2O   
吉布斯自由能
步骤1：先计算O2分子能量 -31.881592786405143  （把优化好的一个O2分子放在一个盒子里，cp2k计算获得能量）
           先计算H2分子能量 -1.1613144308840  （把优化好的一个H2分子放在一个盒子里，cp2k计算获得能量）
          先计算H2O分子能量-17.219987932384019  （把优化好的H2O分子放在一个盒子里，cp2k计算获得能量）
          计算频率  求出他们的 G(O2)   G(H2)   G(H2O)
      不确定是否准确，
步骤2：计算OOH*的能量  *的能量，O*的能量   分别计算他们的频率 求出 G（OOH*） G（*） G（OH*）

计算  1式子    吉布斯自由能= G（OOH*） - G（*）- G(O2) - 1/2G(H2)

          2式子    吉布斯自由能= G（O*）+ G (H2O) - G（OOH*） - 1/2G(H2)

求助：1吉布斯自由能这样计算是否正确？
         2给出了O2 H2 H2O 计算的输入和输出文件，这样计算获得的结果是否可以写论文直接用

sobereva

没具体OOH*的输入文件没法说
O2基态是三重态
孤立体系还用k点莫名其妙
孤立体系用非周期性明显是更好选择
注意对于算能量目的，DZVP-MOLOPT-SR-GTH太烂，这只适合opt freq目的

qzm

sobereva 发表于 2024-5-1 22:23
没具体OOH*的输入文件没法说
O2基态是三重态
孤立体系还用k点莫名其妙

感谢社长，更新了OOH的输入，计算能量选TZ2P-MOLOPT 可以吗， 请问计算步骤是否存在错误呢？