adf:gcmc2020
差别
两侧同时换到之前的修订记录前一修订版后一修订版 | 前一修订版上一修订版两侧同时换到之后的修订记录 |
adf:gcmc2020 [2022/01/20 20:37] – [参数设置] liu.jun | adf:gcmc2020 [2022/08/19 22:35] – [参数设置] liu.jun |
---|
======ReaxFF-GCMC:巨正则系综蒙特卡洛模拟====== | ======ReaxFF-GCMC:巨正则系综蒙特卡洛模拟====== |
应用案例: | =====巨正则系综===== |
| 组成系综的系统与一温度为T、化学势为μ的无限大粒子源相接触,此时系统不仅与热源存在能量交换,而且可以同粒子源有粒子交换,最后温度、离子数目处于动态平衡,这种系综称巨正则系综。 |
| |
| ====应用案例==== |
* [[https://www.fermitech.com.cn/ams/ams_application/reaxffhighlight201803|铂使用巨正则系综反应力场模拟铂电催化剂的氧化与失活(Small 2019, Angew. 2017)]] | * [[https://www.fermitech.com.cn/ams/ams_application/reaxffhighlight201803|铂使用巨正则系综反应力场模拟铂电催化剂的氧化与失活(Small 2019, Angew. 2017)]] |
| |
* μ<sub>O2</sub><sup>ref</sup>(T,P<sub>ref</sub>)是O<sub>2</sub>在温度T、压强P<sub>ref</sub>下,实验测得的化学势,如果查不到化学势,可以近似使用"[[adf:ir2020]]"中的吉布斯自由能作为化学势 | * μ<sub>O2</sub><sup>ref</sup>(T,P<sub>ref</sub>)是O<sub>2</sub>在温度T、压强P<sub>ref</sub>下,实验测得的化学势,如果查不到化学势,可以近似使用"[[adf:ir2020]]"中的吉布斯自由能作为化学势 |
* k为波尔兹曼常数,T为温度 | * k为波尔兹曼常数,T为温度 |
* P是O<sub>2</sub>的分压。如果只有O<sub>2</sub>一种气体,那么kT*ln(P/P<sub>ref<sub>)=kT*ln1=0 | * P是O<sub>2</sub>的分压。如果只有O<sub>2</sub>一种气体,那么kT*ln(P/P<sub>ref</sub>)=kT*ln1=0 |
* E<sub>O2</sub><sup>diss</sup>是分子的解离能,也就是分子解离为原子所需的能量(正值),可以使用ADF结构优化得到的Bond Energy的绝对值 | * E<sub>O2</sub><sup>diss</sup>可以使用ADF结构优化得到的Bonding Energy(<0,用于抵消Gibbs自由能中内能的部分。实际上更严格而言包括电子贡献的Bonding energy以及原子核贡献的内能两个部分,只是后者一般较小,两项在Gibss自由能的计算结果中,均可看到) |
| |
**分子的MC化学势:**一般分子的MC化学势计算公式,可以参考上述公式,只是不用乘以系数1/2 | **分子的MC化学势:**一般分子的MC化学势计算公式,可以参考上述公式,只是不用乘以系数1/2 |
adf/gcmc2020.txt · 最后更改: 2024/02/05 13:06 由 liu.jun