adf:gcmc

差别

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

到此差别页面的链接

--- adf:gcmc [2019/06/28 07:50] – liu.jun
+++ adf:gcmc [2020/05/18 16:12] (当前版本) – [注意1：] liu.jun
@@ 行 1: / 行 1: @@
-======如何进行巨正则系综蒙特卡洛GCMC模拟======
+======2016以前的版本如何进行巨正则系综蒙特卡洛GCMC模拟======
 体系与一温度为T的粒子（分子、原子）源接触，体系不仅同粒子源有能量交换，而且可以同粒子源有粒子的交换，最后达到平衡，这种系综称巨正则系综，可以用于模拟块体表面吸附小分子、原子的反应过程。
-=====新版本ADF=====
-ADF2017.211之后的版本，将GCMC的功能集成到了图形界面，具体可以参考教程[[adf:voltageprofileofdischarging]]中的第3步“对Alpha硫晶体进行巨正则系综蒙特卡洛模拟”
-=====旧版本ADF=====
 以铁表面的吸附O<sub>2</sub>、H<sub>2</sub>O为例来说明。
-<WRAP center round download 30%>
-[[https://www.jianguoyun.com/p/DUpDeS8QmZ2ZBhiR-CQ|计算文件下载（点击）]]
-</WRAP>
 ====建模：====
@@ 行 149: / 行 144: @@
 ===注意2：化学势的计算===
 关于‘cmpot   !chemical potential of molecule’这一项，上面的例子中，化学势都设置为-75kcal/mol。化学势计算公式：
-[μ(T,P)+]+k<sub>b</sub>*ln(P/P')-E<sub>d</sub>]/2
+[μ(T,P)+k<sub>b</sub>T*ln(P/P')-E<sub>d</sub>]/2
-  * μ(T,P)可以通过计算该分子的气相Gibbs自由能得到，参考[[adf:gibbsofgas]]。计算的时候，所取温度、压强与系综温度、压强一致即可。
+  * μ(T,P)是在温度、压强下，实验测得的化学势
-  * k<sub>b</sub>为波尔兹曼常熟。
+  * k<sub>b</sub>为波尔兹曼常数，T为温度
-  * P'是小分子的分压强。如果只有一种气体，那么k<sub>b</sub>*ln(P/P')＝k<sub>b</sub>*ln1＝0
+  * P'是小分子的分压强。如果只有一种气体，那么k<sub>b</sub>T*ln(P/P')＝k<sub>b</sub>T*ln1＝0
-  * E<sub>d</sub>是分子的解离能，也就是ADF计算得到的分子的Bonding Energy。
+  * E<sub>d</sub>是分子的解离能
 ====提交任务====

adf/gcmc.1561679446.txt.gz · 最后更改: 2019/06/28 07:50 由 liu.jun

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