表面能计算公式问题

追枫少年

卡开发发 发表于 2018-1-12 11:50
0.1eV的振荡已经不是很大了，你可以考虑再加上workfunction修正看看会不会好点。

感谢老师的回复
因为我看的文献震荡小于10-3的量级（DOI: 10.1039/c4cp06089j），为我做的和他相比差的过大，所以一直认为那里除了问题，我在尝试计算功函数分析一下结果，再次感谢老师!

zorow

卡开发发 发表于 2018-1-12 10:56
一般情况相同原子slab比bulk能量高，slab之间距离拉开后能量升高还是正常的。

谢谢卡开发发老师，那是不是真空层厚度就不用做这样的收敛性测试了，取一个12以上的厚度即可？

卡开发发

zorow 发表于 2018-1-12 15:00
谢谢卡开发发老师，那是不是真空层厚度就不用做这样的收敛性测试了，取一个12以上的厚度即可？

原则上要做，不过经验上一般15A左右应该足够。

zorow

卡开发发 发表于 2018-1-12 15:58
原则上要做，不过经验上一般15A左右应该足够。

感谢您的回复，我取了11、13、15、17、19的厚度，但是计算出的表面能分别为41.30、41.32、41.34、41.36、41.37(J/m2)，感觉似乎一直在均匀地增加，所以对这个收敛性测试比较困惑。。。
那我就取15埃的厚度好了

卡开发发

zorow 发表于 2018-1-12 16:08
感谢您的回复，我取了11、13、15、17、19的厚度，但是计算出的表面能分别为41.30、41.32、41.34、41.36、 ...

结构的上下表面对称否，不对称的话可能得把功函数修正加上。

zorow

卡开发发 发表于 2018-1-12 16:29
结构的上下表面对称否，不对称的话可能得把功函数修正加上。

结构的上下表面貌似是一样的，但因为不平整所以都是向上凸的，这样的需要功函数修正吗

卡开发发

zorow 发表于 2018-1-12 16:50
结构的上下表面貌似是一样的，但因为不平整所以都是向上凸的，这样的需要功函数修正吗

不需要。按说没有上下表面不对称的情况，15A左右应该差不多了，如果很有必要的话，你可以测试更大的真空层看下和现在有啥区别，测试1x1的表面做单点就行。

zorow

卡开发发 发表于 2018-1-12 19:12
不需要。按说没有上下表面不对称的情况，15A左右应该差不多了，如果很有必要的话，你可以测试更大的真空 ...

好的卡开发发老师，谢谢您的指点！

小潘潘

同学，你好，我是一名量化计算小白，请问你有没有如何用MS计算表面能的文章？最好是步骤详尽些的，以便于我可以模仿，谢谢啦

granvia

zorow 发表于 2018-1-12 10:19
卡开发发老师，我还有一个问题请教您：
在测试真空层厚度时，居然得出了厚度越大表面能越大的结果，请问 ...

可以算算不同厚度下的功函数，看看厚度是否收敛。这个vasp的handout里有算例

YMZoey

卡开发发 发表于 2018-1-12 11:07
两个面是否是一样的没关系，比如一个AB两种原子构成的晶体沿着晶面方向原子的排列有可能是...ABABAB...， ...

请问，如果上下两个表面不同的话，那算出来的表面能是哪一个表面的能量呢？

卡开发发

YMZoey 发表于 2019-6-3 09:03
请问，如果上下两个表面不同的话，那算出来的表面能是哪一个表面的能量呢？

符合计量比的情况应该是这个晶体劈裂形成这两种表面的平均能量。

ZZ爱学习呀

zorow 发表于 2018-1-11 20:43
谢谢卡开发发老师的指点！
假如我的bulk有4个A原子，4个B原子，而slab则为8A8B，则n的值即为8/4=2，是这 ...

您好，请问计算表面能的时候，切面的层数以及是否需要扩胞应该怎么选取呀？

卡开发发

ZZ爱学习呀 发表于 2020-4-3 04:46
您好，请问计算表面能的时候，切面的层数以及是否需要扩胞应该怎么选取呀？

没掺杂的情况没必要扩胞，等同于使用平移对称性。层数的确定需要做层数和表面能的收敛性测试。

ZZ爱学习呀

卡开发发 发表于 2018-1-12 11:07
两个面是否是一样的没关系，比如一个AB两种原子构成的晶体沿着晶面方向原子的排列有可能是...ABABAB...， ...

老师，请问计算表面能的时候，切面的层数以及是否需要扩胞应该怎么选取呀？