用户工具

站点工具


adf:gcmc2020

差别

这里会显示出您选择的修订版和当前版本之间的差别。

到此差别页面的链接

两侧同时换到之前的修订记录前一修订版
后一修订版
前一修订版
上一修订版两侧同时换到之后的修订记录
adf:gcmc2020 [2022/01/20 20:37] – [参数设置] liu.junadf:gcmc2020 [2022/08/19 22:35] – [参数设置] liu.jun
行 1: 行 1:
 ======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)]]
  
行 51: 行 54:
   * μ<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

© 2014-2022 费米科技(京ICP备14023855号