这里会显示出您选择的修订版和当前版本之间的差别。
两侧同时换到之前的修订记录前一修订版后一修订版 | 前一修订版 | ||
adf:gcmc2020 [2021/02/24 10:38] – [ReaxFF-GCMC:巨正则系综蒙特卡洛模拟] liu.jun | adf:gcmc2020 [2024/02/05 13:06] (当前版本) – [结果查看] liu.jun | ||
---|---|---|---|
行 1: | 行 1: | ||
======ReaxFF-GCMC:巨正则系综蒙特卡洛模拟====== | ======ReaxFF-GCMC:巨正则系综蒙特卡洛模拟====== | ||
- | 应用案例: | + | =====巨正则系综===== |
+ | 组成系综的系统与一温度为T、化学势为μ的无限大粒子源相接触,此时系统不仅与热源存在能量交换,而且可以同粒子源有粒子交换,最后温度、离子数目处于动态平衡,这种系综称巨正则系综。 | ||
+ | |||
+ | ====应用案例==== | ||
* [[https:// | * [[https:// | ||
行 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中是经过计算的,因此鼠标选中两个原子,窗口左下角则会显示两个原子之间的键级。 |