这里会显示出您选择的修订版和当前版本之间的差别。
两侧同时换到之前的修订记录前一修订版后一修订版 | 前一修订版 | ||
atk:用准谐近似理论研究反应速率 [2016/12/08 15:50] – [优化 NEB 与 HTST 速率计算] dong.dong | atk:用准谐近似理论研究反应速率 [2016/12/08 15:54] (当前版本) – [参考文献] dong.dong | ||
---|---|---|---|
行 104: | 行 104: | ||
==== 结果分析 ==== | ==== 结果分析 ==== | ||
- | 计算完成后,文件'' | + | 计算完成后,文件 '' |
{{ : | {{ : | ||
- | 势能曲线表明能垒为1.3eV。反应的动画展示了预期的跳跃过程。在初始态和末态之间有一个单能垒。关闭Movie Tool窗口。选择LabFloor中的HTSTEvent对象并点击右侧面板的HTST Rates插件。 | + | 势能曲线表明能垒为 1.3 eV。反应的动画展示了预期的跳跃过程。在初始态和末态之间有一个单能垒。关闭 Movie Tool 窗口。选择 LabFloor 中的 HTSTEvent 对象并点击右侧面板的 HTST Rates 插件。 |
{{ : | {{ : | ||
- | 显示了关于反应的信息如能垒和指前因子。可以通过改变温度值来计算任意温度下的速率,默认温度为300K。 | + | 显示了关于反应的信息如能垒和指前因子。可以通过改变温度值来计算任意温度下的速率,默认温度为 300K。 |
{{ : | {{ : | ||
- | 正反应和逆反应的能垒和前因子在理论上应该相同,然而在实际过程中是不同的。这是由于NEB端点几何构型的微小不同(它们每一个都是从一个略微不同的开始结构被最小化的)。而且,由于这是一个ATK-DFT计算,由于eggbox效应,平移对称性没有被完美的保持。这个效应描述了用以表示电荷密度和势的实空间网格与原子位置之间的关系。这些网格点会随着不同原子而向不同方向对齐,从而打破了平移对称性。在本实例中对于总能影响非常小(< | + | 正反应和逆反应的能垒和前因子在理论上应该相同,然而在实际过程中是不同的。这是由于 NEB 端点几何构型的微小不同(它们每一个都是从一个略微不同的开始结构被最小化的)。而且,由于这是一个 ATK-DFT 计算,由于 eggbox 效应,平移对称性没有被完美的保持。这个效应描述了用以表示电荷密度和势的实空间网格与原子位置之间的关系。这些网格点会随着不同原子而向不同方向对齐,从而打破了平移对称性。在本实例中对于总能影响非常小(< |
- | 我们也可以画一个Arrhenius图来显示HTST速率随着温度的变化情况。一个Arrhenius图是$\log({\rm rate})$与$1/ | + | 我们也可以画一个 Arrhenius 图来显示 HTST 速率随着温度的变化情况。一个 Arrhenius 图是 $\log({\rm rate})$ 与 $1/T$ 的关系。当与温度的关系以这种方式绘出,图应该总是直线。直线的斜率是 $-\Delta E/k_{\rm B}$ 而 y 轴截距是指前因子。从图可以很容易地说明指前因子是 HTST 预测的最大可能速率。 |
{{ : | {{ : | ||
- | Pt的跳跃扩散反应速率在室温下非常低。反应之间的平均时间(速率的倒数)是在世纪量级的。然而,在600K我们可以期望反应每秒发生几十次。 | + | Pt 的跳跃扩散反应速率在室温下非常低。反应之间的平均时间(速率的倒数)是在世纪量级的。然而,在 |
===== 计算多个基元反应步的速率 ===== | ===== 计算多个基元反应步的速率 ===== | ||
- | Pt吸附原子在Pt(100)面扩散的交换机制是一个很有趣的学习案例,因为它不像跳跃机制那样通过一个单元反应来进行。计算的设置细节在另一个[[http:// | + | Pt 吸附原子在 Pt(100)面扩散的交换机制是一个很有趣的学习案例,因为它不像跳跃机制那样通过一个单元反应来进行。计算的设置细节在另一个[[http:// |
- | 收敛的CI-NEB计算显示了如下的反应坐标图表: | + | 收敛的 CI-NEB 计算显示了如下的反应坐标图表: |
{{ : | {{ : | ||
- | 这里你可以看到有两个鞍点因而有两个反应步。第一个能垒大约为1.4eV,第二个能垒约为0.35eV。对于这个反应的HTST速率计算将会涉及将这个计算分成两个分开的NEB构型以确定每一步反应的速率常数。 | + | 这里你可以看到有两个鞍点因而有两个反应步。第一个能垒大约为 1.4 eV,第二个能垒约为 0.35 eV。对于这个反应的 HTST 速率计算将会涉及将这个计算分成两个分开的 NEB 构型以确定每一步反应的速率常数。 |
- | 根据你试图想要构建的模型,忽略中间体态通常会是一个合理的近似。举例来说,在600K初始反应步的速率大约为1.7Hz,而第二步的反应速率大约为1GHz($10^{12}\exp[-0.35 eV/ 600 K k_{\rm B}]$)。所以,忽略反应中间体可能是合理的。 | + | 根据你试图想要构建的模型,忽略中间体态通常会是一个合理的近似。举例来说,在600 K 初始反应步的速率大约为 1.7 Hz,而第二步的反应速率大约为 1GHz($10^{12}\exp[-0.35 eV/ 600 K k_{\rm B}]$)。所以,忽略反应中间体可能是合理的。 |
- | 更细致的研究应当运行一个动力学Monte Carlo模拟来确定系统的动力学。只要涉及的反应步超过一个,反应动力学就不会严格的服从一阶速率方程,尽管它通常是一个有效的近似。 | + | 更细致的研究应当运行一个动力学 Monte Carlo 模拟来确定系统的动力学。只要涉及的反应步超过一个,反应动力学就不会严格的服从一阶速率方程,尽管它通常是一个有效的近似。 |
===== 参考文献 ===== | ===== 参考文献 ===== | ||
- | [1] George H. Vineyard. Frequency factors and isotope effects in solid state rate processes. Journal of Physics and Chemistry of Solids, 3(1):121 – 127, 1957. doi: | + | * [1] George H. Vineyard. Frequency factors and isotope effects in solid state rate processes. Journal of Physics and Chemistry of Solids, 3(1):121 – 127, 1957. doi: |
- | + | | |
- | [2] Henry Eyring. The activated complex in chemical reactions. The Journal of Chemical Physics, 3(2): | + | |
- | + | * 原文:http:// | |
- | [3] Arthur F. Voter and Jimmie D. Doll. Transition state theory description of surface self‐diffusion: | + | |
- | + | ||
- | 本文翻译:王吉章 | + |