求助这种材料上的反应情况如何使用Gaussian计算反应的ΔG

2429783189

我想使用Gaussian算一个反应的ΔG，我研究的反应大致如下：A吸附在材料M上（为简单表示，后续用M-A表示吸附在M上的A），B加进来与A反应，A变成α（仍然吸附在M上），B变成β，此时，我要计算该反应的ΔG，我计算的公式应该是第一种：μ（M-α-β）-μ(M-A-B），即计算过程中只需要算复合物Mαβ（三者放在一个gjf文件里）和复合物MAB（三者放在一个gjf文件里）的吉布斯自由能并相减，还是第二种：μ（M-α）+μ（β）-μ（M-A）-μ（B），即算Mα复合物的吉布斯自由能+单独开一个gjf文件算β的吉布斯自由能-MA复合物的吉布斯自由能-单独开一个gjf文件算B的吉布斯自由能，

KazusaT

周期性结构表面吸附反应用gaussian不太好处理吧，最好考虑用cp2k或vasp等
不过也许可以用簇模型，可以看看sob老师写的簇模型的帖子

2429783189

KazusaT 发表于 2025-11-19 22:06
周期性结构表面吸附反应用gaussian不太好处理吧，最好考虑用cp2k或vasp等
不过也许可以用簇模型，可以看看 ...

对我就是参考的sob老师的帖子用的簇模型，而且收敛很顺利，只是在算ΔG的步骤上犯了难

KazusaT

2429783189 发表于 2025-11-19 22:11
对我就是参考的sob老师的帖子用的簇模型，而且收敛很顺利，只是在算ΔG的步骤上犯了难

我的理解是这样，如果考虑反应的过程是：
1. M + A + B （自由状态）
2. M-A + B （A吸附）
3. M-A + M-B （B吸附）
4. M-alpha + M-beta （反应）
5. M-alpha + beta （beta脱附，也可能alpha先脱附）
6. M + alpha + beta （alpha脱附）
那么第一种计算的就是反应一步的ΔG，第二种计算的包括了B的吸附、反应和beta脱附的ΔG

睡觉不打鼾

不加零点校正能，算出来的应该是电子能量，而不是Gibbs自由能吧

2429783189

KazusaT 发表于 2025-11-19 22:38
我的理解是这样，如果考虑反应的过程是：
1. M + A + B （自由状态）
2. M-A + B （A吸附）

感谢回复。

2429783189

睡觉不打鼾 发表于 2025-11-20 11:05
不加零点校正能，算出来的应该是电子能量，而不是Gibbs自由能吧

对我的计算加了零点校正能，忘记提了