表面计算

小白装大神

想做一些表面掺杂的工作，看了一些资料，还是有很多地方不太明白，希望各位可以给指点一下，请问：
1.原子层数、真空层厚度以及固定的原子层数先对哪个进行测试？
2.k点和截断能的测试是在以上三种测试之前还是之后进行？
3.对于某些晶面在满足化学计量比的情况下只能是极性面，然后我对底层原子加赝氢饱和后进行表面能的计算，这样可以吗？还是得进行偶极修正？
4.计算掺杂表面形成能的公式如下，是否合理？
E（formation）=Etot（掺杂表面）-Etot（未掺杂表面）+n(被替换原子）u(被替换原子）-n（掺杂原子）u（掺杂原子）[n，u分别指原子个数以及对应的化学势]
希望各位大神不吝赐教，谢谢！

卡开发发

2、先测试k和截断能；
1、可以先取一个很大的真空层（可以拿1x1的单胞来测试）把层数测试出来，然后再按照测试得到的层数测试真空；
3、底层用赝氢饱和还能够算表面能吗（从表面能定义出发）？应该可以考虑计算表面能的时候不进行赝氢饱吧。和如果不希望另外一面的表面态出现，通过赝氢饱和，严格来说，两个表面结构不对称都应该进行偶极校正，更重要的看关注什么样的性质，比如work function那就不得不进行校正；
4、如果替换的原子化合价和被替换的相同，应该没有问题，不过要注意怎么定义u.

表面参杂我觉得是很难做的一个工作，从表层到内层，每一层的原子都不是等价的，可能性诸多。

小白装大神

卡开发发 发表于 2015-8-29 17:27
2、先测试k和截断能；
1、可以先取一个很大的真空层（可以拿1x1的单胞来测试）把层数测试出来，然后再按照 ...

谢谢你的回答！您说先应该测试k点和截断能，但是所需要的原子层数和真空层厚度都是未知的，测试应该怎么做，我可以选取体材料优化所测得的截断能，然后根据表面模型的晶格取一个k点，在确定原子层数和真空层厚度之后，进行表面计算之前再测取合理的k点和截断能。
我关注的是表面掺杂形成能这一块，是否也必须偶极矫正，需要的话应该怎么做，CASTEP可以实现吗？在此表示感谢！

卡开发发

小白装大神 发表于 2015-8-31 09:27
谢谢你的回答！您说先应该测试k点和截断能，但是所需要的原子层数和真空层厚度都是未知的，测试应该怎么 ...

截断能在计算Bulk的体系就能够确定了，k可以按照与Bulk相同的倒空间尺寸来取，之后直接取个真空大一点（如30A）的1x1的slab，然后测试层数，之后只要反过来在确定的层数下测试真空大小是否充分就行。

原则上应该偶极校正，在SCF的那个选项卡能够直接找到。

小白装大神

卡开发发 发表于 2015-8-31 10:14
截断能在计算Bulk的体系就能够确定了，k可以按照与Bulk相同的倒空间尺寸来取，之后直接取个真空大一点（ ...

非常感谢

wernea

请问表面弛豫的时候固定底层原子，然后表面的原子是允许改变xyz，还是固定xy，只允许改变z。

小白装大神

wernea 发表于 2015-8-31 20:37
请问表面弛豫的时候固定底层原子，然后表面的原子是允许改变xyz，还是固定xy，只允许改变z。

允许改变xyz

wernea

小白装大神 发表于 2015-8-31 21:53
允许改变xyz

好的，谢谢

小白装大神

wernea 发表于 2015-8-31 22:17
好的，谢谢

有没有接触过极性表面，怎么处理的？

wernea

小白装大神 发表于 2015-9-1 21:28
有没有接触过极性表面，怎么处理的？

没有，我还在学习阶段，还没建模。我现在还没搞懂极性表面如何矫正的问题

小白装大神

wernea 发表于 2015-9-2 10:35
没有，我还在学习阶段，还没建模。我现在还没搞懂极性表面如何矫正的问题

我也一样啊

wernea

小白装大神 发表于 2015-9-2 15:17
我也一样啊

加油

jiewei

wernea 发表于 2015-8-31 20:37
请问表面弛豫的时候固定底层原子，然后表面的原子是允许改变xyz，还是固定xy，只允许改变z。

xyz

wernea

jiewei 发表于 2015-9-4 09:54
xyz

好的，谢谢!

小白装大神

卡开发发 发表于 2015-8-31 10:14
截断能在计算Bulk的体系就能够确定了，k可以按照与Bulk相同的倒空间尺寸来取，之后直接取个真空大一点（ ...

您好！再问你一个初级的问题，原子层数测试确定之后是将relax层数取为测试得到的层数，底部再加固定原子层，还是在测试的这些原子层的底部进行固定去确定relax的原子层数？不知道表达清楚没，谢谢你！