adf:gcmc2020
差别
这里会显示出您选择的修订版和当前版本之间的差别。
| 两侧同时换到之前的修订记录前一修订版后一修订版 | 前一修订版 | ||
| adf:gcmc2020 [2020/11/28 16:35] – [ReaxFF-GCMC：巨正则系综蒙特卡洛模拟] liu.jun | adf:gcmc2020 [2024/02/05 13:06] (当前版本) – [结果查看] liu.jun | ||
|---|---|---|---|
| 行 1: | 行 1: | ||
| ======ReaxFF-GCMC：巨正则系综蒙特卡洛模拟====== | ======ReaxFF-GCMC：巨正则系综蒙特卡洛模拟====== | ||
| - | 应用案例： | + | =====巨正则系综===== |
| + | 组成系综的系统与一温度为T、化学势为μ的无限大粒子源相接触，此时系统不仅与热源存在能量交换，而且可以同粒子源有粒子交换，最后温度、离子数目处于动态平衡，这种系综称巨正则系综。 | ||
| + | |||
| + | ====应用案例==== | ||
| * [[https:// | * [[https:// | ||
| 本教程以Si表面的氧化为例。 | 本教程以Si表面的氧化为例。 | ||
| - | GCMC是AMS驱动的，各个模块使用起来，除了各个方法本身的参数（也就是Main窗口的设置）之外，其他参数设置是一样的。因此其他模块的GCMC模拟，也可以参考ReaxFF的GCMC模拟教程。 | + | GCMC是AMS驱动的，各个模块使用起来，除了各个方法（分子体系DFT、周期边界条件DFT、MOPAC、DFTB等）本身的参数（也就是Main窗口的设置）之外，其他参数设置是一样的。因此其他模块的GCMC模拟，也可以参考ReaxFF的GCMC模拟教程。 |
| =====建模===== | =====建模===== | ||
| 行 42: | 行 45: | ||
| ===MC化学势的计算=== | ===MC化学势的计算=== | ||
| - | **原子的MC化学势：**如果粒子源是原子，那么原子的MC化学势可以用统计力学从第一原理推导出来，或者从文献资料中的热化学表中查询。单个氧原子的化学势等于同温度压强下$O_2$的MC化学势的1/ | + | **原子的MC化学势：**如果粒子源是原子，那么原子的MC化学势可以用统计力学从第一原理推导出来，或者从文献资料中的热化学表中查询。单个氧原子的化学势等于同温度压强下O< |
| O原子的MC化学势计算公式： | O原子的MC化学势计算公式： | ||
| 行 48: | 行 51: | ||
| {{ : | {{ : | ||
| - | * 因为$O_2$是双原子分子，所以系数为1/ | + | * 因为O< |
| - | * μ$_{O2}^{ref}$(T,P$_{ref}$)是$O_2$在温度T、压强P$_{ref}$下，实验测得的化学势，如果查不到化学势，可以近似使用" | + | * μ<sub>O2</ |
| * k为波尔兹曼常数，T为温度 | * k为波尔兹曼常数，T为温度 | ||
| - | * P是$O_2$的分压。如果只有$O_2$一种气体，那么kT*ln(P/ | + | * P是O< |
| - | * E$_{O2}^{diss}$是分子的解离能，也就是分子解离为原子所需的能量（正值），可以使用ADF结构优化得到的Bond Energy的绝对值 | + | * E<sub>O2</ |
| **分子的MC化学势：**一般分子的MC化学势计算公式，可以参考上述公式，只是不用乘以系数1/ | **分子的MC化学势：**一般分子的MC化学势计算公式，可以参考上述公式，只是不用乘以系数1/ | ||
| 行 58: | 行 61: | ||
| {{ : | {{ : | ||
| =====结果查看===== | =====结果查看===== | ||
| - | 在movie中可以看到Si的氧化结果： | + | 在movie中可以看到Si的氧化中间过程： |
| {{ : | {{ : | ||
| + | |||
| + | 能量曲线平稳后，系统也达到平衡状态。 | ||
| + | ====键级==== | ||
| + | 在AMSinput中的键级，未经计算，因此其键级是图形窗口根据简单的价电子规则而猜测出来的。在Movie中是经过计算的，因此鼠标选中两个原子，窗口左下角则会显示两个原子之间的键级。 | ||
adf/gcmc2020.1606552556.txt.gz · 最后更改: 2020/11/28 16:35 由 liu.jun