关于极性面钝化添加赝氢问题！
记不清楚在那个话题里面见到有人问如何添加赝氢原子来饱和极性面或者消除内偶极矩。当时斑竹表示出很大的兴趣，让提供一些这方面文献资源。很遗憾，当时没有高手出现解答这个问题，虽然我不是什么高手，但是这个问题我有些心得。当时由于刚做父亲（我们家闺女今天5个月零12天），再加上我的CASTEP计算出点问题，一直没有好好摸索这个问题。
今天有空，就当作送给大伙的新年礼物！
在二六族、三五族纤锌矿结构半导体中，由于沿c轴方向缺少镜面对称，也就是说在c轴方向上是一层A原子一层B原子堆积而成，这些正负电荷的累积导致沿c轴方向出现一个偶极矩。详见PHYSICAL REVIEW B 67, 035403 （2003）Density-functional study of the structure and stability of ZnO surfaces。块体裁成表面时候会出现多种情况（例如电荷转移、表面原子再构、吸附杂质原子等）消除偶极矩。
于是在第一原理模拟这些极性表面时候要考虑如何消除这个偶极矩，使用第一原理程序包最好最直观的办法就是添加赝氢。
VASP里直接提供了很多赝氢如H0.5，H0.75，H0.25，H1.25，等。也就是说这些氢原子是为了某种目的而人为制造的，它的电荷不再是1，而分别为0.5，0.25，0.75，1.25等。
返回小木虫查看更多
今日热帖
京公网安备 11010802022153号
未完待续，老婆叫我提热水兼洗尿布！回头更新！
问一下LZ， castep如何加分数赝氢原子？
我想看表面重构，加真空层， 底表面最下一层必须钝化吗？
谢谢LZ
回复楼上，并继续！
目前我见到的有关赝氢钝化的计算，都是在极性表面上的。如果是模拟极性表面，我觉得需要钝化！
以ZnO为例， Import ZnO，然后 Build----Cleave Surface， cleave plane 填写0 0 1， Top 为0，Depth 3.0，即总共为六个双层。
继续点击Build---Crystals---Build Vacuum Slab，输入真空层厚度，大约为1.0nm。确定即可。
待Supercell构建成功后，点击Add atoms， Element 选择H，Oxidation state选择吸附原子，输入想放置原子的具体位置。然后点中添加的原子，修改Property中的FormalCharge。
赝氢在VASP中的解释就是artificial-atom即人造原子，也就是为了某种需要而认为假设出来的原子，只不过它的电荷数很接近于H，就如同对H原子作了一个简单修改而成，所以叫赝氢！
欢迎讨论阿！
[ Last edited by amynliou on 2009-1-3 at 11:13 ]，
在添加VAccum构建Supercell的时候，我们也可以直接Cleave face 选项框里最底下一栏的 Cap bands on （）faces with （），两个括号是待选项，第一个可以选择Neither，Both， Bottom， Top；第二个可以选择覆盖表面用的元素。
但是这里添加的H的 Formal Charge都是，双击，把它们修改为需要的分数电荷！
[ Last edited by amynliou on 2009-1-3 at 11:09 ]