今日添加了一个g-C3N4与CeO2 (111)异质结的建模例子。

顺便再放几点相关讨论：
1. 表格里Results旁边有一块区域，原本是希望直接输出变换矩阵的，但在我这里返回HexGrids查找相应单元格位置的Excel函数有点故障，只能留空了。实际上对HexGrid用ctrl+f查找也不完全可靠，哪怕Excel里公式 - 计算 - 计算选项设置为自动来生成整个HexGrid，查找时也有可能返回0个结果。
2. 表格里用的失配率是按两种超胞的晶胞参数之差占晶胞参数的比例来计算的，设置一个如1%的失配率上限对于小晶胞还好，而对于比较大的晶胞则必须留意实际空间距离误差超过0.5埃（随便举的例子，不是严格的阈值）导致原子有异常接触的问题，建议降低失配率上限、增大超胞尺寸上限来找合适的解。
3. 所谓“根号扩胞”其实是英文名 Wood's notation 的记号的不太严谨的简化，而晶格矢量变换矩阵对应的记号就叫Matrix notation，这在讨论表面吸附的维基Overlayer页有简单介绍，用Google搜 Wood's notation 也能找到一些课件。
4. 对twistronics研究的双层二维转角体系如平行错位pi-pi堆积的转角石墨烯，建孤立团簇模型用量子化学软件做计算是比较常见的做法，我估计这可能与建周期性模型不好选取适应实际夹角对应的点阵矢量有关，比如$64 \vec{a} + \vec{b}$与$64 \vec{a} - \vec{b}$的夹角也有1.5627度，$63 \vec{a} + 31\vec{b}$与$63 \vec{a} + 32\vec{b}$的夹角为1.0501度，然而此时超胞边长已经是原胞的五六十倍了，原子数更是接近三四千倍，计算量极大。

好像有个笔误 “旁边的28对应矢量”应该是“旁边的63对应矢量”对不对