求助，超晶胞建模的一些问题

慕容轩逸

老师们好，我现在建立超晶胞遇到一个问题，我有两种做法不知道对不对，一个是先建立2*2*1的晶胞，在进行切面，一个是先切面，在建立2*2的超晶胞，（真空层是加在Z方向的）。不知道这两种做法是否都可以？
还有就是在切面时的thickness和真空层厚度都是需要做收敛性测试的吗？

starrywang

1. 无论如何，晶胞构建后需要优化；
2. 2*2*1已经算是超晶胞了吧；
3. 建立2*2*1超晶胞后，如果不优化的话，那么建立超晶胞与切面的顺序应该不影响；
4. 切面厚度和真空层高度是需要参考文献以及自己测试的。

卡开发发

建议先切面，优化，再建立表面的超晶格。这样做的原因是，参数选择可靠的情况下，切面1x1的表面优化得到的结果建立超晶格和超晶格做优化原则上应该相同，但直接计算1x1的表面计算量却小很多，这等同于使用了平移对称性。如果你觉得不放心，可以在1x1优化完建立超晶格再进行一次优化，原则上1步就能完成。楼上提到的厚度和层高也是在1x1表面进行效率更高一些。

慕容轩逸

卡开发发 发表于 2020-1-7 20:09
建议先切面，优化，再建立表面的超晶格。这样做的原因是，参数选择可靠的情况下，切面1x1的表面优化得到的 ...

您好，谢谢卡开发发老师耐心的解答，我想在请教您几个问题，
1.如果我是按照先切面在建立超晶胞的顺讯的话，我切面之前选择的是晶体库里的结构，然而它并不是原胞，我是不是可以转化成原胞在切面呢？
2.我建立好模型之后，是要添加吸附分子的，我是在吸附之前优化一次，在添加吸附分子在优化一次；还是直接建立吸附模型在优化呢？
3.在优化的时候我看小*虫上之前有的帖子并没有选择优化晶胞参数，因为他选择切面的模型就是优化过晶胞参数的了，不知这样可否？
4，我做的是表面吸附，不知道是否需要固定除去表面几层原子之外的原子？
希望老师您不吝赐教，感谢感谢

慕容轩逸

starrywang 发表于 2020-1-7 10:42
1. 无论如何，晶胞构建后需要优化；
2. 2*2*1已经算是超晶胞了吧；
3. 建立2*2*1超晶胞后，如果不优化的 ...

老师，我想在问一下您几个问题，因为我这个模型目前我还并没有找到做过吸附模型的相关文献，所以恐怕需要我自己做收敛性测试了。
1.需要做测试的参数就是这两张图片上的Thickness（切面厚度）和Vacuum orientation(真空层厚度）两个吧，测试标准应该也是和普通的收敛性测试观察最终能量的变化就可以吗？
2.我在网上看有的人说还要测试原子层数，不知道这个原子层数指的是哪个参数呢？
希望老师您有时间了可以帮忙解释一下。

卡开发发

慕容轩逸 发表于 2020-1-8 10:27
您好，谢谢卡开发发老师耐心的解答，我想在请教您几个问题，
1.如果我是按照先切面在建立超晶胞的顺讯的 ...

1、晶面一般都是按照单胞定义的，优化可以用原胞。
2、建议：1x1的表面优化一下，一般很快就能完成，然后扩展超晶格，加上吸附分子，再进行优化。
3、计算表面也没法优化晶格参数，否则两块平板之间有相互作用，可能真空层塌掉也可能根本收敛不了。
4、我认为没有必要，原因如下：
（1）固定是非物理手段。即便固定下表面，下表面仍然和真空接触，电子态也会受到这个影响，所谓保持和块状固体一致的性质是不成立的。如果要模拟一个实际的固体表面，应该给足够的层高，这样中间的区域才会有类似于块状固体一样的性质。
（2）固定并不保证能够获得加速。是否做固定并不会减小KS方程的求解的计算量，不合适的固定甚至可能导致结构优化每一步SCF的迭代数目被拉长甚至完全收敛不了。

starrywang

慕容轩逸 发表于 2020-1-8 10:35
老师，我想在问一下您几个问题，因为我这个模型目前我还并没有找到做过吸附模型的相关文献，所以恐怕需要 ...

1. 真空层测试和你的吸附分子有关，有的分子很大，而且吸附姿态是“站立”的，那么层高低了的话，实际分子会受上下两个表面的影响。我个人喜欢通过吸附结构和吸附能来判断层高怎么样。这点还需要卡开发发老师来把把关。
2. 原子层数就是thickness。

慕容轩逸

卡开发发 发表于 2020-1-8 12:38
1、晶面一般都是按照单胞定义的，优化可以用原胞。
2、建议：1x1的表面优化一下，一般很快就能完成，然 ...

非常感谢卡开发发老师耐心的解答。

慕容轩逸

starrywang 发表于 2020-1-8 13:38
1. 真空层测试和你的吸附分子有关，有的分子很大，而且吸附姿态是“站立”的，那么层高低了的话，实际分 ...

非常感谢！您的回复对我很有帮助

慕容轩逸

卡开发发 发表于 2020-1-8 12:38
1、晶面一般都是按照单胞定义的，优化可以用原胞。
2、建议：1x1的表面优化一下，一般很快就能完成，然 ...

卡开发发老师，还想请问一下您，关于真空层和切面厚度的测试是以什么为判据的呢？

卡开发发

慕容轩逸 发表于 2020-1-8 19:30
卡开发发老师，还想请问一下您，关于真空层和切面厚度的测试是以什么为判据的呢？

表面能或者功函数可能比较合适。

dhdhdh

卡开发发 发表于 2020-1-8 22:42
表面能或者功函数可能比较合适。

老师，我需要构建氢终止的金刚石和水分子的界面。那我是不是需要在金刚石原胞上进行优化，切面再构建超晶胞。那是要再加赝氢之后再对超晶胞进行优化吗，还是先优化再加氢呢？最后加上水分子之后再进行固定优化这样对吗

卡开发发

dhdhdh 发表于 2021-6-30 20:49
老师，我需要构建氢终止的金刚石和水分子的界面。那我是不是需要在金刚石原胞上进行优化，切面再构建超晶 ...

切面后不急着建超晶格，然后直接加H再优化一下，这样应该更快一些，最后再放水。

dhdhdh

卡开发发 发表于 2021-7-1 00:14
切面后不急着建超晶格，然后直接加H再优化一下，这样应该更快一些，最后再放水。

好的，谢谢老师~

dhdhdh

卡开发发 发表于 2021-7-1 00:14
切面后不急着建超晶格，然后直接加H再优化一下，这样应该更快一些，最后再放水。

卡开发发老师，我还有点疑问。
1.在通过计算表面能对原子层厚度进行测试的时候是否需要固定中间的原子层呢，我看有的说要固定中间原子层，弛豫上下底面，有的说要固定下层。
2.测试完之后加赝氢的话，是必须要严格按照化学键中点的位置加呢，还是可以直接加上氢原子再优化呢